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Abstract. We introduce and analyze a system of relativistic fermions in a space- 
time continuum, which interact via an action principle as previously considered in 
a discrete space-time. The model is defined by specifying the vacuum as a sum of 
Dirac seas corresponding to several generations of elementary particles. The only 
free parameters entering the model are the fermion masses. We find dynamical field 
equations if and only if the number of generations equals three. In this case, the 
dynamics is described by a massive axial potential coupled to the Dirac spinors. The 
coupling constant and the rest mass of the axial field depend on the regularization; 
for a given regularization method they can be computed as functions of the fermion 
masses. The bosonic mass term arises as a consequence of a symmetry breaking 
effect, giving an alternative to the Higgs mechanism. In addition to the standard 
loop corrections of quantum field theory, we find new types of correction terms to the 
field equations which violate causality. These non-causal corrections are too small 
for giving obvious contradictions to physical observations, but they might open the 
possibility to test the approach in future experiments. 
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1. INTRODUCTION 

In [H] it was proposed to formulate physics based on a new action principle in space- 
time. On the fundamental level, this action principle is defined in so-called discrete 
space-time for a finite collection of projectors in an indefinite inner product space (see 
also [IZ!)- An effect of spontaneous symmetry breaking |16J leads to the emergence of 
a discrete causal structure (see [7] for an explanation in simple examples), which for 
many space-time points and many particles should go over to the usual causal structure 
of Minkowski space (for the connection between discrete and continuum space-times 
we also refer to \24\ \W[ [T7] and the survey article [SO])- Furthermore, on a more 
phenomenological level, it is shown in [14^ Chapters 4-8] that the action can also be 
analyzed in Minkowski space in the so-called continuum limit, where the interaction is 
described effectively by classical gauge fields coupled to second-quantized Dirac fields. 
Finally, in [21j the existence of minimizers of our action principle is proved in various 
situations, both in discrete and in continuum space-times. 

Apart from deriving the general formalism of the continuum limit, in Chap- 
ters 4-8] it is shown that for a suitable system involving 24 Dirac seas, the resulting 
effective gauge group as well as the coupling of the effective gauge fields to the Dirac 
fields have striking similarities to the standard model. However, the detailed form of 
the effective interaction so far has not been worked out. 

This work is the first of a series of papers devoted to the detailed analysis of our 
action principle in the continuum limit and to the derivation of the effective field 
equations. In order to make the presentation as clear and easily accessible as possible, 
our procedure is to begin with small systems, which are composed of only a few Dirac 
seas, and then to gradually build up larger and more complicated systems. In the 
present paper, we consider a system of several Dirac seas (corresponding to several 
"generations" of elementary particles) in the simplest possible configuration referred 
to as a single sector. The only free parameters entering the model are the masses of 
the Dirac particles of each generation. However, we do not specify the form of the 
interaction, which is completely determined by our action principle. Also, we do not 
put in coupling constants nor the masses of gauge bosons. The analysis of the model 
in the continuum limit reveals that we get dynamical field equations if and only if the 
number of generations equals three. In this case, the dynamics can be described by 
a massive axial potential Ag, coupled to the Dirac equation. The corresponding Dirac 
and field equations (stated for notational simplicity for one Dirac particle) become 

{i^ + 7'4a -m)^ = 0, Coj^-C2Al = 12^ "^tS"^ , 

where j'^ = d^^A^^ — ^A^ is the corresponding axial current. The coupling constant and 
the rest mass of the axial gauge field are described by the constants Co and C2 , which 
for a given regularization method can be computed as functions of the fermion masses. 
The mass term of the gauge field arises as a consequence of a symmetry breaking effect, 
giving an alternative to the Higgs mechanism. The field equations involve surprising 
corrections which challenge the standard model of elementary particle physics: First, 
the field equations involve additional convolution terms of the form 

-/[0]*Ja' + 6/[2]*^^ 

where /jpj are explicit Lorentz invariant distributions. These convolution terms give 
rise to small corrections which violate causality. Moreover, we get new types of higher 
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order corrections to the field equations. We also find additional potentials which are 
non-dynamical in the sense that they vanish away from the sources. 

In order to make the paper self-consistent, we introduce our fermion systems and the 
continuum limit from the basics. However, to avoid an excessive overlap with previous 
work, we present a somewhat different point of view, where instead of considering a 
discrete space-time or a space-time continuum of finite volume, we work exclusively 
in Minkowski space. Furthermore, we always restrict attention to a single sector. For 
clarity, we omit the more technical aspects of the regularization, relying instead on 
results from the corresponding chapters of the book . 

The paper is organized as follows. In Chapter [2] we introduce our action principle 
in a space-time continuum. In Chapters [3H5] we review and adapt the methods for 
analyzing this action principle in the continuum as developed in [13]. More precisely, 
in Chapter [3] we describe the vacuum by a system of regularized Dirac seas. We list 
all the assumptions on the vacuum state, either motivating them or explaining how 
they can be justified. In Chapter S] we construct more general fermion configurations 
in Minkowski space by modifying and perturbing the vacuum state, also introducing 
particles and gauge fields. We also outline the mathematical methods for analyzing 
the unregularized fermionic projector with interaction. In Chapter [5l we explain how 
interacting systems are to be regularized, and how to treat the regularization in an 
effective way. This leads us to the formalism of the continuum limit, which allows us 
to analyze our action principle in the continuum, taking into account the unknown 
regularization details by a finite number of free parameters. In the following Chap- 
ters [HHlQ] the continuum limit of our action principle is worked out in detail; this is 
the main part of the paper where we present our new results. Chapter [6] is devoted 
to the leading singularities of the Euler-Lagrange equations on the light cone, where 
the vacuum contributions ( §6.ip are modified by phases coming from the chiral gauge 
potentials ( §6.2p . The next lower orders of singularities are analyzed in Chapter [71 
Then the currents of the gauge fields come into play, and we also get a mass term 
corresponding to the axial gauge field ( §7.ip . Furthermore, we find a corresponding 
contribution of the Dirac current ( ^7.2p . A priori, the different current terms are not 
comparable, because the gauge currents have logarithmic poles on the light cone ( ^7.3p . 
But provided that the number of generations is at least three, these logarithmic poles 
can be compensated by a local axial transformation, as is developed in ^7.41 - ^7.61 

Chapter [8] is devoted to the derivation and analysis of the field equations. In ^S.ll we 
show that the Euler-Lagrange equations corresponding to our action principle give rise 
to relations between the Dirac and gauge currents. If the number of generations equals 
three, we thus obtain field equations for the axial gauge potential (see Theorem 18. 2p . 
These field equations involve non-causal correction terms, which are analyzed and 
discussed in ^8.21 and ^8.31 In ^8.41 we explain how the standard loop corrections of 
quantum field theory (QFT) appear in our framework, and how loop corrections of 
the non-causal terms could be obtained. In ^8.51 we get a connection to the Higgs 
mechanism and explain why in our approach no Higgs particle appears. We finally 
compute the coupling constant and the rest mass of the axial field for a few simple 
regularizations ( §8.6p . 

In Chapter [9] we analyze and discuss further potentials and fields, including scalar 
and pseudoscalar potentials, bilinear potentials, as well as the gravitational field and 
a conformal axial field. In Chapter [10] we consider nonlocal potentials, which can be 
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used to satisfy the Euler-Lagrange equations to higher order in an expansion near the 
origin (see Theorem 110. Sp . 

In order not to interrupt the explanations in the main chapters by longer calcu- 
lations, the more technical parts are worked out in the appendices. Appendix 
supplements the estimates needed for the derivation of the Euler-Lagrange equations 
in the continuum limit in ^5.2[ All the calculations in the formalism of the contin- 
uum limit as needed in Chapters [6H9] are combined in Appendix [Bl which also reviews 
the general method as developed in [14^ Appendix G]. All the formulas given in this 
appendix have been obtained with the help of computer algebra. In Appendix [Cl the 
local axial transformation is worked out non-perturbatively, thus putting the analysis 
in ^7.6l on a fully convincing basis. In Appendix[D]we compute and analyze the smooth 
contributions to the fermionic projector as needed in ^8.11 this is done by modifying 
a resummation technique first introduced in [12] . Finally, in Appendix |E] we outline 
how our constructions and results can be extended to the setting where the Dirac seas 
involve weight factors, as was proposed in ^18j and [23j . 

2. An Action Principle for Fermion Systems in Minkowski Space 

In relativistic quantum mechanics, a fermionic particle is described by a Dirac wave 
function ^ in Minkowski space (Af, (.,.)). In order to describe a many-particle system, 
we consider an operator P on the Dirac wave functions and interpret the vectors in the 
image of P as the occupied fermionic states of the system (for a discussion of the Pauli 
exclusion principle and the connection to the fermionic Fock space formalism see [14:\ 
Chapter 3 and Appendix A]). We assume that P has an integral representation 

{P^){x) = f P{x,y)^{y)dS (2.1) 

with an integral kernel P{x,y). Moreover, we assume for technical simplicity that 
P{x,y) is continuous in both arguments x and y; then the integral in (j2.ip is clearly 
well-defined if for the domain of definition of P we choose for example the space 
Cg°(M)^ of smooth wave functions with compact support. Moreover, we assume that P 
is symmetric with respect to the Lorentz invariant inner product 

<^'|^.>= / ^>{x)^{x) (fx , (2.2) 
Jm 

where ^ = ^'^'^^ is the usual adjoint spinor (^^ is the complex conjugate spinor). In 
other words, we demand that 

<P^|$> = <^r|p$> for aU ^, ^> G C^{M)'' . (2.3) 

This condition can also be expressed in terms of the kernel by 

P(x, y)* = 7°P(x, y)^-i^ = P{y, x) for all x, y G M , (2.4) 

where the dagger denotes the transposed, complex conjugate matrix. We refer to P 
as the fermionic projector. The vectors in the image of P are referred to as the wave 
functions of our system. We point out that for the moment, these wave functions do 
not need to be solutions of a Dirac equation. 

For any space-time points x and y, we next introduce the closed chain A^y by 



A^y = P{x,y) P{y,x) 



(2.5) 
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It is a 4 X 4-matrix which can be considered as a hnear operator on the wave functions 
at X. For any such hnear operator A we define the spectral weight |j4| by 

4 

\A\=^\\\, (2.6) 
1=1 

where Ai,...,A4 are the eigenvalues of A counted with algebraic multiplicities. For 
any x,y G M we define the Lagrangian C by 

^xy[P] = l^xyl — ^ l^xyP • (2-7) 

Integrating over space-time, we can furthermore introduce the functionals 

S[P] '""^^^"^ // C^y[P]d'xdS 

J J MxM ^2 g-j 

T[P] [[ \A^y\'d'xd^. 

J J MxM 

These expressions are only formal because the integrands need not decay for large x 
or y, and thus integrals may be infinite (similar as in classical field theory, where 
the space-time integral over the Lagrangian diverges without imposing suitable decay 
properties at infinity). The functional S is referred to as the action. Our variational 
principle is to minimize S under the constraint that T is kept fixed. For this mini- 
mization we vary the fermionic projector in the following sense. In order to prevent 
trivial minimizers, the variation should preserve the normalization of the wave func- 
tions. This normalization should be performed with respect to the Lorentz invariant 
inner product (j2.2p . However, we do not want to assume that this inner product is 
finite for the wave functions ^ in the image of P (indeed, for physical wave functions, 
the inner product <^'|^'> is in general infinite because the time integral diverges). 
Our method for avoiding the divergences in (|2.2p and (j2.8p is to consider variations 
which outside a compact set are the identity. 

Definition 2.1. An operator U on the Dirac wave functions is called unitary in a 
compact region if 

(i) <C/^' I C/^'> = I for all compactly supported . 

(ii) The operator V := U — \ has the representation 

{V^){x) = / v{x,y)^{y)d^y 

JM 

with a smooth integral kernel v{x,y) which has compact support, i.e. there is a 
compact set K C M such that 

v{x, y) = unless x £ K and y £ K . 

Thus introducing a variation of the wave functions by the transformation ^ ^ U^, 
all the wave functions are changed only in the compact region K C M, in such a way 
that all inner products in this region, i.e. all the integrals 

/ ^'(x) $(x) d'^x , 
Jk 

remain unchanged. Having introduced a well-defined notion of "varying the fermionic 
projector while respecting the inner product (j2.2p " . we can now specify what we mean 
by a minimizer. 
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Definition 2.2. A fermionic projector P of the form (j2.ip is a minimizer of the 

variational principle 

minimize S for fixed T (2-9) 
if for any operator U which is unitary in a compact region and satisfies the constraint 

[ d'x [ d''y(\A^y[P]f-\A,y[UPU-']f)=0, (2.10) 

the functional S satisfies the inequality 

[ d^x [ d^y(/:,y[UPU-']-C,y[P]) > 0. (2.11) 

We point out that, since U changes the wave functions only inside a compact set 
the integrands in (j2.10p and (j2.1ip clearly vanish if x and y are outside K. However, 
is is not obvious that the integrals over the region x £ K and y € M \ K (and 
similarly x € M \ K and y G M) exist. By writing (j2.10p and (j2.1ip we implicitly 
demand that the integrand in (j2.10p and the negative part of the integrand in (j2.1ip 
should be in L^{M x M,R). 

Before going on, we briefly discuss this action principle and bring it into the context 
of previous work. We first remark that, in contrast to \14:\ \17\ . we here ignore the 
condition that P should be idempotent. This is done merely to simplify the presen- 
tation, anticipating that the idempotence condition will not be of relevance in this 
paper. The action principle (|2.9p was first introduced in a discrete space-time in [T^ 
§3.5]. Apart from the obvious replacement of sums by integrals, the action here differs 
from that in [14^ §3.5] only by an irrelevant multiple of the constraint T- This has 
the advantage that the Lagrangian (j2.7p coincides with the so-called critical case of 
the auxiliary Lagrangian as introduced in ^17j : this is the case relevant in our setting 
of one sector. Note that this Lagrangian is symmetric (see [l7l equation (13)]) and 
non-negative, 

£^y[P] = Cy^[P] and C^y[P]>0. 

Moreover, the action principle (12. 9p can be regarded as an infinite volume limit of the 
variational principle in [211 Chapter 3] (possibly also in the limit where the number of 
particles tends to infinity). In the special case of homogeneous systems, our variational 
principle is closely related to the variational principle in infinite volume as considered 
in [211 Chapter 4]. Working with unitary transformations in a compact region, we 
can make sense of the action principle even in infinite space-time volume without 
assuming homogeneity; this procedure can be seen in analogy to considering variations 
of compact support in the Lagrangian formulation of classical field theory (like a 
variation 6 A € C^{M) of the electromagnetic potential in classical electrodynamics). 

3. Assuming a Vacuum Minimizer 

Apart from the general existence results in [17^ [2T] and the simple examples in [Tj [21] , 
almost nothing is known about the minimizers of our action principle. Therefore, 
before we can do physics, we need to assume the existence of a special minimizer 
which describes a physically meaningful vacuum. In this chapter, we compile our 
assumptions on this vacuum minimizer, and we outline in which sense and to what 
extent these assumptions have been justified in [181 123j . At the end of this chapter, 
we will explain how to work with these assumptions in practice. 
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Taking Dirac's original concept seriously, we want to describe the vacuum by "com- 
pletely filled Dirac seas" corresponding to the masses mi, . . . ,mg of g generations of 
elementary particles (later we will set g = 3, but for the moment it is preferable 
not to specify the number of generations). Thus our first ansatz for the integral 
kernel of the fermionic projector of the vacuum is the Fourier transform of the projec- 
tors 2m^(^ + "^Z?) ™ Dirac states on the lower mass shells, 

P{x, y)=pj^^(}^ + mp) 6{e - m|) e{-k^) e-'^^^-y) . (3.1) 

(Here is the Heaviside function, and k{x — y) is a short notation for the Minkowski 
inner product {k,x — y). The slash denotes contraction with the Dirac matrices, 
thus }^ = kjj^ . We always work in natural units h = c = 1, and for the signature of 

the Minkowski inner product we use the convention (H ).) We always index the 

masses in increasing order, 

mi < 1712 < . . . < nig . (3.2) 

The Fourier integral ()3.ip is well-defined as a distribution. If the vector y — x is 
spacelike or timelike, the integral (j3.ip exists even pointwise. However, if the vector 
y — X is null, the distribution P{x, y) is singular (for details see |14[ §2.5]). In physical 
terms, these singularities occur if y lies on the light cone centered at x. Thus we 
refer to the singularities on the set where (x — y)^ = as the singularities on the 
light cone. As a consequence of these singularities, the pointwise product in (j2.5p is 
ill-defined, and the Lagrangian (12. 7p has no mathematical meaning. In order to resolve 
this problem, one needs to introduce an ultraviolet regularization. In position space, 
this regularization can be viewed as a "smoothing" on a microscopic length scale. 
It seems natural to identify this microscopic length scale with the Planck length ip. 
Likewise, in momentum space the regularization corresponds to a cutoff or decay on 
the scale of the Planck energy Ep = £p^. Clearly, the Planck scale is extremely 
small compared to the length scale ^macro of macroscopic physics, and thus it seems 
reasonable to expand in powers of ^p/^macm- However, such an expansion would 
not be mathematically meaningful, because Taylor series can be performed only in 
continuous variables (but not in a constant, no matter how small). Therefore, it is 
preferable to denote the regularization length by the variable e, which may vary in 
the range < e ^ ^macro (clearly, e can be always be thought of as being of the order 
of the Planck length). We are thus led to a one-parameter family of regularizations. 
We assume that these regularized Dirac sea configurations are all minimizers. We also 
compile all assumptions on the regularization as introduced in jl4| Chapter 4]. 

Assumption 3.1. (regularized Dirac sea vacuum) There is a family (P^)£>o of 
fermionic projectors whose kernels P^{x,y) (as defined by (j2.ip ) have the following 
properties: 

(i) Every P^{x,y) is a minimizer in the sense of Definition 12.21 

(ii) Every P^{x,y) is homogeneous, i.e. it depends only on the variable ^ := y — x. 

(iii) Taking its Fourier transform, 

P'ix. y) = j P'ip) ^-'"^^-'^ > (3.3) 
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P'^ is a distribution with a vector-scalar structure, i.e. 

P^k) = (v'^ik) 7^' + c^%k) 1) fip) (3.4) 

with a vector field v^, a scalar field (p^ and a distribution f^, which are all 
real- valued. 

(iv) If the regularization is removed, goes over to P (as given by (j3.ip ). i.e. 

9 

lim P'{k) = P{k) := ^(^ + m^) 6{k'^ - m}) G(-A:°) 

/3=1 

with convergence in the distributional sense. 

The assumptions so far seem natural and are easy to state. The assumption (iii) can 
even be derived in the Hamiltonian framework by minimizing the physical energy of a 
regularized Dirac sea configuration, describing the electrostatic interaction in the mean 
field approximation [291 ISO]- To understand the following assumptions, one should 
notice that the singularities of P{x, y) on the light cone arise because its Fourier 
transform P{k) is supported on the mass shells k"^ = m'^, which are hypersurfaces 
being asymptotic to the mass cone k'^ = (for details see [HI §4.2]). Thus in order to 
control the behavior of P^ near the light cone, we need to make suitable assumptions 
on P^{ui, k) for uj ^ —\k\ ~ e~^. 

(v) We assume that the distribution P'^ is supported on hypersurfaces described 
by graphs, i.e. the distribution in (|3.4p should be of the form 

g 

ficu, k) = J2^{^ + \k\+ a^{k)) . (3.5) 

/3=1 

These hypersurfaces should be asymptotic to the mass cone in the sense that 
afs{k) ~ e if I A; I ~ e^^ . 

Except for these singularities, P^{k) is so regular that the singularities as e \ 
of P^{x, y) on the light cone are completely described by the behavior of P^{k) 
on the hypersurfaces (|3.5p . up to corrections of higher order in e. We refer to 
this assumption as the restriction to surfaces states. 

(vi) On the hypersurfaces (13. 5p and for \k\ ~ the vector field in (]3.4p should 
be parallel to /c, up to a small error term. More precisely, decomposing as 

= s^{k)k + w^{k) 

with a scalar function s"^, the vector field w'^ should be bounded by 

^shear 

where Eghear < 1 • (3.6) 

Referring to the effect of this assumption in position space, we say that the 
vector component is null on the light cone. 

(vii) The functions in (j3.4p either vanish, (j)^{k) = = v^{k), or else (j)^{k) > and 
the vector field is time-like and past-directed. Furthermore, 

v%pf = <P%pf . 

For a discussion of the assumptions (v) and (vi) we refer to [141 Chapter 4]. The 
condition (vii) requires a brief explanation. This assumption is clearly satisfied without 
regularization ()3.ip (in which case we choose v{p) = p/{2uj) and (j) a positive function 
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which on the mass shells takes the values ma/{2uj)). A closely related condition was 
first proposed in [TH Chapter 4] as the assumption of half-occupied surface states. 
This condition was motivated by the wish to realize the Dirac sea configurations with 
as few occupied states as possible, noting that the condition (vii) implies that the 
matrix P^{k) has rank at most two. Furthermore, the condition (vii) implies that the 
image of the matrix P^{k) is negative definite with respect to the inner product 
Prom the mathematical point of view, this definiteness ensures that can be regarded 
as a negative definite measure on an indefinite inner product space as introduced in |2H 
Chapter 4], and this assumption is crucial for the general compactness result |2H 
Theorem 4.2]. Thus the physical intuition and the mathematical requirements fit 
together. Moreover, in the case when P'^{k) does not vanish, we can choose a suitably 
normalized orthogonal basis (^'p,i, ^p,2) of the image of P^{k) such that {27r)^ P^{k) = 
— ^fc^i^'fc^i — ^A:,2^A;,2- Substituting this representation into the Fourier integral ()3.3p 
and using (13. 5p . we obtain 




where ^^^^(3.) = ^p,a e for p = {—\k\ — ais{k),k). This representation is helpful 
because it shows that the regularized fermionic projector of the vacuum is composed 
of negative-energy wave functions; the index a can be thought of as describing the two 
spin orientations. 

We next outline the approach taken to justify the above assumptions. In |18j a class 
of regularizations is constructed for which the action remains finite when the regular- 
ization is removed (more precisely, this is done by proving that the constructed regular- 
izations satisfy the so-called assumption of a distributional A^P-product). These regu- 
larizations are spherically symmetric, but they break the Lorentz symmetry. However, 
after suitably removing the regularization, we obtain a well-defined Lorentz invariant 
action principle. This Lorentz invariant action principle is analyzed in [23J, and it is 
shown that for certain values of the masses and the so-called weight factors (which for 
simplicity we do not consider in the main text of this paper; however, see Appendix[E|). 
the Dirac sea configuration (j3.ip is indeed a minimizer, in a sense made precise using 
the notion of state stability. Following these results, "good candidates" for satisfying 
the above assumptions are obtained by regularizing the state stable Dirac sea configu- 
rations from [23] according to to the regularization scheme in [18j . The remaining task 
for giving a rigorous justification of Assumption 13.11 is to use the freedom in choosing 
the regularization such as to obtain a minimizer in the sense of Definition 12.21 This 
task seems difficult and has not yet been accomplished. In [211 Theorem 4.2] the ex- 
istence of minimizers is proved within the class of homogeneous fermionic projectors; 
but this is considerably weaker than being a minimizer in the sense of Definition 12.21 
In technical terms, the main difficulty is to quantify the influence of the spherically 
symmetric regularization on the action, even taking into account contributions which 
remain finite when the regularization is removed. Despite this difficult and technically 
challenging open problem, it is fair to say that the results of [18^, '23] show that Dirac 
sea configurations tend to make our action small, thus explaining why Assumption 13. II 
is a reasonable starting point for the continuum analysis. 

We finally explain how to work with the above assumptions in practice. Ideally, 
the fields v^, (p^ and the distribution in (j3.4p could be determined by minimizing 
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our action (j2.9p . thus giving detailed information on . Such a minimization process 
is indeed possible (see |21[ Theorem 4.2] for a general existence result and [24J for a 
lattice formulation), but so far has not been analyzed in sufficient depth. Thus for the 
time being, there is a lot of freedom to choose the functions in (j3.4p . Our program is 
not to make a specific choice but to consider instead general functions f^, t/)^ and 
Our subsequent analysis will clearly depend on the choice of these functions, and our 
task is to look for conclusions which are robust to regularization details. This so-called 
method of variable regularization (which is worked out in detail in [HI §4.1]) leads to 
the formalism of the continuum limit which will be explained in Chapter [5] below. 



4. Introducing an Interaction 

Our next goal is to generalize the regularized fermionic projector of the previous 
chapter so as to include an interaction. Postponing the treatment of the regularization 
to Chapter [5l we shall now extend the definition of the fermionic projector of the 
vacuum (j3.ip to the case with interaction. We outline the methods developed in [Til 
HailS]; see also [Ml Chapter 2]. 

4.1. A Dirac Equation for the Fermionic Projector. First, it is useful to re- 
cover (j3.ip as a solution of a Dirac equation: Replacing the ordinary sum in (j3.ip by 
a direct sum, we introduce the so-called auxiliary fermionic projector p^"^ by 

^'^"^(^, y) = / (04 + "^^) ^i^' - ^) ^"'^''"'^ (4-1) 

(thus P^^^{x,y) is represented by a 4(7 x 4g-matnx). It is a solution of the free Dirac 
equation 

{i^-mY)P''''''{x,y) = , (4.2) 

where the mass matrix Y is composed of the rest masses corresponding to the three 
generations, 

9 

my = 0m^ (4.3) 

/9=1 

(here m > is an arbitrary mass parameter which makes Y dimensionless and will 
be useful for expansions in the mass parameter; see also [141 §2.3]). The fermionic 
projector of the vacuum is obtained from p^^^ by taking the partial trace over the 
generations defined by 

9 

P= • (4-4) 

a,/3=l 

The obvious idea for introducing an interaction is to replace the free Dirac equa- 
tion (|4.2p by a Dirac equation with interaction, 

(i^ + S-my)P^''^(x,y) = 0, (4.5) 

where !B is a general perturbation operator, and to introduce the fermionic projector 
again by taking the partial trace ()4.4p . In order to ensure that the resulting fermionic 
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projector is again symmetric ()2.3p . we generalize the inner product (j2.2p to the wave 
functions of the auxihary Dirac equation by setting 




and demand that the auxihary fermionic projector is symmetric with respect to this 
new inner product, 

<P"""^aux|$aux> = <^aux|P"""^aux> for aU ^aux, ^-aux G {Mf' . (4.7) 

In order to obtain a coherent framework, we shah always assume that the auxiliary 
Dirac operator is also symmetric, 

<(i^+S-my)^aux|^aux> = <*aux|(i^ + S-my)«>aux> • (4.8) 

This equation gives a condition for the operator !B describing the interaction. Apart 
from this condition and suitable regularity and decay assumptions, the operator S can 
be chosen arbitrarily; in particular, it can be time dependent. In typical applications, 
B is a multiplication or differential operator composed of bosonic potentials and fields. 
The choices of 23 relevant for this work will be discussed in §4.51 below. 

4.2. The Interacting Dirac Sea. Clearly, the Dirac equation (j4.5p has many dif- 
ferent solutions, and thus in order to determine p^"^^^ we need to specify of which 
one-particle states p^"^ should be composed. In the vacuum (j4.ip . this can be done 
by taking all the negative-energy solutions, i.e. all states on the lower mass shells 
{A;^ = m?p,k^ < 0}. Unfortunately, the concept of negative energy does not carry 
over to the situation of a time-dependent interaction ()4.5p . because in this case the 
energy of the Dirac wave functions is not conserved; this is the so-called external field 
problem (see [14^ §2.1]). The clue for resolving this problem is the observation that the 
negative-energy states in (14. ip can be characterized alternatively using the causality of 
the Dirac Green's functions in a specific way. This causal approach generalizes to the 
situation (|4.5p and allows to extend the concept of the Dirac sea to the time-dependent 
setting. It gives rise to a unique definition of a fermionic projector P^'^^ in terms of 
a power series in 23. More precisely, the so-called causal perturbation expansion ex- 
presses P^^^ as sums of operator products 

OO "max (A:) 

psea^^ CaCi,,23C2,a23 ••• SCfc+l,„, (4.9) 

A;=0 a=0 

where the factors Ci^a are the Green's functions or fundamental solutions of the free 
Dirac equation (14. 2p . and the Ca are combinatorial factors (for details see [TTl and [141 
§2.2, §2.3]; for a more recent account on idempotence and unitarity questions see [22j). 
In the language of Feynman diagrams, each summand in (14. 9p is a tree diagram. 
These tree diagrams are all finite, provided that 23 satisfies suitable regularity and 
decay assumptions at infinity (see [14^ Lemma 2.2.2.]). 

4.3. Introducing Particles and Anti-Particles. The fermionic projector p^''^ is 
interpreted as a generalization of completely filled Dirac seas to the interacting situ- 
ation (j4.5p . In order to bring particles and anti-particles into the system, we add the 
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projectors on states ^'i, . . . , which are not contained in the image of the opera- 
tor P^^^ (the particle states) and subtract the projectors on states $i, . . . , which 
are in the image of P^"^^ (the anti-particle states), 

P^^(x,y) = P^^^(x,2/) - _^^,(x)^fc(y) + _ ^ cl.,(x)cl>;(y) . (4.10) 

k=l 1=1 

Finahy, the fermionic projector is again obtained by taking the partial trace over the 
generations (14. 4p . 

The wave functions in ()4.10p are normalized such that they are orthonormal with 
respect to the usual integral over the probability density, 

/ {Wl^^^i){t,x)dx = 6ki= [ (^kl°mt,x)dx. (4.11) 

The factors ±^ in (HTOj) are needed for the proper normalization of the fermionic 
states with respect to the inner product (j2.2p . The fact that this last inner product is 
indefinite has the effect that the normalization factors in (|4.10p have the opposite sign 
from what one would have expected in Hilbert spaces. These normalization conditions 
are derived in |141 §2.6] by considering the system in finite 3-volume and taking the 
infinite volume limit. For the dependence of the normalization constants on the global 
geometry of space-time see [25] . 

4.4. The Light-Cone Expansion and Resummation. We now outline the meth- 
ods for analyzing the fermionic projector in position space (for details see [12\ I13|). 
The following notion is very useful for describing the structure of the singularities on 
the light cone. 

Definition 4.1. A distribution A{x, y) on M x M is of the order 0((y — x)'^^), p £ Z, 
if the product 

{y-x)-^P A{x,y) 

is a regular distribution (i.e. a locally integrable function). It has the light-cone 
expansion 

oo 

A{x,y) = ^ A^^{x,y) 

J =90 

with £ Wj if the distributions A^^\x,y) are of the order 0{{y — x^^) and if A is 
approximated by the partial sums in the sense that for all p > g, 

p 

A{x,y) - ^ A^^^{x,y) is of the order 0{{y - xf'f^'^) . 

j=90 

Thus the light-cone expansion is an expansion in the orders of the singularity on the 
light cone. As the main difference to a Taylor expansion, the expansion parameter (y — 
x)^ vanishes for any fixed x on an unbounded set, namely the whole light cone centered 
at X. In this sense, the light-cone expansion is a nonlocal expansion. 

For a convenient formulation of the light-cone expansion of the fermionic projector, it 
is helpful to work with a generating function, i.e. a power series in a real parameter a > 
whose coefficients are functions in (y — x)^ which are of increasing order on the light 
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cone. The first ansatz for such a generating function is the Fourier transform Ta{x,y) 
of the lower mass sheh with k"^ = a, 

Ta{x, y) = J^, Kk' - a) G(-fcO) e-^'^^-y^ . (4.12) 

Carrying out the Fourier integral and expanding the resulting Bessel functions, one 
obtains 



+ ^ t (lo. we\ + c, + e(,^) eie)) ^ ^ 



(4.13) 



where we again used the abbreviation ^ = y — x, and e denotes the sign function (i.e. 
e{x) = 1 if X > and e{x) = — 1 otherwise). The real coefficients cj are given explicitly 
in [m §2.5]. Unfortunately, due to the factor log|a^^|, the expression (I4.13j) is not a 
power series in a. In order to bypass this problem, we simply remove the logarithms 
in a by subtracting suitable counter terms, 

T-(x,,) := Ux,y) - ^ log \a\ g . (4.14) 

The resulting distribution Ta^^ is a power series in a, and it is indeed the right choice 
for our generating function. We denote its coefficients by 



da 



(n = 0,1,2,...) (4.15) 



a=0 

and also introduce T^^^^ via the distributional equation 

AT(0)(x,y) = ^ (y-x),r(-i)(x,y) . (4.16) 

Combining Fourier techniques with methods of hyperbolic partial differential equa- 
tions, one can perform the light-cone expansion of each summand of the perturbation 
series (|4.9p . After suitably rearranging all the resulting contributions, one can partially 
carry out the infinite sums. This so-called resummation gives rise to an expansion of 
the interacting fermionic projector of the form 

oo 

P^'^^{x,y) = m^'= (phase-inserted nested line integrals) x T^^\x,y) 

n=— 1 k 

+ P^%x,y) + P'^''{x,y) . (4.17) 

Here the n-summands describe the different orders of the singularities on the light 
cone, whereas the A;-sum describes all contributions to a given order on the light 
cone. The phase-inserted nested line integrals involve 23 and its partial derivatives, 
possibly sandwiched between time-ordered exponentials of chiral potentials. Since 
these nested line integrals are smooth functions in x and y, the series in (14.17P is a 
light-cone expansion in the sense of Definition 14. !( provided that the fc-sum is finite 
for every n. This is indeed the case if S is composed of scalar, pseudoscalar and chiral 
potentials |13j , whereas for a more general perturbation operator S this condition still 
needs to be verified. This expansion is causal in the sense that it depends on S and 
its partial derivatives only along the line segment xy. The contributions P^'^ and P^^, 
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on the other hand, are not causal but depend instead on the global behavior of 23 in 
space-time. They can be written as a series of functions which are all smooth in x 
and y. Their different internal structure gives rise to the names non-causal low energy 
contribution and non-causal high energy contribution, respectively. 

For an introduction to the light-cone expansion and the needed mathematical meth- 
ods we refer to [12] and [13] or to the exposition in [14^ §2.5]. The formulas of the 
light-cone expansion needed in this work are compiled in Appendix [Bl 

4.5. Clarifying Remarks. The above constructions require a few explanations. We 
first point out that, although we are working with one-particle wave functions, the 
ansatz for the fermionic projector (|4.10p describes a many-particle quantum state. 
In order to get a connection to the Fock space formalism, one can take the wedge 
product of the wave functions and <1>; to obtain a vector in the fermionic Fock 
space (for details see [HI Appendix A]). We conclude that (|4.10p describes second- 
quantized fermions. For the description of entangled states see [T9] . 

One should keep in mind that at this stage, the form of the potential 23 has not 
been specified; it can be an arbitrary operator. Indeed, we regard the operator 23 only 
as a device for modifying or perturbing the fermionic projector. We do not want to 
preassume which of these perturbations are physically relevant; instead, we want to 
select the relevant perturbations purely on the basis of whether they are admissible for 
minimizers of our action principle (12. 9p . In other words, our action principle should 
decide how the physical interaction looks like, even quantitatively in the sense that 
our action principle should determine the corresponding field equations. Following this 
concept, we should choose 23 as general as possible, even allowing for potentials which 
are usually not considered in physics. We now give a brief overview over the potentials 
which will be of relevance in the present work. The simplest example is to choose an 
electromagnetic potential. 



More generally, one can choose chiral potentials, which may be non-diagonal in the 
generations. 



usual pseudoscalar matrix). To describe a gravitational field, one needs to choose for 23 
a differential operator of first order; more precisely. 



where P is the Dirac operator in the presence of a gravitational field. 

The above choices of 23 are of course motivated by physical fields observed in nature. 
However, we point out that we do not assume any field equations. Thus the electro- 
magnetic potential in (I4.18P does not need to satisfy Maxwell's equations, in (I4.19P 
we do not assume the Yang- Mills equations, and in ()4.20p the Einstein equations are 



For convenience we shall always omit the coupling constant e in the Dirac equation. Our con- 
vention is obtained from the usual choice S = by the transformation A — >■ e~^A. The coupling 
constant clearly reappears in the Maxwell equations, which we write in natural units and with the 
Heaviside-Lorentz convention as djkA^ — DAfc = e^$7fe4'. As usual, the fine structure constant is 
given by a = /{Atv). 



B = 4. 



(4.18) 




'B = V-i^ 



(4.20) 
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not imposed. This is because, as already pointed out above, our goal is to derive the 
classical field equations from our action principle (j2.9p . 

Apart from the above choices of "B motivated from physics, one can also choose other 
physically less obvious operators, like for example scalar or pseudoscalar potentials, 

■B = <S> + i-f^E (4.21) 

with <1> = <1>^, H = and a,/3 = 1, . . . ,g. Furthermore, one can consider a scalar 
differential operator, 

■B = i^^dj , 

or a higher order differential operator. More specifically, we will find a pseudoscalar 
differential potential useful. 

It is worth noting that one does not need to restrict attention to differential operators. 
Indeed, "B can also be an integral operator, in which case we talk of nonlocal potentials. 
Clearly, one can also take linear combinations of all the above operators S. 

Next, it is worth noting that for the moment, we consider 23 as a-priori given, and 
thus at this stage, our system consists of Dirac particles moving in an external field. 
However, our action principle (|2.9|) will give relations between the potentials contained 
in S and the Dirac wave functions in (|4.1U|) . and thus these potentials will be influenced 
by the Dirac wave functions. This leads to a mutual coupling of the potentials to the 
the Dirac wave functions, giving rise to a fully interacting system. We also point out 
that the potentials and fields contained in 23 should be regarded as classical. Indeed, 
in this paper we will always work with classical bosonic fields. However, as is worked 
out in |19j . the framework of the fermionic projector also allows for the description of 
second-quantized bosonic fields. 

4.6. Relation to Other Approaches. Having outlined our approach, we can now 
give a short review of related works. In order to get a connection to our description of 
the Dirac sea in ^4.2| we begin with the construction of quantum fields in an external 
field. Historically, this problem was first analyzed in the Fock space formalism. Klaus 
and Scharf [35\ [36] considered the Fock representation of the electron-positron field in 
the presence of a static external field. They noticed that the Hamiltonian needs to be 
regularized by subtracting suitable counter terms which depend on the external field. 
They also noticed that the electron-positron field operators in the external field form a 
Fock representation on the standard Fock space of free fields only if the external field 
satisfies a certain regularity condition. This regularity condition is quite restrictive 
and excludes many cases of physical interest (in particular a magnetic field |38j and 
a Coulomb potential |34j). These results show that different external fields in general 
give rise to nonequivalent Fock representations of the electron-positron field operators. 
More recently, in |281 [27] the vacuum state was constructed for a system of Dirac 
particles with electrostatic interaction in the Bogoliubov-Dirac-Fock approximation. 
The conclusion of this analysis is that for mathematical consistency, one must take 
into account all the states forming the Dirac sea. Furthermore, the interaction mixes 
the states in such a way that it becomes impossible to distinguish between the particle 
states and the states of the Dirac sea. 

In the time-dependent setting, Fierz and Scharf [lOj proposed that the Fock rep- 
resentation should be adapted to the external field as measured by a local observer. 
Then the Fock representation becomes time and observer dependent. This implies that 
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the distinction between particles and anti-particles no longer has an invariant meaning, 
but it depends on the choice of an observer. In this formulation, the usual particle 
interpretation of quantum states only makes sense for the in- and outgoing scattering 
states, but it has no invariant meaning for intermediate times. This picture has been 
confirmed in a rigorous self-consistent manner in [31]. For a related approach which 
allows for the construction of quantum fields in the presence of an external magnetic 
field see [6]. In all the above approaches, the Dirac sea leads to divergences, which 
must be treated by an ultraviolet regularization and suitable counter terms. 

As an alternative to working with Fock spaces, one can use the so-called point 
splitting renormalization method, which is particularly useful for renormalizing the 
expectation value of the energy-momentum tensor [3j. The idea is to replace a function 
of one variable T(x) by a two-point distribution T{x,y), and to take the limit y — )• x 
after subtracting suitable singular distributions which take the role of counter terms. 
Analyzing the singular structure of the counter terms leads to the so-called Hadamard 
condition (see for example [26j). Reformulating the Hadamard condition for the two- 
point function as a local spectral condition for the wave front set [M] turns out to 
be very useful for the axiomatic formulation of free quantum fields in curved space- 
time. As in the Fock space formalism, in the point splitting approach the particle 
interpretation depends on the observer. This is reflected mathematically by the fact 
that the Hadamard condition specifies the two-point distribution only up to smooth 
contributions, thus leaving the smooth particle wave functions undetermined. For a 
good introduction to free quantum fields in curved space-time we refer to the recent 
book [1]. 

As mentioned at the beginning of §4.51 in our approach the connection to the Fock 
space formalism is obtained by choosing a basis of the image of the fermionic projector 
and taking the wedge product of the basis vectors (for details see [HI Appendix A] 
or [19]). If in this construction the states of the Dirac sea are taken into account, 
we get precisely the framework in [10]. The connection to the Hadamard condition 
is even closer. Indeed, considering the light-cone expansion locally for y near x, the 
summands in (|4.17p coincide precisely with the singular distributions in the Hadamard 
construction. Since the non-causal contributions P^^ and P^'^ are smooth functions, 
we conclude that the integral kernel of the fermionic projector satisfies the Hadamard 
condition, provided that the perturbation expansions for P'^^ and P'^ converge (a subtle 
technical problem which we do not want to enter here). Thus in a given external field, 
P^^^{x, y) can be interpreted as the two-point function, and using the methods of [41' T] 
one could construct the corresponding free QFT. 

A major difference of our approach is that our framework allows for the description 
of an interacting theory, where the coupling of the fermions to bosonic fields and the 
back-reaction of the bosonic fields to the fermions is taken into account. In this setting, 
the interaction is described by our action principle (|2.8p . The mathematical framework 
is no longer equivalent to standard QFT. In particular, P{x, y) cannot be interpreted 
as the two-point function of a corresponding QFT, simply because the notions of QFT 
can no longer be used. But we still get a connection to the Feynman diagrams of QFT 
(as will be explained in ^8.41 below). 

Another major difference of our approach is that the distribution P^^^ as defined 
by the causal perturbation expansion (|4.9p distinguishes a unique state which can be 
interpreted as the fermionic vacuum state where all Dirac seas are completely filled. 
Thus working relative to this distinguished state, there is a unique observer independent 
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particle interpretation, even at intermediate times (see |1H Section 5] for a discussion 
of this point). At first sight, this distinguished particle interpretation might seem of 
purely academic interest, because P^^^ is defined globally in space-time and is thus 
not accessible to a local observer. However, our action principle (j2.8p does have access 
to quantities defined globally in space-time, and in this way the distinguished particle 
interpretation enters the physical equations. More precisely, P^^^ drops out of the 
Euler-Lagrange equations corresponding to our action principle, up to terms which are 
well-defined and explicitly computable, even determining the smooth contributions. In 
this way, the arbitrariness of working modulo smooth contributions (in the Hadamard 
condition) or modulo regular counter terms (in the Fock space formalism) is removed. 
The corresponding smooth contributions to the physical equations will be analyzed 
in ^8.11 and Appendix [Dj They are nonlocal and violate causality, as will be explained 
in go 

A frequently asked question is how our approach relates to Connes' noncommutative 
geometry [5j. In particular, can our approach be thought of as a Lorentzian version of 
noncommutative geometry? Clearly, both approaches have in common that the Dirac 
operator plays a central role. Moreover, the light-cone expansion is the Lorentzian 
analog of local expansions of the resolvent near the diagonal. A major difference is 
that instead of considering the whole spectrum of the Dirac operator, we only consider 
the eigenspaces corresponding to the masses of the Dirac particles of our system. 
Furthermore, we only take "half the eigenspaces" by constructing Dirac seas and build 
in additional particle and anti-particle states (14.10p . Another major difference concerns 
the mathematical structure of our action principle (j2.8p . Namely, this action cannot 
be thought of as a spectral action, because it is impossible to express it in terms of 
spectral properties of the Dirac operator. This is obvious from the fact that in (12. 7p 
and (12. Sp we perform a nonlinear (and even non-analytic) transformation of the ker- 
nel P{x, y) before integrating over x and y. As a consequence, there is no connection 
to a regularized trace or Hilbert-Schmidt norm of P. The specific form of our action 
principle makes it possible to regard the structures of Minkowski space as emerging 
from a self-organization of the wave functions in discrete space time (see [20]), an 
idea which has no correspondence in noncommutative geometry. On the other hand, 
noncommutative geometry has deep connections to Riemannian geometry, index the- 
ory and number theory. We conclude that despite superficial similarities, the aims, 
ideas and methods of our approach are quite different from those in noncommutative 
geometry. 

5. The Continuum Limit 

In Chapter [3] we described the vacuum by a family of regularized fermionic pro- 
jectors P^ . Our next goal is to use the information on the regularized vacuum to 
also regularize the fermionic projector with interaction. We cannot expect that this 
information will suffice to determine the interacting fermionic projector in all details, 
because it is unknown how the interaction affects the fermionic projector on the mi- 
croscopic scale. But as shown in [14^ Chapter 4 and Appendix D], there is a canonical 
method to regularize the formulas of the light-cone expansion (j4.17p . This method also 
gives a meaning to composite expressions as needed for the analysis of the action prin- 
ciple introduced in Chapter [2l In particular, it allows us to analyze the corresponding 
Euler-Lagrange equations in the continuum, taking into account the unknown regu- 
larization details by a finite number of free parameters. We now outline this method, 
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relying for all technical issues on the detailed analysis in [M]. The method in ^5.21 is 
a major improvement and simplification of the techniques in |14[ Appendix F]. 

5.1. Weak Evaluation on the Light Cone. Our method is based on the physically 
reasonable assumption of macroscopic potentials and wave functions which states that 
both the bosonic potentials in (j4.5p and the fermionic wave functions in ()4.10p vary 
only on the macroscopic scale and are thus almost constant on the Planck scale. Then 
the idea is to try to regularize the perturbation expansion (j4.9p in such a way that 
the interaction modifies the fermionic projector also only on the macroscopic scale. 
As exemplified in jl4| Appendix D] in the perturbation expansion to first order, this 
idea can be realized by demanding that the perturbation expansion should be gauge 
invariant and satisfies a causality condition. Performing the light-cone expansion for 
the thus regularized perturbation expansion and using the form of the regularized 
vacuum minimizers as specified in Assumption 13.11 one obtains a simple regularization 
scheme for the continuum fermionic projector (|4.17p . which we now describe. 

The non-causal contributions P^'^ and P^'^, which are already smooth in x and y, are 
not regularized. Likewise, the smooth nested line-integrals are not regularized. Thus 
we only regularize the singularities of the factors T^") on the light cone, and this is 
done by the replacement rule 

^mPrjJj\ (5.1) 

where the factors T^^\^) are smooth functions defined by Fourier integrals involving 
the functions f^, (j)'^ and in the ansatz ()3.4p . If the partial trace is taken, we 
clarify the handling of the generation index by accents. More precisely, we extend the 
replacement rule ()5.ip to (see also [IH §7.1]) 

9 

mf Y^^^--- ^7/ • • • Yj''-' r(") ^ yy • • • y r^Jj^ , (5.2) 

a,/3.7ii---,7p-i = l ^ 



p factors Y 



and in particulai0 
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S'^ ^ m% and ^ m Y^ T^") mY T^^ . 

Q,/3=l o,/3=l 

(n) 

Fortunately, the rather complicated detailed form of the factors Tj^j will not be needed 
here, because these functions can always be treated symbolically using the following 
simple calculation rules. In computations one may treat the Tj^"'* like complex func- 
tions. However, one must be careful when tensor indices of factors ^ are contracted 
with each other. Naively, this gives a factor which vanishes on the light cone and 
thus changes the singular behavior on the light cone. In order to describe this effect 
correctly, we first write every summand of the light cone expansion ()4.17p such that it 
involves at most one factor ^ (this can always be arranged using the anti-commutation 
relations of the Dirac matrices). We now associate every factor ^ to the corresponding 



^ In contrast to the convention in [14], here we always write out the factors g counting the number 
of generations (in [14], the factor g was absorbed into the factors Tj^"' and T^"')- The shorter notation 
in [14| has the disadvantage that reinserting the factors of g in the end is a potential source of confusion 
and may lead to computational errors. In the convention here, the factors ri"-* without regularization 
always coincide with the distributions (|4.15|l and (|4.16|l . 
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(n) 

factor ^[p] • In simple calculations, this can be indicated by putting brackets around 
the two factors, whereas in the general situation we add an index to the factor ^, giving 
rise to the replacement rule 

The factors ^ which are contracted to other factors ^ are called inner factors. The 
contractions of inner factors can be handled with the so-called contraction rules 



'"■•"""■''^■-.fi-;' +-!„"?) (5.3) 



^i>i:^=-^(nTir''+T}:t'')^ (5-5) 



which are to be complemented by the complex conjugates of these equations. Here 



the factors Zj^^ can be regarded simply as a book-keeping device to ensure the correct 
application of the rule (j5.5p . The factors t|"? have the same scaling behavior as 



(n) 

the T^pj , but their detailed form is somewhat different; we simply treat them as a new 
class of symbol^. In cases where the lower index does not need to be specified we 
write To'^\ After applying the contraction rules, all inner factors ^ have disappeared. 
The remaining so-called outer factors ^ need no special attention and are treated like 
smooth functions. 

Next, to any factor tJ"^ we associate the degree degTj"^ by 

degri"^ = 1 - n. 

The degree is additive in products, whereas the degree of a quotient is defined as the 
difference of the degrees of numerator and denominator. The degree of an expression 
can be thought of as describing the order of its singularity on the light cone, in the sense 
that a larger degree corresponds to a stronger singularity (for example, the contraction 
rule (j5.5p increments n and thus decrements the degree, in agreement with the naive 
observation that the function z = vanishes on the light cone). Using formal Taylor 
expansions, we can expand in the degree. In all our applications, this will give rise to 
terms of the form 



r](x,y) ^^^=^^^= with ?7(x, y) smooth . (5.6) 

rp{ci) rp{Cj) rp{dl) rp{ds) 

Here the quotient of the two monomials is referred to as a simple fraction. 

A simple fraction can be given a quantitative meaning by considering one-dimensional 
integrals along curves which cross the light cone transversely away from the ori- 
gin ^ = 0. This procedure is called weak evaluation on the light cone. For our purpose, 
it suffices to integrate over the time coordinate t = for fixed ^ 7^ 0. Using the 
symmetry under reflections ^ — >• — ^, it suffices to consider the upper light cone t ^ |^|. 



''We remark that the functions T^py will be of no relevance in this paper, because they contribute 



to the EL equations only to degree three and lower; see i)9.4l 
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The resulting integrals will diverge if the regularization is removed. The leading con- 
tribution for small e can be written as 



where L is the degree and Creg, the so-called regularization parameter, is a real-valued 
function of the spatial direction ^/|^| which also depends on the simple fraction and 
on the regularization details (the error of the approximation will be specified below). 
The integer r describes a possible logarithmic divergence; we postpone its discussion 
until when we need it (see ^7.3p . Apart from this logarithmic divergence, the scalings 
in both and e are described by the degree. 

When analyzing a sum of expressions of the form (|5.6p . one must know if the cor- 
responding regularization parameters are related to each other. In this respect, the 
integration-by-parts rules are important, which are described symbolically as follows. 

(n) 

On the factors To we introduce a derivation V by 

Extending this derivation with the Leibniz and quotient rules to simple fractions, the 
integration-by-parts rules states that 



° " ° 1 1 " ° 1=0. (5.8) 

J,{ci) _ _ _ rp{c^l) rp{dl) _ _ _ rp{ds) 

These rules give relations between simple fractions (the name is motivated by the fact 
that when evaluating ()5.8p weakly on the light cone ()5.7p . the rules state that the 
integral over a derivative vanishes). Simple fractions which are not related to each by 
the integration-by-parts rules are called basic fractions. As shown in [14, Appendix E], 
there are no further relations between the basic fractions. Thus the corresponding basic 
regularization parameters are linearly independent. 

We next specify the error of the above expansions. By not regularizing the bosonic 
potentials and fermionic wave functions, we clearly miss the 

higher orders in e/^macro • (5-9) 

Furthermore, in ()5.7p we must stay away from the origin, meaning that we neglect the 

higher orders in £/\(,\ ■ (5.10) 

The higher oder corrections in e/|^| depend on the fine structure of the regularization 
and thus seem unknown for principal reasons. Neglecting the terms in ()5.9p and (I5.10p 
also justifies the formal Taylor expansion in the degree. Finally, we disregard the 
higher order corrections in the parameter eshear in (13. 6p . 

The above symbolic computation rules give a convenient procedure to evaluate com- 
posite expressions in the fermionic projector, referred to as the analysis in the con- 
tinuum limit: After applying the contraction rules and expanding in the degree, we 
obtain equations involving a finite number of terms of the form ()5.6p . By applying 
the integration- by-parts rules, we can arrange that all simple fractions are basic frac- 
tions. We evaluate weakly on the light cone (j5.7p and collect the terms according 
to their scaling in ^. Taking for every given scaling in only the leading pole in e, 
we obtain equations which involve linear combinations of smooth functions and basic 
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regularization parameters. We consider the basic regularization parameters as empir- 
ical parameters describing the unknown microscopic structure of space-time. We thus 
end up with equations involving smooth functions and a finite number of free param- 
eters. We point out that these free parameters cannot be chosen arbitrarily because 
they might be constrained by inequalities (see the discussion after |14l Theorem E.l]). 
Also, the values of the basic regularization parameters should ultimately be justified 
by an analysis of vacuum minimizers of our variational principle (as discussed at the 
end of Chapter [3]) . 

In view of the later considerations in ^8.11 we point out that the above calculation 
rules are valid only modulo smooth contributions to the fermionic projector. This 
can be understood from the fact that these rules only deal with the terms of the 
series in (j4.17p . but they do not take into account the smooth non-causal high and low 
energy contributions. But the above calculation rules affect these smooth contributions 
as well. To give a simple example, we consider the distribution T^'^\ which according 
to (liJ3]l - ()il3D is given by 



Multiplying hy z = in the distributional sense gives a constant 



=r<«) = -^. (5.H) 



On the other hand, the contraction rule (15. 5j) yields 



The last relation gives much finer information than the distributional equation (jS.lip . 
which is essential when we want to evaluate composite expressions weakly on the light 
cone ()5.7p . However, the constant term in (j5.1ip does not appear in (|5.12p . The way 
to think about this shortcoming is that this constant term is smooth and can thus be 
taken into account by modifying the corresponding low energy contribution P'^(x,y) 
in (|4.17p . Indeed, this situation is not as complicated as it might seem at first sight. 
Namely, the smooth contributions to the fermionic projector need special attention 
anyway and must be computed using the resummation technique explained in Appen- 
dix |Dj When performing this resummation, we can in one step also compute all the 
smooth contributions which were not taken into account by the formalism of the con- 
tinuum limit. Thus altogether we have a convenient method where we first concentrate 
on the singularities on the light cone, whereas the neglected smooth contributions will 
be supplemented later when performing the resummation. 

We note that the above procedure needs to be modified for the description of gravity, 
because in this case we need relations between terms involving different powers of the 
Planck length e. We postpone this generalization to a future paper. 

5.2. The Euler-Lagrange Equations in the Continuum Limit. We now return 
to the action principle of Chapter [2j Our goal is to bring the conditions for a mini- 
mizer (12.10p and (12. lip into a form suitable for the analysis in the continuum limit. 
We begin by considering a smooth family P{t) of fermionic projectors and compute 
the corresponding first variation of the action. We differentiate (|2.1ip with respect 
to r, treating the constraint (j2.10p with a Lagrange multiplier. For convenience, we 



(V/)^ := _ „; - „: , (5.15) 
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introduce the functional 

S^iP]"'' [[ C^[A,y]d^xd^y with C^[A] = \A^\ - fi\A\\ (5.13) 
J Jmxm 

Choosing ^jl = \ gives precisely the action (j2.8p . whereas by ahowing a general /x G M 
we take into account the Lagrange multiplier. We thus obtain the condition 

= 5cS^[P]= // ReTr|v£/,[.4^y]5P(x,2/)|d^x A, (5.14) 

where 6P := -P'(O). Here we consider P(y,x) via 

P(y,:E) = P(x,y)* =7°P(x,y)V 

as a function of P{x,y), and V denotes the gradient where the real and imaginary 
parts of P{x, y) are treated as independent variables, i.e. 

Of . di_ 

dReP{x,y)a dlm.P{x,y)a 
and a, /3 = 1, . . . , 4 are spinor indices. Introducing the integral operator R with kernel 

R{y,x) ■.= VC^[A,y], (5.16) 

we can write (j5.14p as a trace of an operator product, 

6S^[P] =Retr(i? 6P) . 

In order to get rid of the real part, it is convenient to replace R by its symmetric part. 
More precisely, introducing the symmetric operator Q with kernel 

Q{x,y) = ^{R{x,y) + R{y,xr), (5.17) 

we can write the variation as 

SS^[P] = 2 tr {Q 5P) . (5.18) 

As explained before Definition 12. H we want to vary the fermionic projector by uni- 
tary transformations in a compact region. Thus the family of fermionic projectors P{t) 
should be of the form 

P{t) = U'\t)PU{t) (5.19) 
with a smooth family U{t) of unitary transformations in a fixed compact region K 
(see Definition 12. ip with U{0) = 1. Then the operator B = —iU'{0) has the integral 
representation 

{B^){x) = [ B{x,y)^{y)d^y 

with a smooth compactly supported integral kernel B G C^{K x ET, C*^^). Differen- 
tiating (j5.19p yields that 5P = i[P,B], and substituting this identity into (|5.18p and 
cyclically commuting the operators inside the trace, we can rewrite the condition (|5.14p 
as 

= tr([P,Q]P). 

Since B is arbitrary, we obtain the Euler-Lagrange (EL) equations 



[P,Q]=0, (5.20) 



stating that two operators in space-time should commute. For more details on the 
derivation of the EL equations we refer to |14i §3.5]. 
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When analyzing the commutator (j5.20p in the continuum hmit, the kernel Q{x,y) 
can be evaluated weakly using the formula (j5.7p . The subtle point is that, according 
to (j5.10p . this weak evaluation formula only applies if x and y stay apart. But writing 
the commutator in (j5.20p with integral kernels, 

[P,Q]{x,y)= [ (p{x,z)Q{z,y)-Q{x,z)P{z,y))d^z, (5.21) 

we also integrate over the regions z y and z ~ x where the kernels Q(z, y) and Q(x, z) 
are ill-defined. There are several methods to resolve this difficulty, which all give the 
same end result. The cleanest method is the method of testing on null lines. We now 
explain the ideas and results of this last method, referring for the rigorous derivation 
to Appendix |A] (for other methods of testing see [Ml Appendix F]). The idea is to take 
the expectation value of the commutator in ()5.2ip for two wave functions ^'i and ^^2) 
one being in the kernel and one in the image of the operator P. Thus 

P^-i = and ^2 = -P^ (5.22) 

for a suitable wave function Then, using the symmetry of P with respect to the 
indefinite inner product (12. 2p . we find 

<^'i I [P, Q] = <P^i I Q<^> - <^'i I QP^> =- <^i\ Q^2> ■ (5.23) 

Now the commutator has disappeared, and the EL equations (|5.20p give rise to the 
condition 

= <^'i|Q^'2>= // Q{x,y)^iix)^2iy)d^xdS. (5.24) 
J Jmxm 

The hope is that by choosing suitable wave functions ^'i and 4*2 of the form (I5.22p 
having disjoint supports, we can evaluate the expectation value (|5.24p weakly on the 
light cone (15. 7p . thus making sense of the EL equations in the continuum limit. 

The basic question is to what extent the constraints (j5.22p restrict the freedom in 
choosing the wave functions ^'i and ^'2. For clarity, we here explain the situation in 
the simplified situation where P is composed of one free Dirac sea of mass m, 

P{x, y)= [ (^ + m) (5(A:2 _ ^2) Q{-k^) g-^M^-s/) (5.25) 

J (2vr)'* 

(but Q can be a general operator for which the methods of Chapter [5] apply) . The 
generalization to several generations and a P with general interaction is worked out in 
Appendix El In order to extract information from (I5.24P and ()5.7p . it is desirable that 
the wave functions ^'i and 5*2 are as much as possible localized in space-time. For the 
wave function ^'i, this requirement is easy to fulfill by removing a strip of width Ao; 
around the lower mass shell in momentum space. For example, we can construct a 
wave function supported near the origin by choosing for a given parameter 5 > a 
smooth function r] supported in the ball of radius 5 in Euclidean and setting 




(5.26) 



where 77 is the Fourier transform of r/, and xi is the characteristic function defined 
by x/(a^) = 1 if 2; £ and xii^) = otherwise. In the limit Alo \ 0, the characteristic 
function in ()5.26p becomes the identity, so that ^'i goes over to r]. Moreover, for 
any Au > 0, the function ^'i is indeed in the kernel of the operator P, because it 
vanishes on the lower mass shell. Thus by choosing Aw sufficiently small, we can 
arrange that ^1 is arbitrarily close to rj and satisfies the condition in (j5.22p (indeed. 



AN ACTION PRINCIPLE FOR A FERMION SYSTEM AND ITS CONTINUUM LIMIT 



25 




Figure 1. Intersection of the null line £ with the singular set of Q{x,y). 



in finite space-time volume one cannot choose Ao; arbitrarily small, leading to small 
corrections which will be specified in Appendix [Aj see Remark IA.4p . 

The construction of ^'2 is a bit more difficult because ^'2 must lie in the image 
of P, and thus it must be a negative-energy solution of the Dirac equation (i^ — 
m)^2 = 0. Due to current conservation, it is obviously not possible to choose ^'2 to 
be localized in space-time; the best we can do is to localize in space by considering 
a wave packet. According to the Heisenberg Uncertainty Principle, localization in a 
small spatial region requires large momenta, and thus we are led to considering an 
ultrarelativistic wave packet of negative energy moving along a null line £, which does 
not intersect the ball -65(0) C where ^'i is localized. By a suitable rotation and/or 
a Lorentz boost of our reference frame (t,x), we can arrange that 

£ = {[r,-T + e, 0, 0) with r G M} 

with ^ > 0. For ^'2 we take the ansatz 

= {i^ + m) ( 6"*^^*+"^') (j){t + x- £, y, z)) + (small corrections) , (5.27) 



where the smooth function cj) is supported in 3^(0) C M^, and the frequency < as 
well as the length scales 6 and £ are chosen in the range 

e < ir^r^ < (5 < £,£macro,m"^ . (5.28) 

The small corrections in (j5.27p are due to the non-zero rest mass, the dispersion and 
the condition that ^2 must have no contribution of positive-energy (for details see 
Appendix . 

Except for the small corrections to be specified in Appendix [Aj the support of the 
wave function in (j5.26p lies in Bs{0), and thus it is disjoint from the support Bs{!^) 
of the wave function ^'2 in (I5.27P . Hence the integrals in (I5.24p only involve the region 
X ^ y where Q{x,y) is well-defined in the continuum limit. Furthermore, the null 
line 2, intersects the null cone around x in precisely one point y for which 1^*^! = |^| ~ ^ 
(see Figure [1]). Since this intersection is transverse, we can evaluate the expectation 
value (j5.24p with the help of (j5.7p . In view of the freedom in choosing the parameter £ 
and the direction of £, we conclude that (j5.7p itself must vanish. 



<5(x, y) = if evaluated weakly on the light cone . (5.29) 
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The above consideration is made rigorous in Appendix |Al More precisely, in Propo- 
sition IA.21 the above arguments are extended to the setting involving several gener- 
ations and a general interaction, and the scaling of the correction terms in (j5.27p is 
specified to every order in perturbation theory. This proposition applies to our action 
principle (|2.9p and all interactions to be considered here, thus justifying (j5.29p in all 
cases of interest in this paper. Moreover, in Remark IA.4I we consider the corrections 
to ()5.26p which arise if the lifetime of the universe is finite. Using that this lifetime 
can be estimated by the time from the big bang as known from experiments, we show 
that the correction to (j5.26p can indeed be neglected for our universe. 

To summarize, we saw that within the formalism of the continuum limit, the com- 
mutator in (j5.20p vanishes only if Q{x, y) itself is zero. This result is the strongest 
condition we could hope for, because in view of (j5.18p it implies that arbitrary first 
variations of the action vanish, even if we disregard the constraint that P must be a 
projector. We refer to (j5.29p as the Euler- Lagrange equations in the continuum limit. 

We finally remark that by replacing the null lines by null geodesies, the above 
method could immediately be generalized to situations involving a gravitational field. 
However, the estimates of Appendix Rl would become more demanding. 

6. The Euler-Lagrange Equations to Degree Five 

We proceed with the analysis of the EL equations in the continuum limit ()5.29p 
using the methods outlined in Chapters H] and \E\ For clarity, we begin in the vacuum 
and then introduce more and more interaction terms. Furthermore, we consider the 
contributions to the EL equations to decreasing degree on the light cone. In this 
chapter, we consider the most singular contributions of degree five. The contributions 
of degree four will be analyzed in Chapter [71 whereas the contributions to even lower 
degree are discussed in Chapter [H 

We point out that many results of this chapter are already contained in [lAi Chap- 
ters 5 and 6], albeit in a more general setting, while mainly restricting attention to 
one generation. In order to lay consistent foundations for the new calculations of 
Chapters EHHl we here present all calculations in a self-contained way. 

6.1. The Vacuum. To perform the light-cone expansion of the fermionic projector 
of the vacuum, we first pull the Dirac matrices out of the Fourier integral (14. ip and 
use (|4.12p to obtain 




After removing the logarithmic mass terms by the replacement Ta — s- T^*^^, the light- 
cone expansion reduces to a Taylor expansion in the mass parameter a. Restricting 
attention to the leading degree on the light cone, it suffices to consider the first term 
of this expansion. Using ()4.16p and taking the partial trace ()4.4p . we obtain for the 
regularized fermionic projector (for the factors of g see Footnote [2] on page [T9|) 

P(x,y) = |^r(y^) + (deg<2), (6.2) 

where for notational convenience we omitted the indices of the factor ^, and where 
the bracket (deg < 2) stands for terms of degree at most one. 
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Using this formula for the fermionic projector, the closed chain (|2.5p becomes 

-.2 



A^y = ^ i^T^-'^m^-'^) + ^(deg < 3) + (deg < 3) 



(6.3) 



where ^ := ^j^-'. Its trace can be computed with the help of the contraction rules (j5.4p . 



Tr{A,y) = (CjC^) T^y^^ T^y^^ = ^{z + z) Tl^'' T^^'' + (deg < 3) 
Next we compute the square of the trace-free part of the closed chain, 



(z + z) U + liz + zf 



^(-i)^(-i) 



16 

~ 64 ^ ' V M M 

Combining these formulas, we see that to leading degree, the closed chain is a solution 
of the polynomial equation 

(6.4) 



A 



9 



xy 



{z + z)T-'T, 



z-z) T^-i^T^-i) 



We point out that the calculations so far are only formal, but they have a well-defined 
meaning in the formalism of the continuum, because to our end formulas we will be 
able to apply the weak evaluation formula (15. 7p . Having this in mind, we can interpret 
the roots of the polynomial in ()6.4p 



\ _ 9 , (-1). (-1) 



and A_ = ^r(-^)(.r(-^)) 



as the eigenvalues of the closed chain. Using the contraction rule (jS.Sp . these eigenval- 
ues simplify to 



A, 



52r(°)r(7^) + (deg<3) 



A, 



52r(7^)r(';) + (deg<3). 



(6.5) 



The corresponding spectral projectors, denoted by F±, are given by 

Axy - A- 

^+ - a7^ ' 

a short calculation yields 



Axy A^ 

A_ - A+ 



1 



^ 2 V z-z 



3 +^(deg<0) + (deg<0). 



(6.6) 



Since in the formalism of the continuum limit, the factors z and 1 are treated as two 
different functions, we do not need to worry about the possibility that the denominator 
in (j6.6p might vanish. Similarly, we can treat ^ and ^ simply as two different vectors. 
Then the matrices -F+ and F- have rank two, so that the eigenvalues A+ and A_ are 
both two-fold degenerate. A straightforward calculation yields 



A^y = A+F+ + A_F_ + ^(deg < 3) + (deg < 3) , 



(6.7) 
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showing that our spectral decomposition is indeed complete. An important general 
conclusion from (j6.5p and (j6.6p is that in the vacuum, the eigenvalues of the closed 
chain form a complex conjugate pair, and are both two-fold degenerate. 

We now give the corresponding operator Q which appears in the EL equations of 
the continuum limit (j5.29p . 

Proposition 6.1. For the fermionic projector of the vacuum (|6.2p . the operator Q as 
defined by (|5.17p and (j5.16p takes the form 

Q{x,y) = iig' (1 - 4/^) T^JJ^T^y^) (deg < 5) . (6.8) 

In order not to distract from the main points, we first discuss the consequences of 
this result and derive it afterwards. According to the EL equations in the continuum 
limit ()5.29p . the expression ()6.8p must vanish. This determines the value of the La- 
grange multiplier ^ = \. Thus the action (I5.13j) reduces to the action in (|2.8I1 . and we 
conclude that 

P is a critical point of S , (6.9) 
disregarding the constraint T = const. This result can be understood immediately 
from the form of the Lagrangian ()2.7p and the fact that the eigenvalues of A^y form 
a complex conjugate pair. Namely, writing the spectral weights in ()2.7p via ()2.6p as 
sums over the eigenvalues (both of multiplicity two), we obtain 

£,.,[P] = (|A+|-|A_|)2 . 

The expression |A+| — |A_| clearly vanishes for a complex conjugate pair, and the fact 
that it appears quadratically is the reason why even first variations of Cxy[P] vanish, 
explaining ()6.9p . 

The last argument applies whenever the eigenvalues of A^y form a complex conjugate 
pair, making it possible to show that Q vanishes in a more general sense. First, a 
straightforward calculation shows that the eigenvalues of the closed chain Axy form a 
complex conjugate pair to every degree on the light cone (for details see [HI §5.3]), 
and thus Q vanishes identically in the formalism of the continuum limit. Moreover, 
going beyond the formalism of the continuum limit, in [18] it is shown that there are 
regularizations of the vacuum for which the operator Q vanishes up to contributions 
which stay finite in the limit e \ 0. Furthermore, in [18] it is shown that restricting 
attention to such regularizations does not give any constraints for the regularization 
parameters Creg in (15. 7p . Since we are here interested in the singularities of Q{x,y) in 
the limit e \ as described by the weak evaluation formula ()5.7p . we can in what 
follows assume that in the vacuum, the operator Q vanishes identically. 

The remainder of this section is devoted to deriving the result of Proposition l6.11 For 
the derivation it is preferable to bypass (I5.16P by determining Q directly from (I5.18p . 
For later use, we assume a more general spectral decomposition of Axy with eigenvec- 
tors A'^^, . . . , A^^ and corresponding one-dimensional spectral projectors 
This setting can be obtained from (j6.7p by choosing pseudo-orthonormal bases in the 
degenerate eigenspaces and letting F^^ be the projectors onto the span of these ba- 
sis vectors. It is convenient to choose these bases according to the following general 
lemma. 

Lemma 6.2. Assume that for a one-parameter family of fermionic projectors P{t) and 
fixed x,y G M , the matrices Axy and Ayx are diagonalizable for all t in a neighborhood 
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of T = 0, and that the eigenvalues of the matrix Axy\T=o are all non-real. Then the 
unperturbed closed chain A^y has a spectral representation 

^^.l.=o = EW^ (6-10) 

k=l 

with the following properties. The last two eigenvalues and spectral projectors are 
related to the first two by 

Af = Af, F3"^ = (Ff^)* and Af = Af , F^^ = (F^y)* . (6.11) 

The first order perturbation 5Axy = drAxy\T=o of the closed chain is diagonal in the 
bases of the non-trivial degenerate subspaces, i.e. 

F]^y{5A,y)F^y = ifk^l and Af = Af^ (6.12) 

The closed chain Ayx has a corresponding spectral representation satisfying (|6.10p - 
()6.12p with all indices 'xy ' are replaced by 'yx '. The spectral representations of A^y 
and Ayx are related to each other by 

Xiv = Af and F^^ P{x, y) = P{x, y) Ff . (6.13) 

Proof. By continuity, the eigenvalues of the matrix Axy are non-real in a neighborhood 
of r = 0. Using (|2.4p . one sees that the matrix Axy is symmetric in the sense that A^y = 
^xy = I^Aty^^. Hence, using the idempotence of the matrix 7*^ together with the 
multiplicity of the determinant, we find that 

det{Axy - A) = det{j^{Aly - A)7°) = det{Aly - A) = det(^^j^ - A) . (6.14) 

Hence if A is an eigenvalue of the matrix A^y, so is A. Thus the eigenvalues must form 
complex conjugate pairs. 

We first complete the proof in the case that there are no degeneracies. For any 
eigenvalue A of Axy we choose a polynomial px{z) with px{X) = 1 and p\{^j) = for all 
other spectral points /i. Then the spectral projector on the eigenspace corresponding 
to A, denoted by F^^ , is given by 

Fiy=px{Axy). (6.15) 

Taking the adjoint and possibly after reordering the indices /c, we obtain the rela- 
tions ()6.10p and ()6.1ip . The general matrix relation det(i?C — A) = det(Ci? — A) (see 
for example \17\ Section 3]) shows that the closed chains Axy and Ayx have the same 
spectrum. Multiplying ()6.15p by P{x, y) and iteratively applying the relation 

Axy P{x, y) = P{x, y) P{y, x) P{x, y) = P{x, y) Ayx , 

we find that F^^ P{x,y) = P{x,y) F^^ . Thus we can label the eigenvalues of the 
matrix Ayx such that (16.13P holds. 

In the case with degeneracies, the assumption that Axy is diagonalizable in a neigh- 
borhood of r = allows us to diagonalize 6Axy on the degenerate subspaces. This 
yields (j6.12p . whereas ()6.1ip can be arranged by a suitable ordering of the spectral 
projectors F^^ . In the degenerate subspaces of Ayx we can choose the bases such 
that ()6.10p and (|6.1ip hold (with 'xy' replaced by 'yx') and that ()6.13p is satisfied. 
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It remains to prove that ()6.12p also holds for Ay^: From (j6.12p we know that for any 
pair /, k with A^^ = A^^, 

= F^^6A^y)F^y = {5P{x, y) P{y, x) + P{x, y) 6P{y, x)) F^^ 

= F^y{5P{x,y))Fr P{y,x) + P(x, (5P(y, , 

where in the last line we applied the second equation in (j6.13p . Multiplying by P{y, x) 
from the left and by P{x,y) on the right, we find 

= P{y,x)F'^ymx,y))Fr Ar + XTFy'{5P{y,x))FrP{x,y) . 

After dividing by Xf^ = A^^ (note that the eigenvalues are non-zero because they 
are assumed to form complex conjugate pairs), we can again use the second equation 
in ()6.13p to obtain 

= P{y, x)F^y{6P{x, y))Fr + Fji'^idPiy, y) 

= Fl^{p{y,x) 6P{x,y)+6P{y,x) P{x,y))Fr = Fl^{5Ay,)Fr , 

concluding the proof. □ 

For later use, we next compute the operator Q in the general setting of the previous 
lemma. Noting that the function in (j5.13p depends only on the absolute values of 
the eigenvalues, we can write 

£,.[^.,] = £^(|A^^|i,...,|A^'^^|4). 
The partial derivatives of the function £^j(|A^^|i, . . . , |A^^|4) will be denoted by Dk- 

Lemma 6.3. Under the assumptions of Lemma 1 6. SX the operator Q in ()5.18p is given 
by 

A 

y) = Efc=i DkC,{\X^% . . . , I Af I) ^ Fi^y P{x, y) . (6.16) 

Proof. The relation (j6.12p allows us to compute the variation of the eigenvalues by a 
standard first order perturbation calculation without degeneracies, 

5\iy = i:T{F^^y5A,y). (6.17) 

Using that that 5\\\ = Re(A 5A/| A|), we can compute the first variation of this function 
with the help of ([6Trl) . 

4 

dC^^y] = Re^^_^ DkC,{\\iy\,. . . , |Af I) -1^ Tr(F,^'^ 5A^y) . (6.18) 

In the last trace we substitute the identity 

5A^y = 6P{x, y) P{y, x) + P{x, y) 6P{y, x) 
and cyclically commute the arguments to obtain 

Tr{F^y SA^y) = Tr {F^^ P{x, y) 5P{y, x) + P{y, x)F^y 6P{x, y)) 
= Tr {F^y P{x, y) 5P{y, x) + Ff P(y, x) SP{x, y)) , 
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where in the last step we applied (j6.13p . Substituting this formula into (j6.18p and 
integrating over x and y, we can exchange the names of x and y such that only 6P{y, x) 
appears. We thus obtain 

5S^[P] = 2Re jj^d^xdSTT{Q{x,y)5P{y,x)) (6.19) 

with the integral kernel Q{x,y) given by (j6.16p . Using Lemma 16.21 one sees that the 
operator corresponding to this integral kernel is symmetric (i.e. Q{x,y)* = Q{y,x)). 
As a consequence, the integral in (|6.19p is real, so that it is unnecessary to take the 
real part. Comparing with (j5.18p . we conclude that the operator with kernel (|6.16p 
indeed coincides with the operator Q in (|5.18p . We note that due to the sum in (j6.16p . 
it is irrelevant how the bases were chosen on the degenerate subspaces of Axy □ 

Proof of Proposition \6.1[ Let us specialize the general formula (j6.16p to our spectral 
representation with eigenvalues (|6.5p and spectral projectors (j6.6p . First, from ()5.13p 
we readily obtain that 

4 

I)fc£^(|Af I, . . . , |Af I) = 2|Afc| - 2/i^ |A,| = 2(1 - 4/x) |A_| . 

1=1 

The product F^^ P{x, y) can be computed with the help of (|6.2p and ()6.6p as well as 
the relations 

[tm = '^{LOi--^ei = -{z-z)i, 

where in the last step we treated the factors ^ and ^ as outer factors and applied the 
contraction rules ()5.3p and ()5.4p . We thus obtain 

F7P(x, y) = (deg < 2) , F!^P(x, y) = | ^ T^^y + (deg < 2) . (6.20) 
Substituting these formulas into (j6.16p and using (|6.5p . the result follows. □ 

6.2. Chiral Gauge Potentials. We now begin the study of interacting systems by 
introducing chiral potentials. Thus we choose the operator S in the auxiliary Dirac 
equation with interaction ()4.5p according to ()4.19p with two real vector fields Ai 
and A^. Sometimes it is convenient to write 23 in the form 

a3 = 4, + 7^4a (6.21) 
with a vector potential A^ and an axial potential A^ defined by 

A, = {Al + Ar)/2 and A^ = {Al - Ar) /2 . (6.22) 

To the considered highest degree on the light cone, the chiral gauge potentials merely 
describe phase transformations of the left- and right-handed components of the fermio- 
nic projector (for details see |13j or |14l §2.5]). More precisely, the fermionic projector 
is obtained from (16. 2p by inserting the phase factors 

P{x, y) = j {xL e-*^" + XR e-*^« ) ^ tI-'^ + (deg < 2) , (6.23) 
where the functions are integrals of the chiral potentials along the line segment xy, 

^T/R = r ^L/R ■■= f A^L/R\ryHl-r)^ " (6-24) 

J X Jo 
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Consequently, the closed chain is obtained from ()6.3p by inserting phase factors 
1 



Axy = ^ {XL VL + XR yn) ^\iT[^^ + ^(deg < 3) + (deg < 3) , (6.25) 
where 



VL = yR = e-^^^^'-'^^') = exp ( - 2i 4 . (6.26) 

Prom (j6.25p one sees that the matrix A^y is invariant on the left- and right-handed 
subspaces (i.e. on the image of the operators xl and XR)- On each of these invari- 
ant subspaces, it coincides up to a phase with the closed chain of the vacuum ()6.3p . 
Using these facts, the eigenvalues (Ag)cg{L,_R,},sG{+,-} and corresponding spectral pro- 
jectors Fl are immediately computed by 



)^J^ = vl/r Xs and F^^^ = xl/r F± 



(6.27) 

with and Fg as in (j6.5p and (j6.6p . We conclude that the eigenvalues of the closed 
chain are again complex, but in general they now form two complex conjugate pairs. 
Since the eigenvalues A^ and A^ differ only by a phase, we see that all eigenvalues 
have the same absolute value, 

|A^| = |A^| = |A^| = |A^| . (6.28) 



Writing the Lagrangian (12. 7p as 

C^y[P] = E\>^T-\{E\>^^\y = l E E (|A^|-|A^;|)' (6.29) 

c,s c,s c,c'£{L,R} s,s'e{±} 

(where we sum over c G {L, R} and s £ {i}), we find that C vanishes identically. Since 
the Lagrangian is quadratic in |A^| — |A^, |, also first variations of C vanish, suggesting 
that the operator Q{x,y) should again vanish identically. This is indeed the case, as 
is verified immediately by applying Lemmas 16.21 and 16.31 We conclude that for chiral 
potentials, the EL equations in the continuum limit (I5.29P are satisfied to degree five 
on the light cone. 

We end this section by explaining how the line integrals in ()6.24p and the phase 
factors in (j6.23p and (I6.25P can be understood from an underlying local gauge symme- 
try (for more details in the general context of non-abelian gauge fields see 114^ §6.1]). 
The local phase transformation ^(x) — )• e*^^^^^(x) with a real function A describes a 
unitary transformation of the wave functions (with respect to the inner product (]2.2p ). 
Transforming all objects unitarily, we obtain the transformation laws 

i^ + 'B-mY ^ e*^(^) (i^ + S - mY) e"*^^^) = + 'B - mY + (^A) (6.30) 

P{x, y) ^ e*^(^) P(x, y)e"*^(^) (6.31) 

A^y ^ e*^(") P(x, y)e-*^(j') e^^(^) P{y, x)e-^^(^) = A^y . (6.32) 

The transformation of the Dirac operator corresponds to a transformation of the vector 
and axial potentials by 

A^^A^ + dh and A^^A^. (6.33) 
These are the familiar gauge transformations of electrodynamics. Using the formula 

A(y) - A(x) = ^ ^A\ry+(^i_r)x dT = J\djA) ^ dr , 
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the phases in (j6.3ip can be described similar to (j6.24p in terms of line integrals. This 
explains why the phase factors in (j6.23p describe the correct behavior under gauge 
transformations. According to (j6.32p . the closed chain Axy is gauge invariant. This is 
consistent with the fact that in (j6.25p and (j6.26p only the axial potential enters, which 
according to (j6.33p is also gauge invariant. 

In order to transform the axial potential, one can consider the local transforma- 
tion ^'(x) —7- e'^^^ '^(^^^'(x). In contrast to the above gauge transformation, this trans- 
formation is not unitary (with respect to the inner product (12. 2p ). and the requirement 
that the Dirac operator and the fermionic projector must be symmetric operators leads 
us to the transformations 

i^ + 'B-mY ^ e^^AW (i^ + B - mY) e^^'^^^) (6.34) 

= i^ + e^^'^(^)(B - my)e^^'^(^) + 7^(^A) 

P{x,y) ^ e-^^'^(^) P(x,2/)e-^^'^(^) . 

Thus the vector and axial potentials transform as desired by 

and -> + OA 

(and also the term mY is modified, but this is of no relevance for the argument here). 
The point is that when we now consider the transformation of the closed chain, 

A^y e^^'^(^)p(x, y) e'^'^^y^ e'^'^^^'^ P{y, x) e^^'^(^) , (6.35) 

the local transformations do not drop out. This explains why in (j6.25p phases involving 
the axial potentials do appear. 

For clarity, we point out that the field tensors and the currents of the chiral gauge 
potentials also affect the fermionic projector, in a way which cannot be understood 
from the simple gauge transformation laws considered above. The corresponding con- 
tributions to the operator Q will be of degree four, and we shall consider them in the 
next chapter. 

7. The Euler-Lagrange Equations to Degree Four 

We come to the analysis of the EL equations to the next lower degree four on the 
light cone. In preparation, we bring the EL equations into a convenient form. 

Lemma 7.1. To degree four, the EL equations in the continuum limit (j5.29p are 
equivalent to the equation 

n := / t(°)t(-) ^ = + (deg < 4) , (7.1) 

where A denotes the perturbation of the eigenvalues (|6.27p to degree two. 

Proof. According to (j6.27p . the eigenvalues to degree three are all non-real. Since this 
property is stable under perturbations of lower degree, we can again apply Lemmas l6.2l 
and 16.31 Noting that before (|6.9p . we fixed the Lagrange multiplier to = we 
consider the Lagrangian (|2.7p . which we now write in analogy to (|6.29p as 

C.y[P] = \j2{\K'\-\>H'\f- 
k,l=l 
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Then the relation (j6.16p can be written as 

According to (16.281) . the curly brackets vanish for the unperturbed eigenvalues. This 
has the convenient consequence that to degree four, it suffices to take into account 
the perturbation of the curly brackets, whereas everywhere else we may work with the 
unperturbed spectral decomposition (|6.27p . 

Q{x,y) = i ^ A(|Af I - lAf^l) ^ F^^ P{x,y) + (deg < 4) . 

Using (j6.20p . we see that we only get a contribution if Afc equals or X^. Furthermore, 
we can apply (16. lip , numbering the eigenvalues such that A^ — We thus obtain 

Q{x,y)= A(|Al|-|A;|)^Xc|5r(y'^ + (deg<4). (7.2) 

ce{L,R} ' 

The EL equations (|5.29p imply that the left- and right-handed components of this 
expression must vanish separately. Thus, again applying (j6.1ip . we obtain the sufficient 
and necessary condition 

a(|A^| - IX^l^^^gTl^-'^ + (deg < 4) = . 

The explicit formulas (j6.27p and (j6.5p yield the result. □ 

It is important to observe that the EL equations only involve the difference of the 
absolute values of the left- and right-handed eigenvalues. This can immediately be 
understood as follows. To the leading degree three, the eigenvalues of A^y form two 
complex conjugate pairs (see ()6.27p ). Since this property is preserved under pertur- 
bations, we can again write the Lagrangian in the form (|6.29p . Hence the Lagrangian 
vanishes identically unless the absolute values of the eigenvalues are different for the 
two pairs. This explains the term A(|A:?^| — \\^\) in ()7.ip . 

As explained on page [28l the expression A(|A:?^| — \X^\) vanishes in the vacuum. 
Furthermore, the phase factors in ()6.27p drop out of this expression. But new types of 
contributions to the interacting fermionic projector come into play, as we now explain. 

7.1. The Axial Current Terms and the Mass Terms. An interaction by chiral 
potentials (|4.19p as introduced in §6.21 affects the fermionic projector in a rather com- 
plicated way. For clarity, we treat the different terms in succession, beginning with 
the contributions near the origin ^ = (the contributions away from the origin will be 
considered in Chapter [TOl) . For the Taylor expansion around = we note that when 
evaluated weakly on the light cone (|5.7p . a simple fraction of degree L has a pole |'^|~'^. 
Thus it is reasonable to say that a term of the form (j5.6p is of the order k at the origin 
if the smooth function rj vanishes at the origin to the order k + L. 

Definition 7.2. An expression of the form (j5.6p is said to be o/ order o(|^|^) at the 
origin if the function rj is in the class o((|^^| + ICI)'^^'^)- 
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In the next lemma we specify the contributions to the EL equations to degree four on 
the hght cone, to leading order at the origin. 

Lemma 7.3. For an interaction described by vector and axial potentials (|4.19|) , the 
expression IZ as defined by (|7.1|) takes the form 

^ = -iik (ja iVi - m^Al N2) + (deg < 4) + 0^-^) , (7.3) 
where ja is the axial current 

jt = d)A^,-UAl, (7.4) 



and Ni , N2 are the simple fractions 
-,3 



No 



^[0] [0] -'[I] -'[I] ^[0] [2] -^[0] ^[01 ) 



(7.5) 
• (7.6) 



Here "c.c." denotes the complex conjugate of the preceding simple fraction; the accents 
were defined in (j5.2p . 



In order not to distract from the main ideas, we postpone the proof of this lemma to 
Appendix[B]and proceed right away to the physical discussion. Prom the mathematical 
point of view, the appearance of the axial current js. is not surprising, because the 
light-cone expansion of the fermionic projector involves derivatives of the potentials. 
In physical terms, this shows that the axial potential affects the fermionic projector not 
only via the phases in ()6.23p . but also via the Yang-Mills current. The term —i£,kja^i 
is referred to as the current term. The other term —i^k ni?A^N2 could not appear 
in ordinary Yang-Mills theories because it would not be gauge invariant. However, as 
pointed out after (j6.34p . the axial ?7(l)-transformations do not correspond to a local 
gauge symmetry, because they are not unitary. Instead, they describe relative phase 
transformations of the left- and right-handed components of the fermionic projector, 
thereby changing the physics of the system. Only the phase transformations (|6.32p 
correspond to a local gauge symmetry, and in view of (f03]) . the term -i^t m^yl^iVa 
is indeed consistent with this local U (l)-symmetry. 

Since the direction ^ can be chosen arbitrarily on the light cone, the condition ()7.ip 
implies that the bracket in (j7.3p must vanish, 

jtN^-m^AlN2 = Q. (7.7) 

If A^i and could be treated as constants, this equation would go over to a field 
equation for the axial potential with rest mass m?'N2/Ni. For this reason, we refer 
to the term —i^k m'^A^N2 in (17. 3p as the mass term. It is remarkable that in our 
framework, the bosonic mass term appears naturally, without the need for the Higgs 
mechanism of spontaneous symmetry breaking (for a detailed discussion of this point 
see §8.5p . We also point out that the simple fraction A^2 involves the mass matrix y, 
and thus the mass term in (|7.7p depends on the masses of the fermions of the system. 

In order to make the above argument precise, we need to analyze the simple frac- 
tions Ni and weakly on the light cone. Before this will be carried out in ^7.31 we 
specify how the Dirac current enters the EL equations. 
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7.2. The Dirac Current Terms. As explained in ^4.31 the particles and anti-particles 
of the system enter the auxiliary fermionic projector via (j4.10p . where we orthonor- 
malize the wave functions according to (|4.1ip . Introducing the left- and right-handed 
component of the Dirac current by 

"/ na 
k=l 1=1 

a decomposition similar to (|6.22|) leads us to define the axial Dirac current by 

"/ ria 

4 = 5^^7V^fc-^^7V^«- (7.8) 

k=l 1=1 

The next lemma gives the corresponding contribution to the EL equations, to leading 
order at the origin. 

Lemma 7.4. Introducing the axial Dirac current by the particle and anti-particle wave 
functions in (j4.10p leads to a contribution to TZ of the form 

n X iCfc 4 Ns + (deg < 4) + odel^^^) , 

where 

N, = C= \tL'^ WtL^ - C.C.1 . (7.9) 



Sir j,{0) 



tL''>T,^'^tL-'^-c.c. 



Here the symbol "x" means that we merely give the contribution to TZ by the Dirac 
current, but do not repeat the earlier contributions given in Lemma 17.31 The proof of 
this lemma is again postponed to Appendix [Bl 

7.3. The Logarithmic Poles on the Light Cone. Combining the results of Lem- 
mas I7.H 17.31 and 17.41 the Euler-Lagrange equations give rise to the equation 

Ck(j'.N,-m'AtN2-4Ns)=0, 



which involves the axial potential Ag, (see (|6.2ip ). the corresponding Yang-Mills cur- 
rent (j7.4p and the axial Dirac current (j7.8p . At first sight, this equation resembles a 
bosonic field equation, which describes the coupling of the Dirac spinors to the bosonic 
field and involves a bosonic mass term. However, the situation is not quite so simple, 
because the factors A^i, and N3 (see ()7.5p . ()7.6p and (I7.9p ) have a mathematical 
meaning only when evaluating weakly on the light cone ()5.7p . Let us analyze the weak 

evaluation in more detail. The simple fraction N3 is composed of the functions T'jop 

TjQj and their complex conjugates, which according to (|4.13p - (|4.16p all have poles 

of the order or In particular, no logarithmic poles appear, and thus we may 
apply (15. 7p with r = to obtain 

r\(\+e reg 

I dtri^k Ja N3 = -f^ v{x) 6 Jaix) + (deg < 4) + 0^-^) 

with a regularization parameter c^^^, where we omitted error terms of the form (|5.9p 
and (|5.10p . The simple fractions Ni and N2, on the other hand, involve in addition 

the functions tJ^"* and To^\ which according to (|4.13p - (|4.15p involve a factor log|^^| 
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and thus have a logarithmic pole on the light cone. As a consequence, in ()5.7p we also 
obtain contributions with r = 1, 



dtri^k 



(deg < 4) + o(|er 



+ 



d{'^ log(elCl) - m'A', + dl^^ log(elCI) 



involving four regularization parameters and d^^- Combining the above weak eval- 
uation formulas, the freedom in choosing the radius |^| and the spatial direction (,/\(,\ 
implies that that the logarithmic and non-logarithmic terms must vanish separately. 



k reg 



-3 



(7.11) 

where and are constants depending on the particular regularization. 

Unfortunately, the system of equations ()7.10p into ()7.1ip is overdetermined and is 
thus too restrictive for physical applications. Namely, at least for generic regulariza- 
tions, the constants c^^^, c^^ and d^^^ are non-zero. Thus solving ()7.10p for and 
substituting into ()7.1ip . one obtains an algebraic equation involving Ja and Ag,. This 
means that either Ja must vanish identically, or else the gauge potential Ag^ is fixed to 
a constant times Ja and thus cannot be dynamical. Both cases are not interesting from 
a physical point of view. The basic reason for this shortcoming is that the bosonic 
current and mass terms have logarithmic poles on the light cone, whereas the Dirac 
current terms involve no such logarithms. Our method for overcoming this problem 
is to insert additional potentials into the Dirac equation, with the aim of compen- 
sating the logarithmic poles of the bosonic current and mass terms. Before entering 
these constructions in ^7.41 we now briefly discuss alternative methods for treating the 
logarithmic poles. 

An obvious idea for reducing the system (I7.10p and (I7.1ip to a single equation is 
to restrict attention to special regularizations where the constants c^*^^ and/or d^-'^^ 
take special values. In particular, it seems tempting to demand that d^-^^^ = d^2^ = 0, 
so that ()7.10p is trivially satisfied, leaving us with the field equations ()7.1ip . This 
method does not work, as the following consideration shows. Differentiating ()4.14p 
and using (I4.13p . one sees that 

327r3 tP = log le^l + CO + inOiC^) e{f) . (7.12) 

Evaluating near the upper light cone w |.^|, we can apply the relation log|.^^| = 



log 



+ log 

.(1) 



to obtain 



327r^ rr = log |2^i + log e - lei + i^eie - +co + o{e - lei) • (7.13) 



When evaluating the corresponding simple fraction weakly ()5.7p . the first term in ()7.13p 
gives rise to the log |^|-dependence, the second term gives the log e-dependence, whereas 
all the other terms do not involve logarithms or are of higher order in e. Obviously, 

the same is true for the complex conjugate To^\ Since in (15. 7p the vector ^ is fixed, 
the vanishing of the log |^|-dependent contribution to the integral ()5.7p implies that 

the simple fraction still vanishes if the factors tJ^^ and tJ^'* are replaced by constants. 
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Inspecting the T^^^-dependence of (j7.5p and (j7.6p and comparing with (j7.9p . we find 
that 

d-g = ^ 4'*^ = and <s ^0 ^ cl'^ = 0. 

Thus if the constants (f^^ and (ig'^^ in (|7.1Up vanish, then (|7.10p becomes trivial as 
desired. But then the constant c^^ in (|7.1ip is also zero, so that the Dirac current 
drops out of the field equation. Again, we do not end up with physically reasonable 
equations. 

Sticking to the idea of considering regularizations where the regularization constants 
have special values, the remaining method is to assume that oZZ regularization constants 
in ()7.10p and (j7.1ip vanish. Then the EL equations would be trivially satisfied to degree 
four on the light cone, and one would have to proceed to the analysis to degree three on 
the light cone. This method does not seem to be promising for the following reasons. 
First, it is not clear whether there exist regularizations for which all the regularization 
constants in (j7.10p and (17. lip vanish. In any case, it seems difficult to satisfy all these 
conditions, and the resulting regularizations would have to be of a very special form. 
This would not be fully convincing, because one might prefer not to restrict the class 
of admissible regularizations at this point. Secondly, there is no reason to believe that 
the situation to degree three would be better, at least not without imposing additional 
relations between regularization constants, giving rise to even more constraints for the 
admissible regularizations. 

We conclude that assuming special values for the regularization constants in ()7.10p 
and ()7.1ip does not seem to be a promising strategy. Thus in what follows we shall not 
impose any constraints on the regularization constants, which also has the advantage 
that our constructions will apply to any regularization. Then the only possible strategy 
is to try to compensate the logarithmic poles by additional potentials in the Dirac 
operator. 



7.4. A Pseudoscalar Differential Potential. Our goal is to compensate the loga- 
rithmic poles of the bosonic current and mass terms by inserting additional potentials 
into the auxiliary Dirac equation (|4.5p . In order to get contributions of comparable 
structure, these potentials should involve a vector field v, which should be either equal 
to the axial potential Ag, or the corresponding axial current ja (see Lemma [7. 3p . Since 
contracting the vector index of v with the Dirac matrices would again give rise to chi- 
ral potentials (I4.19p . we now prefer to contract v with partial derivatives. Moreover, 
since we want to compensate contributions which are odd under parity transformations 
(i.e. which flip the left- and right-handed components), the resulting operator must 
involve the pseudoscalar matrix 7^. The requirement that the Dirac operator should 
be symmetric with respect to the inner product (12. 2p leads us to the ansatz involving 
an anti-commutator 

S = 75 I „i ^ Q. } = dj + 7^ {dj ) . (7.14) 

We refer to this ansatz as a pseudoscalar differential potential. Our ansatz seems 
unusual because such differential potentials do not occur in the standard model nor in 
general relativity. We postpone the physical discussion to the last paragraph of ^7.61 
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The corresponding leading contribution to the fermionic projector is of the form (for 
details see equation ()B.32p in Appendix [B|) 



P{x,y)^^j%{v\y) + v\x))T^-'^ 



+ [lim]T^-'^ + (deg<2). 



(7.15) 



This contribution has a pole of order on the light cone and is therefore much 
more singular than the desired logarithmic pole. A straightforward calculation shows 
that ()7.15p does contribute to the expression TZ in Lemma |7. 11 and thus we conclude 
that (|7.15p is not suitable for compensating the logarithmic pole. 

The key for making use of the pseudoscalar differential potential (j7.14p is to observe 
that the required logarithmic poles do appear to higher order in a mass expansion. 
More precisely, to leading order at the origin, the cubic contribution to the fermionic 
projector is 



3 



P{x,y) X — 7^ [vf\x) + 0{\e\ + |e1)J (^e^ - 27^ r(i)) + (deg < -1) , (7.16) 

where v^^^ is a Hermitian matrix composed of v and Y, 

t;(3) = i (vYYY - YvYY + YYvY - YYYv) (7.17) 

(for details see equation ()B.33p in Appendix [B|). Thus there is hope that the logarith- 
mic poles can be compensated, provided that we can arrange that the contributions 
by ()7.14p to TZ of order tvP, m and tti^ in a mass expansion vanish. The last re- 
quirement cannot be met if we consider one Dirac sea, because the term (|7.15p does 
contribute to TZ. But if we consider several Dirac seas, we have more freedom, as the 
pseudoscalar differential potential (j7.14p can be chosen differently for each Dirac sea. 
For example, we can multiply the potentials acting on the different Dirac seas by real 
constants Qa, 

{'Br^=Qa6a(3 7'{v\dj} with o, /3 = 1 , . . . , 5 • (7.18) 

Using this additional freedom, it is indeed possible to arrange that the contribu- 
tion (j7.15p drops out of TZ. This consideration explains why we must consider several 
generations of elementary particles. 

The critical reader might object that there might be other choices of the operator 23 
which could make it possible to compensate the logarithmic poles without the need for 
several generations. However, the following consideration shows that (j7.14p is indeed 
the only useful ansatz. First of all, since we want to compensate a contribution to the 
fermionic projector at the origin, it seems unavoidable to consider local operators (for 
nonlocal operators see Chapter llOp . The only zero order operator are the chiral poten- 
tials (j4.19p . which were already considered in ^6.2[ Apart from (j7.14p . the only first 
order differential operator involving the vector field v and the pseudoscalar matrix 7^ 
is the operator 

where a^^ = §[7-' , 7^^] are the bilinear covariants. This ansatz can be shown to be use- 
less, basically because the calculations in the continuum limit give rise to contractions 
with the vector ^, which vanish (see also §9.2p . Differential operators of higher order 
must involve the wave operator □, which applied to the Dirac wave functions gives 
rise to lower order operators. This shows that it is not useful to consider differential 
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operators of order higher than one. We conclude that (j7.14p and its generalizations to 
several generations (like ()7.18|) ) are indeed the only possible ansatze for compensating 
the logarithmic poles. 

We end the discussion by having a closer look at the matrix v^^\ (j7.17p . Note 
that the ansatz (|7.18|) is diagonal in the generation index and thus commutes with 
the mass matrix Y. As a consequence, the matrix v^^^ vanishes. This means that for 
compensating the logarithmic poles, the ansatz (j7.18p is not sufficient, but we must 
allow for non-zero off-diagonal elements in the generation index. Thus we replace the 
factors ba in (j7.18p by a Hermitian matrix g = (0^)a,/3=i,...,g, the so-called generation 
mixing matrix. Later on, the generation mixing matrix will depend on the space-time 
point X. This leads us to generalize (|7.19p by the ansatz 

i15)$ = j^Q'^^ix)v^ix),dj} , (7.19) 

thus allowing that the pseudoscalar differential potential mixes the generations. 

7.5. A Vector Differential Potential. Modifying the auxiliary Dirac equation (|4.5p 
by a first oder operator (j7.14p or (j7.19p changes the behavior of its solutions drastically. 
In particular, it is not clear whether the operator 23 can be treated perturbatively ()4.9p . 
In order to analyze and resolve this problem, we begin by discussing the case where 
the potential v in (|7.14p is a constant vector field, for simplicity for one Dirac sea 
of mass m. Then taking the Fourier transform, the Dirac equation reduces to the 
algebraic equation 

(^ - 2i-f^v^kj - m) ^(/c) = . (7.20) 

Multiplying from the left by the matrix {I^ — 2ij^v^kj +m), we find that the momentum 
of a plane-wave solution must satisfy the dispersion relation 

k^ - 4:{v^kj f - = . 

Rewriting this equation as 

g'^^kikj - = with g"-^ := rj"^ - Av'v^ , 

where rf^ = diag(l, —1, —1, —1) is the Minkowski metric, we see that the new disper- 
sion relation is the same as that for the Klein-Gordon equation in a space-time with 
Lorentzian metric g'^K In particular, the characteristics of the Dirac equation become 
the null directions of the metric g^K In other words, the light cone is "deformed" to 
that of the new metric g^^ . 

This deformation of the light cone leads to a serious problem when we want to 
compensate the logarithmic poles, as we now discuss. Suppose that we introduce a 
pseudoscalar differential potential which according to ()7.18p or (I7.19P depends on the 
generation index. In the case (j7.18p . the Dirac seas feel different dispersion relations. 
In particular, the singularities of the fermionic projector P{x, y) will no longer be 
supported on one light cone, but will be distributed on the union of the light cones 
corresponding to the Lorentzian metrics g]l = rf^ — 4:vl^Va. The ansatz (I7.19P leads to 
a similar effect of a "dissociation of the light cone." In the EL equations, this would 
lead to large additional contributions, which are highly singular on the light cone and 
can certainly not compensate the logarithmic poles. 

Our method for bypassing this problem is to introduce another differential potential 
which transforms the dispersion relation back to that of the Klein-Gordon equation in 
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Minkowski space. In the case of a constant vector field v and one generation, tliis can 
be achieved by choosing matrices which satisfy the anti-commutation relations 

[G\G^] = 2rf^ + Sv'v^ and {^^ , (x) } = , 

and by modifying (|7.20p to 

[G^kj - 2ir>v^kj -m)i!{k) = Q. 

This modification of the Dirac matrices can be interpreted as introducing a constant 
gravitational potential corresponding to the metric rf^ + 4vl^Va. This construction is 
extended to the general case (j7.19p as follows. We choose {Ag x 4(7)-matrices G^{x) 
which are symmetric with respect to the inner product on the Dirac spinors and 
satisfy the anti-commutation relations 

{G\x),G^{x)] =2rf^ +SQ{xfv\x)v\x) and {7^G^(x)}=0. (7.21) 

In the auxiliary Dirac equation (14. Sp we insert the additional operator 

'B = i{G^ {x) - -i^) dj + G^ (x) Ej (x) , (7.22) 

where the matrices Ej involve the spin connection coefficients and are not of importance 
here (for details see for example [m §1.5]). We refer to (|7.22p as a vector differential 
potential. In the case (|7.18p . this construction can be understood as introducing for 
each Dirac sea a gravitational potential corresponding to the metric r/*-^ -|- 40^t>^t>o, 
whereas in case (|7.19p . the interpretation is bit more complicated due to the off- 
diagonal terms. 

7.6. Recovering the Differential Potentials by a Local Axial Transformation. 

By introducing the differential potentials (|7.19p and (j7.22p with G^ according to (j7.2ip . 
we inserted differential operators into the auxiliary Dirac equation (j4.5p . We will now 
show that the effect of these operators on the solutions of the auxiliary Dirac equation 
can be described by a local transformation 

^aux(2:) ^ U{x) ^'aux(x) , (7.23) 

which is unitary with respect to the inner product (14. 6p . This simplification also allows 
us to compute the resulting logarithmic pole in detail. We conclude this chapter by a 
physical interpretation of the differential potentials and of the transformation ()7.23p . 

Recall that we introduced the vector differential potential ()7.22p with the goal of 
transforming the dispersion relation back to the form in the vacuum. Thus if u is a 
constant vector field, the combination (I7.19p -l- (l7.22p leaves the momenta of plane-wave 
solutions unchanged. This suggests that the sum (|7.19p -l- (l7.22p might merely describe 
a unitary transformation of the Dirac wave functions. Thus we hope that there might 
be a unitary matrix U{x) such that 

U{i^ - mY)U~^ = i^+ (TM + . 

Let us verify whether there really is such a unitary transformation. The natural ansatz 
for U is an exponential of an axial matrix involving the vector field v and the generation 
mixing matrix, 

U{x) = exp (-ig(x) ^^^^Vj{x)) . (7.24) 
Writing out the exponential series and using that (7^7-' fj)^ = —v^, we obtain 

U{x) = cos(g(^) 1 - i ^"^^^'^^ 7^^ , U{x)-^ = cos(g(^) 1 + i ^"^^^'^^ 7^^ , (7.25) 
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where the angle ip := \/ —v^ is real or imaginary (note that (17.250 is well-defined even 
in the limit (/? — t- 0). A short calculation yields 

and thus 

U{i^ - mY) U'^ = iU-f^U'^dj + Uj^ {idjU''^) - mUYU'^ 
= i^ + 2^^v^ ^M^al [/-la- + U-i\idjU-^) - mUYU-^ 

= ii^ + ^5 ^,Q^ + 2i Ei!!^ ^ yjg^ + U^^iidjU-') - mUYU-' . (7.26) 

In order to verify that the resulting Dirac operator allows us to recover both (I7.19P 
and (I7.22[) . we assume that v"^ is so small that sm{2Qip) w 2Qip and sin^(gv3) 
0^(/?^. Then the second summand in (|7.26p reduces precisely to the difi^erential op- 
erator in (j7.19p . The third summand in (I7.26P gives precisely the differential operator 
in (IL22]), noting that (1731]) has the solution = 7^' + 2q^^v^ + 0(7;^). Likewise, a 
direct calculation shows that the multiplication operators in (I7.19P and (I7.22P are con- 
tained in the fourth summand in (I7.26p . Writing out the fourth and fifth summands 
in (]7.26p . one finds a rather complicated combination of additional chiral, scalar, pseu- 
doscalar and even bilinear potentials. These additional potentials do not cause any 
problems; on the contrary, they guarantee that the total transformation of the Dirac 
wave functions simply is the local transformation (j7.23p . We conclude that with ()7.26p 
we have found a Dirac operator which includes the differential potentials in (j7.19p 
and ()7.22p . It has the nice property that it can easily be treated non-perturbatively by 
the simple local transformation ()7.23p . We refer to the transformation (j7.23p with U 
according to ()7.24p as the local axial transformation. 

To avoid confusion, we point out that we shall always perform the transforma- 
tion (j7.23p after all the other potentials have been inserted into the auxiliary Dirac 
equation. For example, to combine the chiral potentials ()4.19p with the differential 
potentials, we consider the Dirac operator [/(z^+Xi^ij + XR^L""*-^)^"^- The corre- 
sponding auxiliary fermionic projector has the form U{x)P{x, y)U~^{y), where P{x, y) 
is defined by the perturbation series (14. 9p with 23 composed only of the chiral potentials. 
It might seem puzzling that the transformation U also changes the zero order terms 
in the Dirac operator, which thus no longer seem to be under full control. However, 
the unitary transformation poses no restriction for the zero order terms in the final 
Dirac operator, because any undesirable multiplication operator B in the final Dirac 
operator can be compensated simply by inserting the operator —U~^BU into the Dirac 
operator before the unitary transformation. In view of the exhaustive discussion of 
different multiplication operators in ^6.21 and ^9.1f^9.31 the zero order terms in the 
Dirac equation can be changed arbitrarily, and the effect on the fermionic projector is 
well-understood. 

We next work out how the local axial transformation can be used to compensate the 
logarithmic poles. For simplicity, we only consider a perturbation expansion to first 
order in v and remark that the corresponding non-perturbative treatment is given in 
Appendix [Cj To first order in v, the transformation (|7.24p simplifies to 

U{x) = l-ig-/^^{x) + 0{v^) . (7.27) 
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Transforming the auxiliary fermionic projector (j4.ip by U and taking the partial 
trace (j4.4p . we obtain for the perturbation of the fermionic projector the expression 

P >i -i-f^^Qp + iPQ-f^'^ + (D{v'^) , (7.28) 

where we denoted the partial trace similar to (j5.2|) by accents. Here we always sum 
over one index of the generation mixing matrix. Thus it is convenient to introducing 
real functions Cq, and da by 

9 9 

= cp + idp and ^ = - ida , (7.29) 

a=l /3=1 

where the last equation is verified by taking the adjoint of the first and using that q 
is Hermitian. Combining these equations with the fact that the auxiliary fermionic 
projector of the vacuum is diagonal on the generations, we can write ()7.28p as 

9 

P>^Y{-ihl'i,Pp]+{dp7''^,Pp})+0{v''), (7.30) 

(3=1 

where the Pg stand for the direct summands in (|4.ip . The next lemma shows that the 
functions C/j drop out of the EL equations; the proof is again given in Appendix [Bl 

Lemma 7.5. The perturbation of the fermionic projector by the functions cp in ()7.30p 
drops out of the EL equations to degree five and four on the light cone. 

Noting that the diagonal elements of g only contribute to the functions cp (as is obvious 
from (j7.29p and the fact that the diagonal elements of g are real), this lemma again 
shows (in analogy to the consideration leading to the ansatz (|7.19p ) that the mixing 
of the generations is essential for our constructions to work. 

Using the result of the last lemma, it remains to analyze the effect of the functions d^. 
Computing the leading degrees on the light cone, we obtain 



P(x,y)x ^ 
/3=1 



|d^7'^, I + mp T(o)| + (deg < 1) 



\ 7V(x) i r(-i) ( ^ dp{x)) - \ i,{y) r(-i) ( dp{y)) (7.31) 



=1 

9 



+ (7'^(^) + ^'^(2/) Mv)) T^"^ + (deg < 1) . (7.32) 

13=1 (3=1 

Since q is Hermitian, summing the first equation in (I7.29P over /3 gives a real number. 
Hence 

9 

^d/3 = 0, (7.33) 

(3=1 

so that (j7.3ip vanishes. The term in (17.32P has a similar form as the contribution by 
the axial current to P{x,y) (for details see equation (IB.ISP in Appendix [B]), but on 
the light cone it has instead of a logarithmic pole a stronger singularity ~ For 
general regularizations, this term does contribute to the EL equations even to degree 
five on the light cone. This is made precise in the next lemma, which is again derived 
in Appendix [Bl 
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Lemma 7.6. Under the assumptions (j7.33p . the perturbation of the fermionic projector 
by the functions dp in (j7.30p leads to a contribution to TZ of the form 

n^^i,v\Y. dp) ^ r^y^)^;^ [ri^^^ - c.c.) + (deg < 5) + am-') . 

Since we do not want to put any conditions on the regularizations, the EL equations 
yield the constraint 

9 

Y,mpdp = {). (7.34) 

/3=1 

The remaining contributions to the fermionic projector are of degree at most zero on 
the hght cone. They affect the EL equations to degree four, as specified in the next 
lemma, which is again proved in Appendix [Bl 

Lemma 7.7. Under the assumptions (j7.33p and (|7.34p . the perturbation of the fermio- 
nic projector by the functions dp in (|7.30p leads to a contribution to TZ of the form 

TZ X iCk v'' iV4 + (deg < 4) + o(\^\-^) , 

where 

9 



2g' ' 



+ = E Ti.%-y-y^^ - C.C. . (7.36) 



We point out that the simple fraction A'4 involves factors TjjjJ^ and Tj^^ , which in view 
of (j7.12p have the desired logarithmic pole on the light cone. This indeed makes it 
possible to compensate the logarithmic poles of the current and mass terms, as will be 
explained in the next chapter. 

We conclude this chapter by a brief physical discussion. With the previous con- 
structions we recovered the differential potentials of ^7.41 and ^7.51 by a local axial 
transformation of the auxiliary Dirac operator and the auxiliary fermionic projector 
of the form 

{i^ + 'B- mY) U{i^ + B - mY)U-^ , P^'^^(x, y) U{x) P^'^^(x, y) U{y)'^ , 

where U{x) is unitary with respect to the inner product (|4.6p . Here the operator S 
is composed of all bosonic potentials and fields except for the differential potentials. 
It is important to observe that the local axial transformation is performed before the 
partial trace is taken (14. 4p . In particular, it cannot be described by a local unitary 
transformation of the fermionic projector 

P(x, y) U{x) P{x, y) U{yy^ with U G U(2, 2) . 

Indeed, the last transformation would simply lead to a unitary transformation of the 
closed chain A^y — )■ U {x)AxyU (x)~^, leaving the spectrum of the closed chain and thus 
also the Lagrangian unchanged. The local axial transformation, however, does affect 
the Lagrangian, because the spectrum of the closed chain is computed after the partial 
trace has been taken. In view of this fact, the local axial transformation cannot be 
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interpreted merely as a local gauge transformation, but it really changes the physical 
system. Nevertheless, due to the local nature of the transformation, the local axial 
transformation does not change the dynamics of the physical system and is thus not 
directly observable. The physical picture is that the local axial transformation slightly 
modifies the fermionic wave functions at every space-time point, thus affecting the 
local form of the fermionic projector. 

We finally remark that the local axial transformation also influences the inner prod- 
uct (j2.2p between fermionic states, an effect which could in principle be physically 
observable. However, since the integration in ()2.2p extends over all space-time, it is 
not clear how the inner product (j2.2p could be determined in experiments. It thus 
seems that the local axial transformation does not lead to any conceivable measurable 
effects. 

8. The Field Equations 

Having developed a method for compensating the logarithmic poles on the light 
cone, we are now in the position to derive and analyze the field equations. 

8.1. The Smooth Contributions to the Fermionic Projector at the Origin. 

We add the contributions from Lemmas 17.31 17.41 and 17.71 and collect all the terms 

which involve factors of Ti^^ or To^\ Using (j7.12p . we find that the contribution to TZ 
involving factors of log |.^^| has the form 

/3=1 J [0] 

+ (deg<4) + o(|C"|-3). (8.1) 

As explained after (I7.13p . this term must vanish. According to our definition (I7.24p . we 
are free to multiply q{x) by any non-zero number if we divide v by the same number. 
Using this freedom, we may assume that the coefficients (defined by (I7.29P ) satisfy 
the relations 

a 

Y,m'pdp = l. (8.2) 

0=1 

Then the expression in (18. ip vanishes if and only if we set 

v = -^^ + m^YYA^. (8.3) 
b 

With ()7.33p . (I7.34P and ()8.2p we have an inhomogeneous system of three linearly 
independent equations for the unknowns di,...,dg. If 5 < 3, this system has no 
solutions. We thus conclude that the logarithmic poles can be compensated only if we 
have at least three generations. In the case 5 > 3, on the other hand, the system has 
solutions, and in AppendixOit is shown that these solutions can indeed be realized by 
a suitable mixing matrix q. Thus from now on, we assume that g > 3, and we choose g 
as well as v such that the conditions (|7.33p . (|7.34p . (18. 2p and (|8.3p hold. Then the 
logarithmic poles of TZ have disappeared. 

Before analyzing the remaining contributions to 7^, we must have a closer look at 
the non-causal low- and high energy contributions P'*' and P^'^ in the light cone expan- 
sion (j4.17p . These smooth contributions to the fermionic projector were disregarded 
in the formalism of the continuum limit as outlined in Chapter [5j This is justified 
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as long as singular contributions to the fermionic projector are considered. In par- 
ticular, contributions to P{x,y) involving the functions To ^\ Tq^'^ or their complex 
conjugates have poles on the light cone, and therefore smooth corrections would be 
of lower degree on the light cone, meaning that the corresponding contributions to 
the EL equations would be negligible corrections of the form (|5.9p . However, the fac- 
tors To^^ and Tq^^ only have a logarithmic pole, and after the above cancellations of the 
logarithmic poles, the remaining leading contributions are bounded functions. Thus 
smooth corrections become relevant. We conclude that it is necessary to determine 
the smooth contributions to the fermionic projector P{x, y) at the origin x = y. This 
analysis is carried out in Appendix [D] to first order in the bosonic potentials using a 
resummation technique. In what follows, we use these results and explain them. 
To introduce a convenient notation, we write the factors T^-^-^ in generalization 



of (|7n^ as 



.(1) 



iog\e\+i7rQ{e)e{e: 



4) 



W 327r3 

where the real-valued functions which may depend on the masses and the bosonic 
potentials, will be specified below. Taking the complex conjugate of (j7.12p . we get a 

Substituting these formulas into TZ, the factors log |^^| 



similar representation for T^^^ 



■[p] 

cancel each other as a consequence of (|8.2p and (j8.3p . Moreover, a short calculation 
shows that the factors in 0(^^) e(^'^) in ()7.12p also drop out. Applying Lemma [7. H the 
EL equations to degree four yield the vector equation 

jaN5-m^AaNe = JaNs (8.5) 

with A''3 as in (j7.9p and 

-,3 

c.c. (8.6) 



9 



qt: 



(0) L 



(0) {0)^(0)^(-l) 
^[0] ^[0] ^[0] 



3t: 



(0) 



7n(0)^(0)^(0)^(-l) _ 
-^[2] -^[1] -^[0] -'[O] 



mY ^ 

13=1 



(8.7) 



+ 5(.[0]-.[3])^Ti7M7^)M)-T(°) 



3T, 



29 Y' 



(0) 



29 YY 



(-1) (0) (0)^(0) _ 

^[0] ^[o]^[i]^[i] 



r(o) 
-[0] 



(0) (0)^(0)^(-l) _ 
-^[2] ^[1] -^[0] ^[0] 



mY ^ dfs 

13=1 



+ [S[2] - S[3] j 1 [0] ^ [0] I J r., - J r 



8.9) 
.10) 
.11) 



By direct inspection one verifies that the integration-by-parts rules (j5.8p do not yield 
relations between the simple fractions. In other words, the appearing simple fraction 
are all basic fractions. When evaluating weakly on the light cone ()5.7p . all basic frac- 
tions are of degree four, thus producing the same factor e^'^(i|^|)^^. Using furthermore 
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that no logarithmic divergences appear, we conclude that (|8.5p must hold if the basic 
fractions are replaced by the corresponding regularization parameters. 

In the case of more than three generations, the parameters di, . . . ,dg are not uniquely 
determined by the three conditions (|7.33p . (|7.34p and (j8.2p . Thus by choosing the 
generation mixing matrix g appropriately, we can give the quantity Yl$=i '^P ™ (|8.7p 
and (jS.lOp an arbitrary value. In particular, by making this quantity large, we can 
arrange that (jS.Sp holds for an arbitrarily small gauge field. In other words, by a 
suitable choice of the local axial transformation we can make the coupling constant as 
small as we like, so that the field equations become trivial. 

In the remaining case of three generations, the parameters dp are uniquely deter- 
mined by the conditions (|7.33p , ()7.34p and (|8.2p . A short calculation gives 

9 



mY 



TJif, d. 



13 "/3 



1 . 



/3=1 



Substituting this relation into (18.71) and (|8.10p , we obtain the condition 



Co - Cl(s[o] - S[3]) jja 



i C2Y^ + c^YY 



Cl 



2ci(s[2] - m)^Y) Aa = ^Ja, (8.12) 



where the constants cq, . . . , C3 are the four regularization parameters corresponding to 
the following basic fractions: 

11] "[2] ;"[o] "[0] 



Co 



6T, 



(0) 



c.c. 



Cl 



C2 



7^(0) 
2-3 



.(0) 



.(0) 



(-1) (0) (0)^(0) 
^[0] ^[0] ^[1] ^[1] 



c.c. 



i.l3) 



^ 2-3 [^(o)^(o)^(o)^(-i) 

^3 • -Z^l m [2] ^[0] ^[0] 



c.c. 



Since we are free to multiply (j8.12p by a non-zero constant, our field equations (|8.12p 
involve three regularization parameters. For a given regularization scheme, these pa- 
rameters can be computed to obtain numerical constants, as will be explored further 
in ^8.61 Alternatively, these three parameters can be regarded as empirical constants 
which take into account the unknown microscopic structure of space-time. Apart from 
these constants, all the quantities in (|8.12p are objects of macroscopic physics, defined 
independent of the regularization. 

It remains to determine the quantities spj , S[2] and S[3] . As in (jS.ip . we again consider 
the leading order at the origin, and thus it suffices to compute the functions S[p](x,y) 

at X = y. Let us begin with the calculation of 5^3]. Since the factors T^^-^ and T^^-^ only 
appear in (17.36p . the function sp] is obtained by computing the local axial transforma- 
tion of the fermionic projector of the vacuum. According to (I4.13p - f4.16p . the relevant 
contribution to the (3^^ Dirac sea involving the logarithmic pole and the constant term 
at the origin is 



m'r. 



3 log(m^|e'|) + co. 
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Since in (j8.12p only the differences of the functions S[p] appear, we may always disregard 
the constant cq. Computing the corresponding contribution to ()7.30p and comparing 
with (18.41). we find that srqi is the constant 



we find that s\: 

-3 J _3 



5[3] 



[3] 

1 E)^=lf^/3"l^log(^) 



1 



327r3 



log(mp 



n 



1 



/3 



327r3 ^ rni + m2 + m3 mp 



(8.14) 



where in the last step we again used the relations (|7.33p . ()7.34p and ()8.2p . The func- 
tions S[o] and sp] are more difficult to compute. Therefore, we first state the result 
and discuss it afterwards. 



Lemma 8.1. The operators s 
1 



and S[2] appearing in ()8.4p and (|8.12p have the form 



5[0]Ja 



3 

3 -32713 X] 

/3=1 

1 

327r3 m?YY 



log(m|)j, + /g,*ja)+0(^2) 
3 

^ m} ( log(m2 ) + * A,) + 0(^2 



/3=1 



.15) 



?.16) 



where the star denotes convolution, i.e. 



f[p\{x - y) h{y) d^y . 



The convolution kernels are the Fourier transforms of the distributions 



/fol (9) = 6 / {a-a ) log 



1 - (q - a^) 



da 



log 



1 — (q — Q^) 2 



m 



da . 



(8.17) 
(8.18) 



Postponing the proof of this lemma to Appendix |Dl we here merely discuss the result. 
For clarity, we first remark that the convolution operators in ()8.15p and (j8.16p can also 
be regarded as multiplication operators in momentum space, defined by 



/i 



-iqx 



■ 18p . Next, we note that the integrands 



with the functions f^^ as in (I8.17P and 
in (|8.17p and (|8.18p only have logarithmic poles, so that the integrals are finite. In 
Appendix |Dl these integrals are even computed in closed form (see Lemma ID. 21 and 
Figure [6]). Next, we point out that these integrals vanish if = 0, because then 
the logarithm in the integrand is zero. Therefore, the the convolutions by /j^j can be 
regarded as higher order corrections in to the field equations. Thus we can say 
that sjo] and sj2] composed of constant terms involving logarithms of the Dirac 
masses, correction terms /j^j taking into account the dependence on the momentum q'^ 
of the bosonic potential, and finally correction terms of higher order in the bosonic 
potential. 

The constant term in sp] can be understood from the following simple consideration 
(the argument for sp] is similar). The naive approach to determine the constant 
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contribution to the Dirac sea is to differentiate (j4.13p at o = m| to obtain 



1 

32^3 

where we again omitted the irrelevant constant cq. Taking the partial trace and com- 
paring with (j8.4p . we obtain the contribution 

^[o]^^:^ENK) + i) . (8.19) 

/3=1 

This naive guess is wrong because there is also a contribution to the fermionic projector 
of the form ~ ^^fc Ja ^[o|^ ) which when contracted with ^ yields another constant term 
which is not taken into account by the formalism of Chapter [5] (see the term (|B.14p in 
Appendix [B]) . This additional contribution cancels the summand +1 in (|8.19p . giving 
the desired constant term in (|8.15p . 

Next, it is instructive to consider the scaling behavior of the functions s^j in the 
fermion masses. To this end, we consider a joint scaling ^ Lmp of all masses. 
Sine the expressions (|8.14p . (j8.15p and (j8.16p have the same powers of the masses in 
the numerator and denominator, our scaling amounts to the replacement log(m|) — )• 
log(m|) + 21ogL. Using the specific form of the operators S[p], one easily verifies that 
the transformation of the constant terms can be described by the replacement S[p] — )• 
S[p] + 2/(327r3) logL. We conclude that for differences of these operators as appearing 
in (j8.12p . the constant terms are indeed scaling invariant. In other words, the constant 
terms in the expressions S[o] — S[3] and S[2] — S[3] depend only on quotients of the 
masses mi, m2 and m^. 

Before discussing the different correction terms in ()8.15p and ()8.16p . it is conve- 
nient to combine all the constant terms in ()8.15p . (j8.16p and (j8.12p . More precisely, 
multiplying (f8l^ by 967rVco gives the following result. 

Theorem 8.2. The EL equations to degree four on the light cone give rise to the 
condition 

[Co - /[0]*)ja - (C2 - 6/[2]*)Aa = 12^^ ^ o(^2) (g_20) 

involving the axial bosonic potential A^, the corresponding axial current ja o.Tt'd the 
axial Dirac current J^, (see (j6.22p . ()7.4p and ()7.8p ). Here the convolution kernels are 
the Fourier transforms of the distributions 
3 3 ^1 

1 — (a — 

/3=1 /3=1 
3 3 ^1 

{a — a ) — 2 



/3=1 /3=1 ''^ 

3 ^ 3 ^1 

/[2] {q) = ^^1 f[o] (?) = Y^} log 



ml 



da 



^2 

,2 , 

m2 



da 



(with the functions /j^j as defined by (j8.17p and ()8.18p ). The constants Cq and mfC2 
depend only on the regularization and on the ratios of the masses of the fermions. 

With this theorem, we have derived the desired field equations for the axial poten- 
tial ^a- They form a linear hyperbolic system of equations involving a mass term, 
with corrections in the momentum squared and of higher order in the potential. It is 
remarkable that the corrections in the momentum squared are described by explicit 
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convolutions, which do not involve any free constants. In order to make the effect of 
the convolution terms smaller, one must choose the constants Cq and C2 larger, also 
leading to a smaller coupling of the Dirac current. Thus the effect of the convolution 
terms decreases for a smaller coupling constant, but it cannot be arranged to vanish 
completely. 

We proceed by explaining and analyzing the above theorem, beginning with the 
convolution operators /[^j ( ^8.2p and the higher orders in the potential ( ^8.3p . In ^8.41 
we explain how the standard loop corrections of QFT appear in our approach. In ^8.51 
we explain why the Higgs boson does not appear in our framework. Finally, in ^8.61 
we compute the coupling constants and the bosonic rest mass for a few simple regu- 
larizations. 



8.2. Violation of Causality and the Vacuum Polarization. In this section we 
want to clarify the significance of the convolution operators in the field equations ()8.20p . 
Our first step is to bring the convolution kernels into a more suitable form. For any a > 
0, we denote by Sa the following Green's function of the Klein-Gordon equation, 

Sa[X,y) = I TT—r — e 



-ik{x-y) 



(27r)4 

1 ...2n , a Ji(v^) 



^21 



(8.22) 



where in the last step we again set ^ = y — x and computed the Fourier integral 
using the Bessel function Ji. This Green's function is obviously causal in the sense 
that it vanishes for spacelike E,- Due to the principal part, it is indeed the mean of 
the advanced and retarded Green's function; this choice has the advantage that Sa 
is symmetric, meaning that Sa{x,y) = Sa{y,x). Expanding the Bessel function in a 
power series, the square roots drop out, showing that Sa is a power series in a. In 
view of the explicit and quite convenient formula ()8.22p . it seems useful to express the 
convolution kernels in terms of Sa- This is done in the next lemma. 



Lemma 8.3. The distributions /j^j as defined by ()8.17p 

- y) 



an. 



d (I8.18P can he written as 



f[2]i^-y) 



4mi 



Sa{x,y) + 



4m| (a + 2m'n' 



da 

3" 

02 



Sa{x,y) + 



5^{x - y) 



4m^ 



da 



^.23) 



^24) 



Proof. We first compute the Fourier transform of the distribution log|l — /h\ for 
given 6 > 0. Using that lima-5.00 log |1 — /a\ = with convergence as a distribution, 
we have 



log 



1 



d_ 

da 



log 



1 



da 



h - Ma 



— a a 



Now we can compute the Fourier transform with the help of (j8.2ip . Setting b 



rni/ (a — a^), we obtain 



log 



1 — [a — a 



2^ r 



rrir. 



-iq{x-y) 



^4 — y\\ 

.^2 iSa{x,y)^ jda. 
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Srnix,z)V{z)tmiz,y) ^ . ^ .^ . ^ 

i(x,Z) V(Z) Srn[Z,y) 



Figure 2. The support of the integrand in (I8.25p . 
We finally integrate over a, interchange the orders of integration, 

1 roo fOO , TfP'. 



da {a — a'^Y / ^2 da 
J — ^ JAtuI 



2 f \ 
da (■ ■ ■ ) / da (a — a Y Q(a tt , 

^ Jo ^ a-a^J 

and compute the last integral. □ 



Qualitatively speaking, this lemma shows that that the distributions can 
be obtained by integrating the Green's function Sa over the mass parameter a and by 
subtracting a suitable counter term localized at ^ = 0. The interesting conclusion is 
that the convolution kernels /[p] (x, y) in the field equations ()8.20p are weakly causal 
in the sense that they vanish for spacelike ^. But they are not strictly causal in the 
sense that the past influences the future and also the future influences the past. 

Before discussing whether and how such a violation of causality could be observed in 
experiments, we give a simple consideration which conveys an intuitive understanding 
for how the non-causal contributions to the field equations come about. For simplicity, 
we consider the linear perturbation AP of a Dirac sea of mass m by a potential "B, 

AP{x,y) = - J z) V{z) tm{z,y) + tm{x, z) V{z) 

y)) , (8.25) 

where Sm is the Dirac Green's function and tm denotes the Dirac sea of the vacuum, 
i.e. 

tm = {i^ + 'm)T^2 and = (i^ + 771)5^2 , (8.26) 

and T^2 and 5„2 as defined by (j4.12p and (j8.2ip (for details see [HI equations (2.4) 
and (2.5)]). Let us consider the support of the integrand in (j8.25p . The Green's func- 
tion Sm vanishes outside the light cone (see (j8.22p ). whereas the distribution tm is 
non-causal (see (j4.13p ). Thus in (j8.25p we integrate over the union of the double light 
cone (meaning the interior of the light cones and their boundaries) centered at the 
points X and y; see the left of Figure [2j In the limit x — )■ y, the integral in ()8.25p will 
diverge, as becomes apparent in the poles of light cone expansion. But after subtract- 
ing these divergent contributions, we can take the limit a; — )• y to obtain a well-defined 
integral over the double light cone centered at the point x = y; see the right of Figure[2l 
Indeed, the finite contribution at the origin described by this integral corresponds pre- 
cisely to the smooth contribution to the fermionic projector as considered in ^8.11 



This consideration explains why the distributions fl^p-^{x — y) vanish for spacelike ^. 

We even see that the distributions /j^j(x,?/) are closely related to the pointwise prod- 
uct in position space of the Bessel functions appearing in the distributions T^2(x,y) 
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and Sjn2{x,y) for timelike ^. Going into more details, this argument could even be 
elaborated to an alternative method for computing the convolution kernels. However, 
for actual computations this alternative method would be less convenient than the 
resummation technique of Appendix [Dj 

One might object that the above violation of causality occurs simply because in (|8.25p 
we are working with the wrong Green's functions. Indeed, if in (j8.25p the first and 
second factors Sm were replaced by the regarded and advanced Green's function, re- 
spectively, the support of the integral would become strictly causal in the sense that z 
must lie in the causal past of x or y. However, modifying the Green's functions in this 
way is not admissible, as it would destroy the property that the Dirac sea is composed 
only of half of the solutions of the Klein-Gordon equation. More generally, the unique- 
ness of the perturbation expansion of the fermionic projector follows from a causality 
argument (see j'Hl §2.2]). Thus there is no freedom in modifying the perturbation 
expansion, and thus the above violation of causality cannot be avoided. 

The violation of causality in the field equations breaks with one of the most fun- 
damental physical principles. The immediate question is whether and how this effect 
could be verified in experiments. We conclude this section by discussing this ques- 
tion. Before beginning, we point out that the present paper is concerned with a simple 
fermion system, and one should be careful to draw physical conclusions from this over- 
simplified physical model. Also, the author has no expertise to address experimental 
issues. Nevertheless, it seems worth exploring the potential consequences of the causal- 
ity violation in a few "Gedanken experiments," just to make sure that we do not get 
immediate contradictions to physical observations. In order to be closer to everyday 
physics, let us consider what happened if we inserted the nonlocal convolution term 
into Maxwell's equations. For simplicity, we consider one Dirac wave function ^ of 
mass m. Thus dropping the mass term in ()8.20p and choosing for convenience the 
Lorentz gauge, the modified Dirac-Maxwell equations become 



(i^ + 4- m)* = , 



12^2 



.27) 



where chose the constant Cq such that without the convolution terms, the Maxwell 
equations take the familiar form —0^4^ = e^^'7fe^' (note that we again use the con- 
vention where the Dirac equation involves no coupling constants; see also Footnote [1] 
on page [TSl) . In view of Lemma 18.31 the square bracket is an integral operator which 
vanishes for spacelike distances. Furthermore, we see from ()8.17p and (j8.18p (for more 
details see Lemma lD.2p that the functions f^ig) diverge for large q"^ only logarithmi- 
cally. Thus in view of the smallness of the fine structure constant e^/47r ~ 1/137, for 
the energy range accessible by experiments the square bracket in ()8.27p is an invertible 
operator. Thus we may write our modified Dirac-Maxwell equations as 



+ = , 



-aA, 



12^2 



-1 



showing that the convolution term can be regarded as a modification of the source 
term. Alternatively, one may write the Maxwell equation in the standard form 



(8.29) 
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with a so-cahed auxiliary potential A and take the point of view that the convolution 
term only affects the coupling of the electromagnetic potential to the Dirac equation, 



+ = with A 



e2 



i-Y^/[o] 



-1 



A (8.30) 



(note that the wave and convolution operators commute, as they are both multiplica- 
tion operators in momentum space). Both the "source form" (j8.28p and the "coupling 
form" ()8.29p and ()8.30p are useful; they give different points of view on the same sys- 
tem of equations. We point out that, as the inverse of a causal operator, the operator 
on the right of ()8.28p and (j8.30p is again causal in the sense that its integral kernel 
vanishes for spacelike distances. Moreover, for large timelike distances the kernel /[q] 
is oscillatory and decays. More specifically, writing the Green's function Sa in ()8.23p 
with Bessel functions and using their asymptotic expansion for large one finds that 



/gj(x -y)^m^ (e^)-! cos [^4m}e + V') if » m} (8.31) 

(where (p is an irrelevant phase). 

The formulation (j8.29p and (|8.30p reveals that our modified Dirac-Maxwell equations 
are of variational form. More precisely, they can be recovered as the EL equations 
corresponding to the modified Dirac-Maxwell action 

o2 



5dm = / i^{i^ + 
Im 



1 



m 



where F is the field tensor corresponding to the auxiliary potential. Hence by applying 
Noether's theorem, we obtain corresponding conserved quantities, in particular the 
total electric charge and the total energy of the system. Thus all conservation laws 
of the classical Dirac-Maxwell system still hold, but clearly the form of the conserved 
quantities must be modified by suitable convolution terms. 

The simplest idea for detecting the convolution term is to expose an electron to 
a laser pulse. Then the convolution term in the Dirac equation (|8.30p might seem 
to imply that the electron should "feel" the electromagnetic wave at a distance, or 
could even be infiuenced by a laser beam fiying by in the future, at a time when the 
electron may already have moved away. However, such obvious violations of causality 
are impossible for the following reason: An electromagnetic wave satisfies the vacuum 
Maxwell equations OA = (see (I8.29P ). Thus the momentum squared of the elec- 
tromagnetic wave vanishes, implying that /[q] * A = 0, so that the convolution term 
in (j8.30p drops out. In more general terms, the convolution terms are constant if the 
bosonic field is on-shell. We conclude that the convolution terms can be detected only 
by off-shell bosonic fields, which according to (I8.29P occur only at the electromagnetic 
sources. 

Another idea for observing the convolution term is that, according to (I8.28p . it 
modifies the way the Dirac current generates an electromagnetic field. Due to the 
prefactor e^/(127r^) and in view of the fact that the kernel f^ decays and has an 
oscillatory behavior (I8.31j) . this effect will not be large, but it could nevertheless be 
observable. In particular, one may ask whether the positive and negative charges of 
protons and electrons still compensate each other in such a way that macroscopic ob- 
jects appear neutral. If this were not the case, this would have drastic consequences, 
because then the electromagnetic forces would dominate gravity on the large scale. To 
analyze this question we consider for example a crystal containing exactly as many 
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positive and negative charges. Then the corresponding auxihary potential A vanishes 
outside the crystal (except for dipole effects, which fall off fast with increasing dis- 
tance). As a consequence, the potential A defined by (j8.3Up also vanishes outside the 
crystal, and thus there are no observable electrostatic forces outside the crystal, in 
agreement with physical observations. 

More generally, the above considerations show that the convolution term can lead to 
observable effects only if the sources of the electromagnetic field and the Dirac particles 
on which it acts are very close to each other, meaning that the whole interaction must 
take place on the scale of the Compton length of the electron. One conceivable way 
of measuring this effect is by considering electron- electron scattering. In order to 
concentrate on the violation of causality, it seems preferable to avoid the noise of 
the usual electromagnetic interactions by considering two wave packets which stay 
causally separated, but nevertheless come as close as the Compton length. In this 
case, an electron in the future could even affect the motion of an electron in the past. 
However, due to the Heisenberg uncertainty principle, localizing a wave packet on the 
Compton scale implies that the energy uncertainty is of the order of the rest mass, so 
that pair creation becomes a relevant effect. Therefore, arranging such wave packets 
seems a very difficult task. 

Another potential method for observing the convolution term is to get a connection 
to the high-precision measurements of atomic spectra. Thus we conclude the discussion 
by considering the static situation. Integrating the Green's function (j8.2ip over time, 
we can compute the remaining spatial Fourier integral with residues to obtain the 
familiar Yukawa potential, 

dk e-^^^ 



(2^)2 70 k^ + a 



3 (27r)3 |fc|2 + a 



-1 



h 1 f,-\fa,r 



(27r)2 ir j k'^ + a 47r r ' 

where we set r = |^|. Hence in the static case, the convolution operator reduces to the 
three-dimensional integral 

(/[o] * h) {x) = / /[o] {x - y) h{y) dy 

involving the kernel 



1 . , 



da 



a - (a + 2m^) ^ . (8.32) 
a2 



/3=r 

We now consider a classical point charge Ze located at the origin. In order to com- 
pute the corresponding electric field we consider the corresponding Maxwell equa- 
tion (18:28]) . 



AAo(f 



e2 



l-I^/[0] 



Ze' 5'{x) = Ze' 6-\x) + /[o] (f) + 0{e') . 



In order to solve for Aq, we convolute both sides with the Newtonian potential Vb(^) = 
— l/(47rr). To compute the resulting convolution of the Newtonian potential with /[g]. 
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we first observe that (j8.32p involves the Yukawa potential Va- Since convolution cor- 
responds to multiplication in momentum space, we can use the simple transformation 

1 1 if 1 1 \ 

+ — 



to conclude that 



We thus obtain 



|A;|2 |A;|2 + a « 1 |A:|2 + a \k\^ 



{Vo*Va){x) = ^{Va{x)-Vo{x)) 

a 



Ze2 Ze^ ^ 



-v^\i\ 1 



Airr a 



Ami (a + 2ml 



da 

3 

a2 



Here the first summand is the Coulomb potential, whereas the second summand is an 
additional short-range potential. This is very similar to the situation for the relativistic 
correction described by the Darwin term (a relativistic correction to the Schrodinger 
equation; see for example [l2l Section 3.3]). Concentrating the short range potential 
at the origin by the replacement 

6Hx) , 



47rr a 

we can carry out the a-integral to obtain 

. Ze^ Ze^ 1 rS/^N 

Aq{x) = 7T > 5 (x) . 

^ Anr 607r2 ^ ml ^ ^ 

13=1 P 

We thus end up with a correction to the Dirac Hamiltonian of the form 

Ze^ ^ 1 

-f^noncausal = o / ^ — T 5^(x) . (8.33) 

bUvr^ nip 

This coincides precisely with the Uehling potential, which describes the one-loop 
vacuum polarization (see for example [391 eqns (7.94) and (7.95)]). At first sight, it 
might be surprising that we get the same result as in perturbative QFT, although we 
did not consider a fermionic loop diagram and did not encounter the usual ultraviolet 
divergences. In order to see the connection, it is preferable to reconsider the origi- 
nal derivation by Uehling and Serber [47\ I45j . which was based on earlier papers by 
Dirac [8] and Heisenberg [32]. Similar to p.7p . Dirac considers the sum over all sea 
states, 

R{t,x;t',x') = Yl '^i{t,x)Wt\¥) , 

I occupied 

where the wave functions are solutions of the Dirac equation 

(i^ + e4{t, x) - m) {t, x) = 0. 

Thus up to an irrelevant overall constant, R coincides precisely with the kernel of 
the fermionic projector P{x,y). Dirac realizes that R has singularities on the light 
cone and discusses their form. Heisenberg pushes the calculation a bit further and, 
using physical conservation laws and the requirement of gauge invariance, he brings 
the singular contribution to R into a canonical form. The he argues that this singular 
contribution to R should simply be dropped. Uehling and Serber took the remaining 
regular contribution to R as the starting point to derive the corresponding correction 
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to the Maxwell equations. Since the regular contribution to R coincides precisely 
with the non-causal contributions in (|4.17p (albeit in a less explicit form where the 
underlying causal structure is not apparent), it is clear that (|8.33p coincides with the 
usual Uehling potential. 

We conclude that the non-causal correction reproduces the usual vacuum polariza- 
tion as described by the Uehling potential. The main difference in our approach is that 
the singular contributions to the fermionic projector are not disregarded or removed, 
but they are carried along in our analysis. These singular contributions then drop out 
of the Euler-Lagrange equations corresponding to our action principle (j2.9p . In this 
way, all divergences disappear. The remaining finite contributions to the fermionic 
projector give rise to the bosonic current and mass terms in the resulting field equa- 
tions (|8.20p . and also yield the convolution terms which describe the vacuum polar- 
ization. The main advantage of the fermionic projector approach is that no counter 
terms are needed. The back-reaction of the Dirac sea on the electromagnetic field is 
finite, no divergences occur. Moreover, as we do not need counter terms, the setting 
immediately becomes background independent. It is to be expected (although it has 
not yet been worked out in detail) that the singularities of the fermionic projector 
will also drop out of the Euler-Lagrange equations if one sets up the theory in curved 
space-time. 

We finally remark that the connection to Feynman diagrams will be explained in 
more detail in Section 18.41 

8.3. Higher Order Non-Causal Corrections to the Field Equations. The non- 
causal convolution terms in the previous section were obtained by computing the non- 
causal contributions in (I4.17P at the origin, considering the first order of the pertur- 
bation expansion (|4.9p . Likewise, the higher orders of this expansion also contribute 
to P^"" and P^"", giving rise to higher order non-causal corrections to the field equations. 
In this section we briefly discuss the structure of these correction terms (computing 
them in detail goes beyond the scope of this paper). 

It is natural to distinguish between the low and high energy contributions. The 
non-causal low energy contribution P'*^ in (j4.17p can be computed at the origin to 
every order in 23 by extending the resummation technique of Appendix [D] to higher 
order (more precisely, according to the residual argument, we again get sums of the 
form ()D.12p . but with nested line integrals and multiple series, which are to be carried 
out iteratively) . Similar as explained to first order after Lemma 18.31 the resulting cor- 
rections to the field equation are weakly causal in the sense that they can be described 
by convolutions with integral kernels which vanish for spacelike distances. Thus they 
have the same mathematical structure, but are clearly much smaller than the convo- 
lution terms in ^8.2[ 

The non-causal high energy contribution P^'^ in ()4.17p is more interesting, because it 
gives rise to corrections of different type. For simplicity, we explain their mathematical 
structure only in the case of one generation and only for the leading contribution to P^'^ 
(see [22] for details) 
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{i^ + m) T„2 



and tm = {i(p + fn) T^p- , 
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and Ta is the complex conjugate of the distribution Ta, ()4.12p . Thus the distribu- 
tions tjn and tm are supported on the lower and upper mass shell, respectively. Eval- 
uating this expression at the origin gives 

P^''{X,X) = [ d^Zi [ (fz2{tm{x,Zi)'B{zi)t^{zi,Z2)'B{z2)tm{z2,x) 

4 Jm Jm ^ (8.34) 

-t;;i{x,Zi)'B{zi)tm{zi,Z2)'B{z2)t;;i{z2,x)^+0{'B'^) . 

This is similar to a second order tree diagram, but instead of Green's functions it in- 
volves the projectors onto the lower and upper mass shells, which appear in alternating 
order. The expression is well-defined and finite (see [14| Lemma 2.2.2]). Similar to the 
correction terms in Theorem 18.21 our expression is a convolution, but now it involves 
two integrals, each of which contains one factor of 23. Consequently, the integral kernel 
depends on two arguments zi and Z2- The interesting point is that this integral kernel 
does not vanish even if the vectors zi — xotz2 — x are space-like. Thus the correspond- 
ing corrections to the field equations violate causality even in the strong sense that 
in addition to an influence of the future on the past, there are even interactions for 
spacelike distances. This surprising result is in sharp contrast to conventional physical 
theories. However, since for space-like separation the kernels tm decay exponentially 
fast on the Compton scale, the effect is extremely small. In particular, describing this 
exponential decay by the the Yukawa potential, this effect could be described similar to 
the correction (I8.32p and (IS.SSp . But compared to the latter first order correction, the 
second order correction by P^'^ would be smaller by a factor . In view of the discus- 
sion in §8.21 measuring this correction is at present out of reach. Thus it seems that the 
only promising approach for detecting an effect of the high energy contribution is to 
look for an experiment which is sensitive to interactions between regions of space-time 
with spacelike separation, without being disturbed by any causal interactions. 

8.4. The Standard Quantum Corrections to the Field Equations. We now 

explain how the quantum corrections due to the Feynman loop diagrams arise in our 
model. We will recover all the standard quantum corrections. Moreover, we will obtain 
quantum corrections of the previously described non-causal terms (see §8.21 and §8.3p . 
For clarity, we proceed in several steps and begin by leaving out the non-causal con- 
volution terms in the field equations (I8.20p . Furthermore, we consider only one Dirac 
particle of mass m and disregard the interaction of this particle with the states of the 
Dirac sea. Under these simplifying assumptions, the interaction is described by the 
coupled Dirac- Yang/Mills equations 

[i^ + -m)^> = Q, dkiA' - UAk - M^Ak = e^ ^T^Tfc* , (8-35) 

where A is the axial potential, and the bosonic rest mass M and the coupling constant e 
are determined from (I8.20p by setting = C2/C0 and = 127r^/Co. We point out 
that the wave function ^ and the bosonic field A in (I8.35P are classical in the sense 
that no second quantization has been performed. 

The equations (I8.35P form a coupled system of nonlinear hyperbolic partial differ- 
ential equations. For such a system, one can get local existence and uniqueness results 
(see for example [33 j Section 5.3] or j46l Chapter 16]), but constructing global solutions 
is a very difficult task. Therefore, we must rely on a perturbative treatment, giving 
a connection to Feynman diagrams. Although this connection is quite elementary, 
it does not seem to be well-known to mathematicians working on partial differential 
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equations. In physics, on the other hand, Feynman diagrams are usually derived from 
second quantized fields, where the connection to nonlinear partial differential equa- 
tions is no longer apparent. Therefore, we now explain the procedure schematically 
from the basics, hopefully bridging a gap between the mathematics and physics com- 
munities. In order to be in a simpler and more familiar setting, we consider instead 
of (j8.35p the Dirac-Maxwell equations in the Lorentz gauge, as considered in quantum 
electrodynamics (see for example 

{i^ + e4 - m)"^ = , -□^fc = e^7fe^'. (8.36) 

The natural question in the theory of hyperbolic partial differential equations is the 
Cauchy problem, where we seek for solutions of (j8.36p for given initial values 

^{t,x)\t=o = ^o{x) , A{t,x)\t=o = Ao{x) , dtA{t,x)\t=o = Ai{x) . (8.37) 

In preparation, we formulate the equations as a system which is of first order in time. 
To this end, we introduce the field $ with components 

/ ^ \ 

$ = A , (8.38) 

\idtAj 

and write the system in the Hamiltonian form 

idt^{t, x) = H{^{t, x)) := Ho<i> + eB{^) , (8.39) 

where in the last step we decomposed the Hamiltonian into its linear and non-linear 
parts given by 

Ho= I 1, B{<^) = _0 . (8.40) 

V -AO/ V ^'7^' / 

In the case e = 0, we have a linear equation, which is immediately solved formally by 
exponentiation, 

$(t) = e-**^«$o , 

where we set $o = *^|t=o- This equation is given a rigorous meaning by writing the 
so-called time evolution operator e"**^" as a spatial integral operator. 

Lemma 8.4. For any t > 0, the operator e"**^" can be written as 



(e-**^°$)(x) = / Rt{x-y) ^y) dy , (8.41) 



^ In order to bring the system (|8.35p into a comparable form, one first takes the divergence of the 
Yang/Mills equation to obtain 

where in the last step we used the Dirac equation. In particular, the divergence of A in general does 
not vanish. It seems convenient to subtract from A the gradient of a scalar field $, 

Bk ■- Ak - 9fe$ , 

in such a way that the new potential B becomes divergence- free. This leads to the system of equations 
(j^-m + 7^^+7'O))* = 0, = -=^^'*7'^* , (-□-M^)Bfe = e'*7Sfc* + M^afc-I>. 

This system has the same structure as (|8.36p . and it can be analyzed with exactly the same methods. 
For the handling of the factors e see Footnote [T] on page 1151 
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where the integral kernel is the distribution 

/.^(t,x)(z7°) \ 

Rt(.x)=\ -dtS(^{t,x) iS(^{t,x) \, (8.42) 

\ -iAS^{t,x) -dtS^{t,x)J 

which involves the retarded Green's functions defined by 

/Mh. -ikx 
-J—- -r (8.43) 
(27r)4 k^-m? + iekO ^ ^ 



stit,x) = {^l^^ + m)S^,ix). (8.44) 

Proof. Using that for any t > 0, the Green's function S"^ is a solution of the Klein- 
Gordon equation (— □ — a)S'a(x) = 0, a short calculation using (|8.44p shows that (|8.4ip 
is a solution of the equation {idt — i/o)(e~**^"*^*) = 0. In order to verify the correct ini- 
tial conditions, we differentiate with respect to time and carry out the t-integration 
with residues to obtain 

1 f ^ f°° (—ii A" 



limd?S^(t,x) = ——J dke"'^' lim / e"*^' 

t\0 ^ (27r)*7R3 e,t\oJ_ 

— + ieuj 

- J— f rlh , ^ .^ (-^^)" _ / if n = 

" (27r)4 7^3 ^ ' 2uj a.=±^^l^ " 1 -53(f) ifn = l. 

Using this result in ()8.44p and (|8.42p shows that indeed \\m.t\j) Rt{x) = 5^{x). □ 

In the nonlinear situation e 7^ 0, it is useful to work in the so-called "interaction 
picture" (see for example [SI Section 8.5]). We thus employ the ansatz 

m = e-'*^«1>i„t(t) , (8.45) 
giving rise to the nonlinear equation 

idt^,^t = e5int($int(t)) , (8.46) 

where 

Ant(fint(t)) = e^*^«i?(e-^*^oci>.^,(t)). 
We regard ()8.46p as an ordinary differential equation in time, which in view of (j8.4ip is 
nonlocal in space. From (|8.45p one sees that <&int comes with the initial data <I>int|t=o = 
$0- Taking a power ansatz in e, 

$int(t) = + e$S(t) + ^^{t) + . . . , 

a formal solution of the Cauchy problem for $int is obtained by integrating ()8.46p 
inductively order by order, 

it) = <f , ^[11 {t) = -^ 1^ Ant (r)) dr 
^SW = VA„t($L°tV)) • ^S(r) dr 

= {-if dr Vi?int {<^[^ (r)) £ da B,^, {a, 

(here Vi? denotes the Jacobi matrix of B, where as in (I5.15P we consider the real 
and imaginary parts of the arguments as independent variables). Substituting these 
formulas into (j8.45p . we obtain the desired solution <^ of the original Cauchy problem 
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Figure 3. Feynman tree diagrams (left) and bosonic loop diagrams (right). 

expressed as a sum of iterated time integrals, involving intermediate factors of the free 
time evolution operator e~*^^°. In particular, we obtain to second order 



We remark that in the case when B{^) is linear in $, this expansion simplifies to 
the well-known Dyson series (also referred to as the time-ordered exponential). In 
view of (j8.38p . we have derived a unique formal solution of the Cauchy problem ()8.36p 



Combining the above expansion for <I>(t) with the formula for the time evolution op- 
erator in Lemma 18.4^ one can write the above perturbation expansion in a manifestly 
covariant form. Namely, when multiplying the operators Rt with B (or similarly 
or higher derivatives) , the factors 7" in the first component of (I8.42p and in the formula 
for B in (I8.40p cancel each other, giving the Lorentz invariant expression s^^. Likewise, 
the Dirac current in (I8.40p multiplies the retarded Green's function Sq . Moreover, we 
can combine the spatial and time time integrals to integrals over Minkowski space. 
In this way, we can identify the contributions to the perturbation expansion with the 
familiar Feynman diagrams. More precisely, every integration variable corresponds 
to a vertex of the diagram, whereas the bosonic and fermionic Green's functions Sq 
and are written as wiggled and straight lines, respectively. Denoting the argument 
of the solution $(t,y) by y, whereas x = (0,x) stands for the argument of the initial 
values, we obtain different kinds of diagrams. First, one gets all the tree diagrams as 
exemplified in Figured (left). However, the bosonic loop diagrams are missing as a 
consequence of the strict time ordering in the solution of the Cauchy problem. But 
if one relaxes the time-ordering and replaces the retarded Green's functions by more 
general Green's functions (as is motivated physically by the fact that the elecromag- 
netic "cloud" generated by the self-interaction surrounds the electron already at initial 
time), then one gets all bosonic loop diagrams (as for example in Figure [3] (right)). 
We come to the following conclusion: 

• All bosonic loop diagrams can be obtained from the nonlinear system of par- 
tial differential equations (]8.36p . working purely with classical fields. For the 
derivation of the bosonic loop diagrams, there is no need for second quantiza- 
tion. 




and ([8:371) . 
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In order to make the connection to QFT clearer, we point out that in quantum 
physics one usually does not consider the initial value problem (j8.37p . Instead, one is 
interested in the n-point functions, which give information about the correlation of the 
fields at different space-time points. The two-point function is obtained by choosing 
initial values involving ^^-distributions. Similarly, all the n-point functions can be 
recovered once the solution of the Cauchy problem is known. Thus from a conceptual 
point of view, the only difference between our expansion and the Feynman diagrams 
in QFT is that, since in quantum physics the Feynman diagrams do not come from 
an initial value problem, there is a freedom in choosing the Green's function. Note 
that in the setting of the Cauchy problem, one necessarily gets the retarded Green's 
function (see (|8.42p ). In contrast, in QFT one is free to work instead with any other 
Green's function. Indeed, different choices lead to different approaches for handling 
the perturbation series. The most common choice is the so-called Feynman propagator 
(see for example [2]), which is motivated from the physical picture that the positive 
frequencies (describing particles) move to the future, whereas the negative frequencies 
(corresponding to anti-particles) move to the past. In this standard approach, the 
loop diagrams diverge. This problem is bypassed in the renormalization program by 
first regularizing the diagrams, and then removing the regularization while simulta- 
neously adjusting the masses and coupling constants (see for example [1]). A QFT is 
called renormalizahle if this renormalization procedure works to all orders in perturba- 
tion theory, involving only a finite number of effective constants. There are different 
equivalent renormalization procedures, the most common being dimensional renormal- 
ization (see for example |39j). But the Feynman propagator is not a canonical choice, 
and indeed this choice suffers from the problem of not being invariant under general 
coordinate transformations (for more details see |14l §2.1]). An alternative method, 
which seems natural but has not yet been worked out, would be to extend the choice of 
Green's functions in the causal perturbation expansion (14. 9p (see also |22j ) to the loop 
diagrams. Yet another method is the so-called causal approach based on ideas of Ep- 
stein and Glaser [9], which uses the freedom in choosing the Green's function to avoid 
the divergences of QFT (see also 03]). We also mention that our above derivation 
of Feynman diagrams is certainly not the most sophisticated or most elegant method. 
Maybe the cleanest method for the formal perturbation expansion is obtained in the 
framework of path integrals (see for example |371 140]). 

Recall that one simplification of the system ()8.35p was that we considered only 
one Dirac particle and disregarded the interaction of this particle with the states of 
the Dirac sea. In particular, we did not allow for the creation of a particle/anti- 
particle pair. This shortcoming is reflected in our perturbation expansion in that the 
fermionic loop diagrams are missing (see Figure [3] for a few examples). However, we 
already encountered the loop diagram in Figure H] (left) when discussing the vacuum 
polarization in ^8.2[ This shows that the fermion loops appear once we take into 
account the convolution terms in ()8.20p and the other non-causal corrections mentioned 
in ^8.31 In order to include these corrections in the above expansion, we simply add 
them to the perturbation operator eB in ()8.39p . We conclude that 

• The framework of the fermionic projector in the continuum limit yields all the 
Feynman diagrams of standard QFT, including all bosonic and fermionic loop 
diagrams. 

We point out that by this statement we do not mean that the framework of the 
fermionic projector is equivalent to perturbative QFT. As a major difference, the 
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Figure 4. Typical Feynman diagrams involving fermion loops. 



perturbation expansion of the fermionic projector involves a non-trivial combinatorics 
of operator products involving different types of Green's functions and fundamental 
solutions (see |22j for details, or see ()8.34p for an example of a contribution which 
is absent in standard QFT). This difference has no influence on the singularities of 
the resulting Feynman diagrams, and thus we expect that the renormalizability of the 
theory is not affected. But the higher-loop radiative corrections should depend on the 
detailed combinatorics, giving the hope to obtain small deviations from standard QFT 
which might be tested experimentally. The systematic and detailed analysis of these 
corrections will be carried out in a forthcoming publication. 

We close this section with two remarks. First, we again point out that in the above 
derivation of the Feynman diagrams, we worked with classical bosonic fields. This 
raises the question why a quantization of the bosonic fields is at all needed, and what 
this "quantization" actually means. Here we shall not enter a discussion of this point, 
but refer the reader to [15\ Section 4] and to the constructions in [19]. The second re- 
mark is that the described method of first taking the continuum limit, then expanding 
the resulting equations in terms of Feynman diagrams and renormalizing these dia- 
grams should be considered as preliminary. This method has the great advantage that 
it gives a simple connection to Feynman diagrams and to the renormalization program, 
making it easier to compare our approach to standard QFT. But ultimately, a fully 
convincing theory should work exclusively with the regularized fermionic projector, 
thereby completely avoiding the ultraviolet divergences of QFT. Before one can attack 
this program, one needs to have a better understanding of our variational principle in 
discrete space-time. 

8.5. The Absence of the Higgs Boson. In this section we compare the mechanism 
leading to the mass term in the field equations (see (I7.7P and (I8.20j) ) with the Higgs 
mechanism of the standard model. Clearly, our framework is considerably different 
from that of the standard model, so that no simple comparison is possible. But in order 
to make the connection between the formalisms as close as possible, we can consider the 
effective action of the continuum limit and compare it to the action of a corresponding 
model involving a Higgs field. For simplicity leaving out the non-causal convolution 
terms and considering only one particle of mass m, the field equations ()8.20p coupled 
to the Dirac equation are recovered as the EL equations corresponding to the action 
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where A denotes the axial potential and Fij = diAj — djAi is the corresponding field 
tensor. The coupling constant e and the bosonic mass M are related to the constants 
in ([8:20]) by = 12tt'^/Co and = C2/CQ. We point out that this action is not 
invariant under the axial gauge transformation 

^(x) ^ e-*^''^(^)^(x) , A^A + dK, (8.48) 

because both the fermionic mass term m'^^ and the bosonic mass term M'^AjA^ /(^e^) 
have no axial symmetry. As explained in §6.21 the absence of an axial symmetry can 
be understood from the fact that the transformation of the wave function in (|8.48p is 
not unitary, and thus it does not correspond to a local symmetry of our functionals 
in (I2SD (see also dOi]) and ([635])). 

The axial gauge transformation (|8.48p can be realized as a local symmetry by adding 
a Higgs field in complete analogy to the procedure in the standard model. More 
precisely, we introduce (/> as a complex scalar field which behaves under axial gauge 
transformations as 

(t){x) ^ e-2^^(^)</.(x) . (8.49) 

The fermionic mass term can be made gauge invariant by inserting suitable factors 
of (f). Moreover, in view of (|8.48p . we can introduce a corresponding gauge-covariant 
derivative D by 

Dj = dj + 2iAj . 

Thus the Dirac- Yang/Mills-Higgs action defined by 

is invariant under the axial gauge transformation (I8.48P and (I8.49p . We now follow 
the construction of spontaneous symmetry breaking in the standard model. For V we 
choose a double well potential having its minimum at = 1. Then the Higgs field (p 
has a non-trivial vacuum with |(/>| = 1. Thus choosing an axial gauge where (j) is real 
and positive, we can write cp as 

(f){x) = 1 + h{x) 

with a real-valued field h. Since h vanishes in the vacuum, we may expand the action 
in powers of h. Taking the leading orders in /i, we obtain the action after spontaneous 
symmetry breaking 

^DYMH = 5dYM + ^Higgs 

with 

^Higgs = ^ |-m/i^^ + ^ AjA^ + ^ {djh){&h) - 2V"{1) /i^l d*x . (8.51) 



We conclude that for the action (I8.50p . the Higgs mechanism yields an action which 
reproduces the effective action of the continuum limit (j8.47p , but gives rise to additional 
terms involving a real Higgs boson h. The Higgs boson has a rest mass as determined 
by the free parameter V"{1). It couples to both the wave function ^ and the axial 
potential A. 
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The Higgs field h can also be described in the setting of the fermionic projector, as 
we know explain. Note that the coupling terms of the Higgs field to ^ and A can be 
obtained from ()8.47p by varying the masses according to 

m {1 + h{x)) m , M {1 + h{x)) M . 

Taking into account that in our framework, the bosonic masses are given in terms of 
the fermion masses (see (j8.12p ). this variation is described simply by inserting a scalar 
potential into the Dirac equation. Likewise, for a system involving several generations, 
we must scale all fermion masses by a factor 1 + h. This is implemented in the auxiliary 
Dirac equation (j4.5p by choosing 



■B = -mh{x)Y. (8.52) 

The remaining question is whether scalar perturbations of the form (j8.52p occur in 
the setting of the fermionic projector, and whether our action principle (j2.9p repro- 
duces the dynamics of the Higgs field as described by (|8.5ip . A-priori, any symmetric 
perturbation of the Dirac equation is admissible, and thus we can certainly consider 
the scalar perturbation (|8.52p . Since (|8.52p is even under parity transformations, the 
corresponding leading perturbations of the eigenvalues and will be the same. In 
view of (17. ip and the formulas for the unperturbed eigenvalues (I6.27p . we find that the 
leading perturbations corresponding to (|8.52p drop out of the EL equations. We thus 
conclude that, although a Higgs field can be described in our framework, the action 
principle (12. 9p does not describe a dynamics of this field, but instead predicts that 
the Higgs field must vanish identically (for more details on scalar perturbations see 
Lemma IB.ip . 

8.6. The Coupling Constant and the Bosonic Mass in Examples. The reg- 
ularization parameters co,...,C3 in the field equations (|8.12p are given in terms of 
simple fractions (j8.13p . For a given regularization method, we can evaluate these sim- 
ple fractions and compute the coupling constant and the bosonic rest mass. We now 
exemplify the procedure by considering the two simplest methods of regularization: 

(A) An exponential factor in momentum space: We define the distribution in (13. Sp 

by 

9 

p^{k) = ^(^ + mp) 6{k^ - m^) Q{-k^) exp(eA:°) . (8.53) 

/9=1 

To the considered leading degree on the light cone, this regularization corresponds to 
the simple replacements 

7^(0) , 1 1 T^-^) . _^^t(°) = ^ (SKA) 

[Pl 87r3 2r{t-r- ie) ' [p] ^ r dt W 87r3 {t - r - ief ' ^ ^ 

and similarly for the complex conjugates. Using these formulas in (|8.12p . the basic 
fractions all coincide up to constants, giving the equation 

- ^ - 967r3 (s[o] - S[3])) + Irr? (f^ _ yy + 967r3(s[2] - s^^^YY) A = 12^^ _ 

According to Lemma 18.11 and ()8.14p . the functions involve the masses of the 
fermions and also the convolution terms /j^j. For clarity, we here leave out the con- 
volution terms (which are analyzed in detail in ^8.21 and Appendix |D]) . Then we can 
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Figure 5. The coupling constant and the bosonic mass for the regularization (j8.54p . 
write the field equation in the usual form 

where the coupling constant e and the mass M are given functions of the ratios 7112 /rrii 
and m^/mi of the fermion masses. In Figure El these constants are plotted as functions 
of the mass ratios. The coupling constant is of the order one; it is largest if the fermion 
masses are close to each other. The term is always positive, so that the mass term 
has the correct sign. The bosonic mass has the same order of magnitude as the fermion 
masses. 

(B) A cutoff in momentum space: In (A) we considered a regularization which was 
"soft" in the sense that it was smooth in momentum space. To give a complementary 
example, we now consider the "hard" regularization obtained by inserting a Heaviside 
function into the integrand of (|3.ip . Thus we define the distribution in ()3.3p by 

a 

P^k) = + mp) Sik^ - w?) Q{-k^) 9(1 + ek^) . 

13=1 

This regularization is described in analogy to (I8.54p by 

[P] 167r3 2r {t - r) ' W r dt W 

Thus as expected, the cutoff in momentum space gives rise to rapid oscillations in 
position space. Using these formulas in (j8.12p . the resulting basic fractions are no 
longer multiples of each other. But the weak evaluation integrals ()5.7p can still be 
computed in closed form. We thus obtain the field equation 

( - ^ - 96vr3 (s[o] - S[3]))ia + (y^ -YY + 64^3(s[2] - S[3])yy) A = 12^' Ja • 

This equation can be analyzed exactly as in example (A), giving the same qualitative 
results. It is remarkable that the constants in the field equations in example (A) 
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and (B) differ at most by a factor 3/2, indicating that our results do not depend 
sensitively on the method of regularization. 

9. The Euler-Lagrange Equations to Degree Three and Lower 

The wave functions in ()4.10p do not only have a vector and axial component as con- 
sidered in m.2\ but they also have scalar, pseudoscalar and bilinear components. We 
will now analyze the effect of these contributions on the EL equations ( ^9.11 and ^9.2p . 
Moreover, we will insert further potentials into the Dirac equation and analyze the 
consequences. More precisely, in ^9.21 we consider bilinear potentials, whereas in ^9.31 
we consider scalar and pseudoscalar potentials and discuss the remaining possibilities 
in choosing other potentials and fields. We conclude this chapter with a discussion of 
the structure of the EL equations to degree three and lower ( ^9.4|) . 

9.1. Scalar and Pseudoscalar Currents. In analogy to the Dirac currents in m.2\ 
we introduce the scalar Dirac current Jg and the pseudoscalar Dirac current Jp by 

Js = Y,^'^k-^^l'^l and Jp = J2^^^^^k-^^lil^'^l- (9.1) 

fc=l 1=1 k=l 1=1 

According to (j4.10p . these currents lead to a perturbation of the fermionic projector. 
In view of the fact that the scalar and pseudoscalar currents involve no Dirac matrix 
which could be contracted with a factor of ^, one expects that the resulting contribution 
to the EL equations should be one degree lower than that of the axial current (see 
Lemma l7.4p . Thus to leading order at the origin, one might expect contributions of 
the form 

Q{x,y) X Jp (monomial of degree three) + (deg < 3) . (9-2) 

However, perturbing the fermionic projector of the vacuum by (j9.ip . one sees that the 
corresponding contribution to Q{x,y) of degree three on the light cone vanishes (see 
Lemma IB.H where it also explained how the cancellations come about). Taking into 
account the axial potentials, we do get contributions to Q{x,y) of degree three on the 
light cone, which are of the form 

Q{x, y) X A^^ Js (monomial of degree three) + (deg < 3) . (9-3) 

However, these contribution have the same tensor structure as the contributions of the 
axial potentials and currents of degree three, which will be discussed in §9.41 below. 

9.2. Bilinear Currents and Potentials. The particles and anti-particles in ()4.10p 
also have a bilinear component, leading us to introduce the bilinear Dirac current Jt, 

by 

ria 

k=l 1=1 

(here o"-^^ = §[7-^,7'^] are again the bilinear covariants). Likewise, one may want to 
insert a bilinear potential 

■B = Hij{x)a'^ 

into the auxiliary Dirac equation (|4.5p (where H is an anti-symmetric tensor field). 
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Let us briefly discuss the effect of a bilinear current and a bilinear potential on the 
EL equations. The bilinear current corresponds to a perturbation of the fermionic 
projector of the form 

AP(x,y) = -i^J^a,, + o(|eT). 

The bilinear potential, on the other hand, gives rise to different types of contributions 
to the light-cone expansion which involve H and its partial derivatives (for details 

see |12l Appendix A. 5]). When computing the perturbations of the eigenvalues A^''^ 
of the closed chain (cf. [I4i Appendix G] or Appendix [B]), all Dirac matrices in AP 
are contracted with outer factors ^. As a consequence, the contribution of the bi- 
linear current drops out. Moreover, due to the anti-symmetry of H, all the bilinear 
contributions of the bilinear potential to the EL equations vanish. What remains are 
terms involving the divergence djH^^ of the bilinear potential or derivatives of the 
divergence. All these terms can be interpreted as effective vector or axial potentials 
or their derivatives. Furthermore, the resulting contributions to the EL equations are 
of degree three on the light cone. Such contributions will again be discussed in ^9.41 
below. 

9.3. Further Potentials and Fields. Having considered many different perturba- 
tions of the Dirac operator, we are now in the position to draw a few general conclu- 
sions, and to discuss a few potentials and fields which are not covered by the previous 
analysis. First of all, we point out that in Lemmas 17.31 and 17.41 the vector components 
dropped out, so that we got no contribution by an electromagnetic field. This cancel- 
lation can be understood from the general structure of our action principle, namely 
from the fact that the Lagrangian (j6.29p involves the differences of the absolute values 
of the left- and right-handed eigenvalues. As a consequence, any perturbation which 
affects the eigenvalues and A^ in the same way necessarily drops out of the EL 
equations. In other words, the EL equations are only affected by perturbations of the 
fermionic projector which have odd parity. 

The fact that the electromagnetic field does not enter the EL equations does not 
necessarily imply that the electromagnetic field must vanish. But it means that no 
electromagnetic fields are generated in the system, so that the only possible electro- 
magnetic field must be radiation coming from infinity. Having isolated systems without 
incoming radiation in mind, we conclude that our system involves no electromagnetic 
field. 

The above consideration for the electromagnetic field also applies to the gravita- 
tional field, as we now explain. Introducing gravitational fields (see for example |14^ 
§1.5]) gives rise to contributions to the light-cone expansion which involve the met- 
ric, the curvature and the derivatives of curvature (see [121 AppendixB]). The main 
effect of the gravitational field can be understood as a "deformation" of the light cone 
corresponding to the fact that the light cone is now generated by the null geodesies. 
The corresponding contributions to the light-cone expansion drop out of the EL equa- 
tions if we make the action principle diffeomorphism invariant simply by replacing the 
measure d'^x in (12. Sp by det gij \ d'^x. In addition, there are terms involving the cur- 
vature of space-time, whose singularity on the light cone is of so small degree that the 
corresponding closed chain can again be treated perturbatively. Since these curvature 
terms are even under parity transformations, they drop out of the EL equations. We 
conclude that our model involves no gravitational field. 
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From the physical point of view, it might seem disappointing that our model in- 
volves no electromagnetic and gravitational fields. However, the simple explanation 
is that our system of one sector is too small to involve these fields. If one considers 
systems of several sectors, the equation analogous to (j6.29p will involve the differences 
of the eigenvalues and A^, in different sectors. Then potentials no longer drop 
out even if they have even parity, provided that they are not the same in all sectors. 
Only a detailed analysis of systems involving several sectors will show whether the 
electromagnetic and gravitational fields will appear in the physically correct way. 

Having understood why gravitational fields drop out of the EL equations, one might 
want to consider instead an axial gravitational field as described by a perturbation of 
the form 

25 = ij^ji h}^ dj + (lower order terms) . (9-4) 

Such a field cannot occur for the following reason. As just mentioned, a gravitational 
field describes a deformation of the light cone, which due to the diffeomorphism in- 
variance of our action principle does not enter the EL equations. Similarly, an axial 
gravitational field (j9.4p describes a deformation of the light cone, but now differently 
for the left- and right-handed components of the Dirac sea. Thus the light cone "splits 
up" into two light cones, leading to highly singular contributions to the EL equations. 
Thus in order to satisfy the EL equations, the axial gravitational field must vanish. 

In order to avoid the problem of the axial deformations of the light cone, one may 
want to consider a so-called axial conformal field 

!B = iE{x)j^j^dj + (lower order terms) . (9.5) 

But this field is of no use, as the following consideration shows. Let us consider for 
a given real function A the so-called axial scaling transformation U = e'^ M^) _ This 
transformation is unitary, and thus it has no effect on the EL equations. Transforming 
the Dirac operator according to 

(i^ - m) ^U{i^ - m)U'^ = iU'^^ -m + i-f^{^A) 

= cosh(2A(2;)) i^ + sinh(2A(2;)) ij^^ -m + i-f^{^A) , 

the summand sinh(2A(x)) 27^^ can be identified with the first-order term in (j9.5p . The 
summand cosh(2A(x)) i^, on the other hand, is a conformal gravitational field, and we 
already saw that gravitational fields do not enter the EL equations. Thus in total, we 
are left with a perturbation of the Dirac operator by chiral potentials, as considered 
earlier in this paper. 

Next, we briefiy consider a scalar or pseudoscalar potential (I4.2ip . The leading 
contributions to the fermionic projector involve the potentials H and whereas to 
lower degree on the light cone also derivatives of these potentials appear. Since ^ has 
even parity, its leading contribution to the EL equations vanishes. For the potential H, 
the leading contribution cancels in analogy to (19. 2p (see also Lemma IB.ip . But to 
degree four on the light cone, one gets cross terms similar to (j9.3p which also involve 
the axial potential. In order for these additional terms to vanish, we are led to setting 
the scalar and pseudoscalar potentials equal to zero. Thus there seems no point in 
considering scalar or pseudoscalar potentials. Nevertheless, scalar or pseudoscalar 
perturbations might enter the EL equations to degree three and lower, as will be 
discussed in §9.41 below. 

In the analysis of the axial potential we made one assumption which requires a brief 
explanation. Namely, when introducing the axial potential Aa in ^6.21 we assumed that 
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it couples to all generations in the same way (see (j6.2ip ). In view of the constructions 
of ^7.4l - ^7.6l where it was essential that the potentials were different for each genera- 
tion, the ansatz (j6.2ip seems rather special, and one might wonder what would happen 
if we replaced the potential Ag, in (j6.2ip by a matrix potential acting non-trivially on 
the generations. Indeed, this scenario was already discussed in |14l Remark 6.2.3], and 
thus here we briefly repeat the main idea. Suppose that the potential Aa in (j6.2ip 
were a matrix. Then the exponentials in ()6.23p would have to be replaced by ordered 
exponentials of the form 



This is a unitary matrix whose eigenvalues can be regarded as different phase factors. 
Thus when taking the partial trace, we do not get a single phase factor, but instead 
a linear combination of different phases. As a consequence, the relations (j6.28p will in 
general be violated, so that we get a contribution to the EL equations to degree five on 
the light cone. Going this argument backwards, we can say that the EL equations to 
degree five imply that the eigenvalues (j6.27p should involve only one phase, meaning 
that the axial potential must in fact be of the form (j6.2ip with a vector field Ag,. 

9.4. The Non-Dynamical Character of the EL Equations to Low^er Degree. 

In our analysis of the EL equations we began with the leading degree five on the 
light cone. The analysis to degree four revealed the field equations and thus described 
the dynamics. Generally speaking, to degree three and lower on the light cone, we 
get many more conditions, but on the other hand, we also get much more freedom 
to modify the fermionic projector. Namely, to degree three the Dirac currents as 
well as vector and axial potentials give rise to many terms, which in general do not 
cancel each other, even if the field equations of Chapter [8] are satisfied. For example, 
the factors T^^^ involving curly brackets come into play (see (15. 5[) ). and we also get 
contributions involving higher derivatives of the potentials. To degree two and one on 
the light cone, we get even more terms, involving the cross terms of different potentials 
and the terms generated by the mass expansion of the fermionic projector. In order 
to satisfy the EL equations to lower degree, all these terms must cancel each other. 
The good news is that we also get more and more free parameters. For example, we 
can consider scalar and pseudoscalar potentials (I4.2ip . or bilinear potentials, or the 
other potentials discussed in §9.31 All these potentials can be chosen independently 
on the generations. Taking into account that to lower degree on the light cone, more 
and more cross terms between different perturbations come into play, we obtain a very 
complicated structure involving a large number of free parameters to modify the EL 
equations to degree three and lower. 

In view of this complexity, it is not clear whether the EL equations can be satisfied 
to every degree on the light cone or not. The analysis becomes so complicated that 
it seems impossible to answer this question even with more computational effort. A 
possible philosophy to deal with this situation would be to take a pragmatic point 
of view that one should simply satisfy the EL equations as far as possible, but stop 
once the equations can no longer be handled. Since we do not find this point of view 
convincing, we now go one step further and explain why the analysis of the EL to 
lower degree (no matter what the results of this analysis would or will be) will have 
no influence on the dynamics of the system. 
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The field equations of Chapter [8] were dynamical in the sense that they involved 
partial differential equations of the potentials, and by solving these equations one finds 
that the potential is non-trivial even away from the sources. This property of the field 
equations is a consequence of the fact that the leading contributions to the fermionic 
projector (namely the phase factors in (j6.23p ) dropped out of the EL equations. This 
made is possible that the derivative terms became relevant (although they were of 
lower degree on the light cone), leading to dynamical field equations. The set of 
perturbations which can lead to dynamical field equations is very limited, because 
this requires that the potential itself must drop out of the EL equations, meaning 
that the potential must correspond to a local symmetry of the system. All in this 
sense dynamical perturbations have been considered in this paper. This implies that 
all further potentials and fields will be non- dynamical in the sense that the potentials 
themselves (and not their derivatives) will enter the EL equations. This would give rise 
to algebraic relations between these potentials. In particular, these potentials would 
vanish away from the sources. Thus they do not describe a dynamical interaction, also 
making it difficult to observe them experimentally. 

A possible idea for avoiding non-dynamical potentials is to choose the potentials 
differently for each generation, in such a way that the most singular contribution 
vanishes when taking the partial trace over the generation index. This would open the 
possibility that the leading contribution to the EL equations might involve derivatives, 
thereby giving rise to dynamical field equations. We now give an argument which 
explains why this idea does not seem to work: Suppose that the potential is chosen as 
a matrix on the generations. Then for the potential to drop out of the EL equations, we 
must impose that a certain partial trace involving the potential must vanish at every 
space-time point. But this implies that this partial trace also vanishes if the potential 
is replaced by its derivatives. In particular, these derivative terms again drop out of 
the EL equations, making it impossible to get dynamical field equations. 

This concludes our analysis of local potentials. Our treatment was exhaustive in 
the sense that we considered all multiplication operators and considered all relevant 
first order operators. Our results gave a good qualitative picture of how the fermionic 
projector is affected by different kinds of perturbations of the Dirac operator, and what 
the resulting contributions to the EL equations are. Clearly, the present paper cannot 
cover all possible perturbations, and some details still remain unsettled. Nevertheless, 
as explained above, our analysis covers all perturbations which should be of relevance 
for the dynamics of our fermion systems. 



So far, we only analyzed the EL equations at the origin, i.e. for the leading contri- 
bution in an expansion in the parameter ^ = y — x. But of course, the EL equations 
should also be satisfied away from the origin. In this chapter, we will explore whether 
this can be accomplished by introducing nonlocal potentials into the Dirac equation. 
We shall see that this method will indeed make it possible to satisfy the EL equations 
to every order in an expansion around = (see Theorem 110. 5p . But it will not 
become possible to satisfy the EL equations globally for arbitrary x and y (see ^10. 4p . 

In order to introduce the problem, we consider the perturbation of the fermionic 
projector by a particle wave function i.e. in view of (j4.10p 
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In the variable y — x, this contribution oscillates on the scale of the Compton wave 
length, whereas in the variable y + x, it varies typically on the larger atomic or macro- 
scopic scale. For this reason, it is appropriate to begin by analyzing homogeneous 
perturbations which depend only on the variable y — x. In ^10.31 we shall extend our 
constructions to build in an additional dependence on the variable y + x. 

10.1. Homogeneous Perturbations of the Fermionic Projector. In order to 
explain the basic idea in the homogeneous setting, we again consider a system of one 
vacuum Dirac sea (I5.25|) for a given rest mass m > 0. This system is composed of 
states ^k{x) = x{k)^~^^^ which are plane waves of momentum k on the lower mass 
shell. A perturbation of the system which has not been considered so far is to vary the 
momenta of the states, dropping the mass shell condition. Considering a variation 5k 
of the momentum leads us to replace [x) by the plane wave 

e-*(^+'^'=)^ . (10.2) 
Then to first order in the individual states are varied by 

5^>k{x) = -i (x, 5k) ^>k{x) + 0{{5kf) , 
and this yields the following perturbation of the fermionic projector, 

5P{x,y) = i j -^{t5k){^ + m)5{k^ -m^)e{-k^)e-''^^-yK (10.3) 

The vector field 5k can be chosen arbitrarily on the lower mass shell, giving us a lot 
of freedom to vary P{x,y). However, due to the factor ^ inside the inner product, the 
variation necessarily vanishes at the origin x = y {at least under the reasonable as- 
sumption that the function 5k has suitable decay properties at infinity). Our idea is to 
use perturbations of the form (110. 3p to compensate for the effect of the fermionic wave 
functions (110. ip in the EL equations away from the origin. At the origin, where (110. 3p 
has no effect, our previous analysis remains valid, and the EL equations still lead to the 
field equations as worked out in Chapter [8l But the perturbations (|10.3p will justify 
why it is unnecessary to consider the higher order terms in an expansion around ^ = 0. 

The aim of this chapter is to work out this idea quantitatively, and to explore in 
increasing generality what the potential and the limitations of the resulting methods 
are. First of all, we remind the reader that only the axial component of the currents 
enters the EL equations to degree four (whereas the scalar, pseudoscalar and bilinear 
components drop out; see ^9. II and ^9.2p . But the perturbation (jlO.Sp does not have an 
axial component, and thus we must generalize (|10.3p such as to include a perturbation 
which is odd under parity transformations. To this end, we choose a vector field q in 
momentum space with the properties 

{k,q{k))=0 and q{kf = -1 . (10.4) 

Then the operators 

n±(A;) :=^{l^j'4ik)) 

are projectors which commute with the fermionic projector of the vacuum; they project 
onto the two spin orientations in direction of q and —q, respectively. Thus multiplying 
the fermionic projector by II-i-, 

P±{k) = U±{k) + m) (5(A;2 - m^) Q{-k^) , 
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we decompose the Dirac sea into two "subseas" P±, which are still composed of solu- 
tions of the free Dirac equation. We remark that this decomposition was already used 
in |14[ §C.l], where it was shown that in a suitable limit ni\,0, the projectors P±{k) 
go over to chiral Dirac seas composed of left- or right-handed particles. Here the above 
decomposition gives us the freedom to vary the momenta of each subsea independently. 
Thus we generalize (jlO.Sp by 

6P{x,y) = ^ E / ^A^^^''s)Ps{k) e-'=(^-^) , (10.5) 

where 6k^ and 5k- are two vector fields on the lower mass shell. For our purposes, 
it will be sufficient to always assume that these vector fields are smooth and have 
rapid decay. Thus we can extend them to Schwartz functions in momentum space, 
6k± G S{M). As a consequence, the function 6P{x,y) is smooth (but due to the 
restriction to the mass shell in (jl0.5p . it will in general not have rapid decay). Clearly, 
more general functions 5k± could be realized by approximation. 

In the EL equations, 6P{x,y) is contracted with a factor ^. The corresponding 
vector and axial components are computed by 

^Tr(OP(x,y)) =^ I ^{^,5k++6k^){^,k)T^.{k)e^^^ (10.6) 

^Tr(7^^5P(x,y)) J ^ 6k+ - 6k.) {^,mq) T^^ik) e''^ , (10.7) 

where the distribution T^2{k) := 5{k'^ — m?) Q{—k^) is supported on the lower mass 
shell. Collecting the factors of ^, these expressions can be written in the form ^j^iA^^ 
with a symmetric tensor field A>\£^). We first want to eliminate the tensor indices, 
leaving us with scalar Fourier transforms. Thus suppose that for a given smooth 
tensor field A^\$^) we want to find the corresponding vector fields 5k± and q in ()10.6p 
or ()10.7p . In the vector component ()10.6p . we can rewrite the factor {£,,k) as a ^- 
derivative, leading to the equation 

where // is the Fourier integral 

Integrating the ordinary differential equation (jlO.Sp gives the solution 

Jo 

being a smooth vector field. Thus for every choice of the vector index / we must solve 
the equation 

^f(k)T^.{k)e^'^, (10.9) 

where we set f = fi and / = {6k^ + 5k^)i/2. In this way, we have reduced (jl0.6p to 
scalar Fourier integrals of the form (jlO.Op . The same can be accomplished for the axial 
component (jlO.Tp with the following construction. We write the factor ^ in (jlO.Tp . 
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which is contracted with q, as a /c-derivative of the factor e*^^ and integrate by parts. 
Using the relation 

where in the last step we applied the orthogonality relation in (|10.4p . we obtain the 
equation 

1 f d^k 



2tt) 



d 

-^[{5k+- 5k^)imq^ 



T^2{k) e 



iki 



We again fix the index /, but now set /(^) = ^-^ Aji{^). Furthermore, we introduce 
the vector field = {Sk+ — Sk-.)i mq^ . Suppose that the smooth scalar function / 
can be represented in the form (jlO.Op with a suitable Schwartz function /. Then the 
remaining task is to satisfy on the lower mass shell the equation 

-^^'(k) = f- (10-10) 

In order to verify that this equation always has a solution, it is useful to rewrite it as a 
geometric PDE defined intrinsically on the hyperbola Ti = {k £ M \ k'^ = rm?, k^ < 0}. 
Namely, the orthogonality condition in ()10.4p implies that the vector v is tangential 
to Ti, and the derivatives in (jlO.lOp can be rewritten as the covariant divergence on Ti, 

Vjv^ = f e Sin) (10.11) 

(where V is the Levi-Civita connection on T-L). Conversely, for a given vector field v 
on Ti satisfying (jlO.lip , extending t> to a vector field in M which is everywhere orthog- 
onal to k gives the desired solution of (jlO.lOp . A simple solution of (jlO.lip is obtained 
by first solving the Poisson equation A-^(/) = / (for example using the explicit form of 
the Green's function on the hyperbola) and setting v = V<j). However, this solution 
has the disadvantage that v has no rapid decay at infinity. In order to do better, we 
must exploit that our function /(^) = Aji{^) vanishes at the origin, and thus the 
integral of / vanishes, 

^' fdpn = 0. (10.12) 



H 

Combining this fact with the freedom to add to v an arbitrary divergence-free vector 
field, one can indeed construct a solution v of (jlO.lip within the Schwartz class, as is 
shown in the following lemmsd. 

Lemma 10.1. For every function f £ S{T-L) whose integral vanishes (110. 12p . there is 
a vector field v G S{7i) which satisfies (jlO.lip . 

Proof. We parametrize the hyperbola Ti by 

^— \An^~+/j^, p cos i), p sin -d cos ip, p sin i9 sin 93^ G M, 

where p := \k\, and {'d,(p) are the standard polar coordinates on the 2-sphere. Then 
the metric on Ti is diagonal, 

gij = diag ( — o^P^^P^ sin^ ^ 
\p + m'^ 



^\ thank Bernd Ammann for the idea of solving the equation on the leaves of a foliation, after 
subtracting the mean value of /. 
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For the vector field v we take the ansatz as the sum of a radial part = (u^, 0, 0) and 
an angular part = {0,w'^ ,w'^). We also regard w as a vector field on the sphere. 
Then the equation (jlO.lip can be written as 

div-H(n) + divg^iw) = f . (10.13) 

Taking the average of / over spheres defines a spherically symmetric Schwartz func- 
tion 7, 

lip) [ fiP^^^v) d^dcos^ e Sin) . 

47r Js2 

Writing the integration measure as d/x-^ = ^/detg dp d-d dip and using Fubini, the con- 
dition (|10.12p implies that 

7(P) /'^^ =0- (10.14) 

Solving the Poisson equation l^-u(j) = f, the resulting function (j) is smooth and again 
spherically symmetric. Setting u = 'V<j3, we obtain a smooth radial vector field, being 
a solution of the equation 

div^(n)=7. (10.15) 

Writing the covariant divergence as 

d-^W = -7^^ (Vd^gnA = ( -^££=] , (10.16) 

the divergence condition (jlO.lSp becomes an ordinary differential equation, having the 
explicit solution 

^ (P) = J2 / fiV 



P^ Jo \/r2 + m2 ' 

Using (jl0.14p . one immediately verifies that the vector field u and all its derivatives 
have rapid decay at infinity. Hence u is in the desired Schwartz class. 
Using (|10.15p in (|10.13p . it remains to consider the differential equation 

divs<w) = f-7. (10.17) 

Since the function f — f has mean zero on every sphere, it can be expanded in terms 
of spherical harmonics starting at / = 1, 

oo I 

1=1 k=-i 

Since / — / is a Schwartz function, the coefficients cik are all smooth in p, and these 
coefficients together with all their p-derivatives have rapid decay in both p and I, 
uniformly in k. Hence the Poisson equation As2(l) = f — f can be solved explicitly by 

i=ik=-i H' + iJ 

defining again a Schwartz function on Ti. Introducing the vector field w hy = 
^^2'/' = [0, d^cj), cos~'^ dip(j)) gives the desired solution of (|10.17p in the Schwartz 
class. □ 
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10.2. The Analysis of Homogeneous Perturbations on the Light Cone. Fol- 
lowing the above arguments, it remains to consider the scalar Fourier transform (jl0.9p . 
Having a weak evaluation on the light cone (j5.7p in mind, we may restrict attention 
to the light cone L = | = 0}. Thus our task is to analyze the Fourier integral 

fiO = / />) '^(^' - "^') ®(-^°) e*'^ for e G ^ • (10.18) 

More precisely, for a given smooth function / G C°°{M) we want to find a Schwartz 
function / G 5(iVf) such that the Fourier integral (jlO.lSp coincides on the light cone 
with /. The question is for which / such a function / exists. We begin the analysis 
in the simple case that the function / when restricted to the mass shell depends only 
on the variable oj = (by linearity, we can later realize more general functions / by 
superposition). Then the resulting Fourier integral is spherically symmetric, so that 
the function /(^) will only depend on the time and radial variables t = and r = |^|. 
Restricting attention to the light cone t = itr, we end up with a one-dimensional 
problem. More precisely, setting p = and denoting the angle between ^ and k by ■!?, 
the Fourier integral (jlO.lSp becomes 

f{t,r) = ^ dojf{u;)e''^' dp diu"^ - - rr?) dcos i? 6"^^"™^'' 

8vr^ J~oo Jo J-i 

/-oo 



87r3; 



oo 
— m 



doj f{uj) e"^' / p dp 6{u}^ - p' - 
Jo 



167r3 

Hence on the light cone t = itr we obtain the representation 



167r3 

where we set 



^*/(*) = TTTT / dco f{u) (e-+* - e^^*) , (10.19) 



uj± := uj ± \pJ^ — fn^ 

The right side of (jlO.lOp differs from an ordinary Fourier integral in two ways: First, 
the integrand does not involve one plane wave, but the difference of the two plane 
waves e*'^^*, whose frequencies are related to each other by w+W- = vr? . Second, 
in (110. 19p only negative frequencies appear. Let us discuss these two differences af- 
ter each other. The appearance of the combination (6**^+* — e*'^~*) in (jlO.lOp means 
that the coefficients of the plane waves cannot be chosen arbitrarily, but a contri- 
bution for a frequency a;_ < —m always comes with a corresponding contribution 
of frequency a;+ > — m. This frequency constraint makes it impossible to represent 
a general negative-frequency function /; for example, it is impossible to represent a 
function tf{t) whose frequencies are supported in the interval [— oo,— m). However, 
this frequency constraint can be regarded as a shortcoming of working with a single 
Dirac sea. If we considered instead a realistic system of several Dirac seas ()3.ip . the 
Fourier integral ()10.19p would involve a sum over the generations, 

it fit) = -^Y.J du (e-+* - e--*) (10.20) 
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with uj^{uj) := a;ib(a;^— m|)~ 2 and g > 1. Then the freedom in choosing g independent 

functions would indeed make it possible to approximate any negative-frequency 
function, as the following lemma shows. 

Lemma 10.2. Assume that the number of generations g > 2. Assume furthermore that 
the ordinary Fourier transform f of a given Schwartz function f £ S{M) is supported 
in the interval (— oo,0). Then there is a sequence of Schwartz functions fn G 5(M) 
such that the corresponding functions fn{t) defined by the Fourier integrals ()10.20p as 
well as all their derivatives converge uniformly to f{t), 

sup|af(/n(t)-/(t))| forallK>0. (10.21) 

Proof Let us try to find functions f^ in (jl0.20p such that the right side of (jl0.20p 
gives the plane wave e*^* with < 0. Again ordering the masses according to ()3.2p . 
we choose wi such that uj^{uji) = ri or uj^{uji) = Q, i.e. 

Q"^ + ml 

Then choosing P{u}) = ±^(0; — wi), we obtain the desired plane wave e*^*, but as an 
error term we get the plane wave — e*^i* with Q.i = rn^/Q. In order to compensate the 

error, we next choose 002 such that u}\{oj2) = or ^"^{002) = ^^i. Choosing f^{uj) = 
6{u! — UJ2), the plane wave — e*^i* drops out, but we obtain instead the plane wave e*^^* 
with 0,2 = rnf/Qi = m\Q./m'^. We proceed by compensating the plane waves in turns 
by the last Dirac sea and the first Dirac sea. After n iteration steps, the functions 
and take the form 

where 

_ ml'' m2"+2 

- ^ ^2n+, - ^ . 

The Fourier integral (jl0.20p gives rise to the plane waves 

^l 
,2 



^int _ ^ix-nt ^i^g^g A := ^ < 1 . 



In order to form superpositions of these plane waves, we next multiply by a Schwartz 
function h{fi) and integrate over Q. Then the Fourier integral (jl0.20p becomes 

itf^{t) = f h{n) (e*™ - e*^"^*) dn . 

Choosing h{uj) = —di^f{uj)/{2Tr), we can extend the integration to the whole real line. 
After integrating by parts, we can carry out the Fourier integral to obtain 

fn{t) = ^ fi^) {e'""' - A"e*^(^"*)) dn = fit) - X^fiX^t) . 

From this explicit formula it is obvious that the functions /„ converge in the limit n — 
00 as desired (|10.2ip . □ 
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As is immediately verified, the functions /„ as well as all their derivatives also converge 
in L^(M). However, we point out that for the functions and the convergence is 
a bit more subtle. Namely, from (jl0.22p one sees that for large n, these functions 
involve more and more contributions for large oj. A direct calculation shows that in 
the limit n — oo, these functions converge to smooth functions which decays at infinity 
only ~ l/w. 

Let us now discuss the consequences of the fact that (|1U.19|) only involves negative 
frequencies. This restriction is already obvious in the Fourier integral (110. 5p . before 
the reduction to scalar Fourier integrals (jlO.lSp or ()10.19p . Since contractions with 
factors ^ merely correspond to differentiations in momentum space which preserve the 
sign of the frequencies, the following considerations apply in the same way before or 
after the contractions with ^ have been performed. We point out that the contribution 
by a Dirac wave function (jlO.ip can be composed of positive frequencies (=particles) or 
negative frequencies (=anti-particles), and thus in (jlO.ip we cannot restrict attention 
to negative frequencies. This raises the question whether a contribution to (110. ip of 
positive frequency can be compensated by a contribution to (110. ISp of negative fre- 
quency. The answer to this question is not quite obvious, because the EL equations 
involve both P{x,y) and its adjoint P{y,x) = P{x,y)*. Since taking the adjoint re- 
verses the sign of the frequencies, a negative-frequency contribution to P{x, y) affects 
the EL equations by contributions of both positive and negative energy. Thus one 
might hope that perturbations of P{x, y) of positive and negative energy could com- 
pensate each other in the EL equations. However, such a compensation is impossible, 
as the following lemma shows. 

Lemma 10.3. Assume that f and g are the Fourier transforms of chiral perturbations 
of the fermionic projector, such that f has a non-vanishing contribution inside the 
upper mass cone, whereas g is supported inside the lower mass cone. Then the linear 
contributions of f and g to the EL equations to degree four cannot compensate each 
other for all 

Proof. By linearity, we may restrict attention to the spherically symmetric situation, 
so that / and g restricted to the light cone are functions of one variable t. Since the 
negative-frequency component of / can clearly be compensated by g, we can assume 
that / and g are composed purely of positive and negative frequencies, respectively. 
The perturbation g affects the EL equations to degree four by (see Lemma (|7.4p and 
its proof in Appendix [B]) 

n X c (Mi 5(0 + Ms + (deg < 4), 

where 



Ml = 4]'%'^ and M, = -T^.-'^'^ =^ , (10.23) 

and c > is an irrelevant constant. Evaluating weakly on the light cone (|5.7p . we 
obtain the contribution 

cig{t) + C2^, (10.24) 

where ci and C2 are real regularization parameters (real because the degree is even). 
Similarly, the perturbation / yields the contribution 



cif{t) + C2fit). (10.25) 
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In order for (jl0.25p to compensate (|10.24p . both the negative and positive frequencies 
must cancel each other, leading to the conditions 



ci g{t) = C2 f{t) and C2 g{t) = ci f{t) 



(10.26) 



Taking the complex conjugate of the first equation, multiplying it by C2 and subtracting 
ci times the second equation, we get 



We thus obtain the condition 



{4 - 4) fit) = 



Ci = ±C2 



(10.27) 



If this condition holds, we can indeed satisfy (jl0.26p by setting g{t) = ib/(t). 

We conclude that / and g can compensate each other if and only if we impose 
the relation (jl0.27p between the regularization parameters corresponding to the basic 
fractions Mi and M2. Imposing relations between the regularization parameters was 
not used previously in this paper, and one could simply reject ()10.27p by saying that 
we do not want to restrict the class of admissible regularization by introducing such 
relations. However, this argumentation would not be fully convincing, as it would not 
allow for the possibility that the microscopic structure of space-time on the Planck 
scale corresponds to a regularization which does have the special property (|10.27p . 
This possibility is ruled out by the following argument which shows that there are in 
fact no regularizations which satisfy (jl0.27p : As only the real parts of basic fractions 
enter ()5.7p , it suffices to consider the real parts of Mi and M2 ■ Using the specific form 
of these monomials in (jl0.23p . we obtain 



Re(Mi + M2) 



.(-1) 



.(0) 



In,T<»')- 



Re(Mi - M2) 



(-1) 



T, 



(0) 



ReTr 



(0) 



Both these expressions are non-negative, and thus there cannot be cancellations be- 
tween positive and negative contributions, no matter how we regularize. Without a 
regularization, we know from ()4.13p - (j4.15p that 

1 PP 



RerW 



and 



ImT, 



(0) 



Regularizing these terms, we find that for any regularization, both Mi + M2 and Mi — 
M2 are non-zero to degree four on the light cone. Hence (|10.27p is violated. □ 



We come to the definitive conclusion that using perturbations of the form (jlO.Sp . it 
is in general impossible to satisfy the EL equations to degree four globally for all ^. 
But, as we will now show, it is possible to satisfy the EL equations locally near ^ = 0, 
in the sense that we can compensate all contributions in a Taylor expansion in ^ to 
an arbitrarily high order. We consider the obvious generalization of (|10.5p to several 
generations 

5P{x,y) = iy:y: \ ^(tSks) P±{k) e-^^^--y^ . (10.28) 



a 

■'EE 

13=1 s=± 



(2vr)^ 



Proposition 10.4. Suppose that the number of generations g > 2. Then for any 
given smooth functions /ivj^a £ C°°(M) and every parameter L > 2, there are vector 
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fields 5k± and q in the Schwartz class such that for all multi-indices k with 2 < |k| < L, 
Tr U SP{x, y)) - h,{0] _ =0 = 0^ [Xr (7^^ 6P{x, y)) - h,{0] _ • (10.29) 



Before coming to the proof, we point out that, since the right side of p0.6p and (|1U.7|) 
involves two factors of it is in general impossible to satisfy p0.29|) in case \k\ < 2. If 
the functions hv or are formed by taking the vector or axial component of a contri- 
bution to the fermionic projector (for example the contribution by a particle (|lU.ip ). 
then (I10.29P is trivially satisfied in the case |k| =0, simply because all functions vanish 
at the origin ^ = 0. However, in this case the linear terms in ^ of and /la will in 
general be non-zero. Indeed, these linear terms are responsible for the field equations 
as worked out in Chapter [HI Thus (|10.29p means that all higher order terms in an 
expansion in ^ can be compensated by perturbations of the form (I10.28p . 

Proof of Proposition \10.4\ Following the arguments after (jlO.Sp and applying Lem- 
ma [inill it again suffices to analyze scalar Fourier integrals. Furthermore, using the 
polarization formula for the multi-index k, we may restrict attention to the spherically 
symmetric situation (jl0.20p with Schwartz functions G 5(M). Keeping track of 
the factors ^ in the arguments after (jlO.Sp . it remains to show that for every smooth 
function h{t) there are functions G S{M) such that the corresponding function f{t) 
defined by (|10.20p satisfies the conditions 

^(/(0-Mi))L=o = for alU = 0,...,L. (10.30) 

We choose a test function fj G Cq°((— 1, 1)) and denote its Fourier transform by r/. 
Furthermore, we choose a parameter Qq < —Arug and introduce the function r^f^g G 
C^{{-hmg,-?,mg)) by 

Thus ryf^Q is supported for large negative frequencies. For any such frequency VL G 
suppr^Hp, we want to construct the plane wave e*^*, with an error term which is 
again of large negative frequency. To this end, we proceed similar as in the proof of 
Lemma llO. 21 bv iteratively perturbing the first and last Dirac seas, but now beginning 
with the first sea. Thus we first construct the plane wave e*^* by perturbing the first 
sea, and compensate the error term by perturbing the last Dirac sea. This gives in 
analogy to (110.221) 

giving rise to the plane wave 

^iQt _ ^int ml/ml _ (10.31) 

Multiplying by i)^^ and integrating over fi, we find that the function 

1 / \ rn^ 

/(t) := _(e*^o*r?(40oi) -e*^^*r/(4S7it)) with = ^ 

it\ J mf 

has the desired Fourier representation (|10.20p . Since differentiating (|10.20p with 
respect to t merely generates factors of a;£, the functions /^'^(t) := t~^dl{tf{t)), 
I = 1, . . . ,L, can again be represented in the form (jl0.20p . By a suitable choice of 
the function fj, we can clearly arrange that the parameters /(O), f^^Ho), f^^\0) 
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are linearly independent. Thus by adding to / a suitable linear combination of the 
functions we can arrange (|lU.30p . □ 

For clarity, we point out that the function / in (jl0.30p will in general not be a good 
global approximation to h. In particular, it is impossible to pass to the limit L — )■ oo. 
We also remark that this proposition also holds in the case g = 1 of only Dirac 
sea. However, in this case it would not be possible to compensate the error term, so 
that instead of (|10.3ip we would have to work with the combination e*^* — e**"* 
This has the disadvantage that in the limit — t- — oo, we would get contributions of 
low frequency, making it impossible to generalize the result to the non-homogeneous 
situation (see the proof of Theorem 110.51 below). This is why Proposition 110.^ was 
formulated only in the case g > 2. 

10.3. Nonlocal Potentials, the Quasi-Homogeneous Ansatz. We now want to 
extend the previous results to the non-homogeneous situation. Since we only consider 
perturbations which are diagonal on the generations, all constructions immediately 
carry over to several generations by taking sums. For notational simplicity, we will 
write all formulas only for one generation. Our method is to first describe our previous 
perturbations of the fermionic projector (jl0.5p by homogeneous perturbations of the 
Dirac operator. Replacing this perturbation operator by a nonlocal operator will 
then make it possible to describe the desired non-homogeneous perturbations of the 
fermionic projector. We begin by noting that the plane wave (|10.2p is a solution of the 
Dirac equation (i^ — 6Ijt — m)^' = 0{{Sk)'^). Thus the linear perturbation (|10.3p can 
be described equivalently by working with the perturbed Dirac operator in momentum 
space 

]ji - 61ji{k) - m . 

Similarly, for the perturbation (jl0.5p . we must find a perturbation n of the Dirac 
operator which is symmetric and, when restricted to the image of the operators P±{k), 
reduces to the operators —5^^- In order to determine n, it is convenient to decompose 
the vector fields 5k± as 

5k± = V ±{(t)q + A) where {A,q) = Q (10.32) 

(thus V is the vector part, whereas (j) and A describe the components of the axial 
part parallel and orthogonal to q, respectively). Then using (|10.4p together with the 
relations 

we obtain 

-5J^±p± = (-f T (0^ + 4)) P± = -fP± + i<P4 + A){iH)P± 



m 



where in the last step we used that (^^ — m)P± = 0. The square bracket has the desired 
properties of n. Thus the perturbation (jlO.Sp is equivalently described by the Dirac 
operator in momentum space 

^ + n + m with n{k) = -f{k) + ^!^-f^I^ + -f^i{k)4{k) , (10.33) 

m 

and the fields V, A and cp as defined by (110. 32p . Note that n is composed of three terms, 
which can be regarded as a vector, an axial and a bilinear perturbation of the Dirac 
operator. The perturbations in (jl0.32p are all homogeneous. The reader interested 
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in the perturbation expansion for the corresponding Dirac solutions is referred to [14:\ 
§C.l]. 

In order to generahze to non-homogeneous perturbations, we write the operator n 
as the convolution operator in position space 

(n*)(x) = / n(x,y)V(y) A (10.34) 
Jm 

with the integral kernel 

Now we can replace n(x, y) by a general nonlocal kernel. We refer to the operator n 
as a nonlocal potential in the Dirac equation. For technical convenience, it seems 
appropriate to make suitable decay assumptions at infinity, for example by demanding 
that 

n(x,y) e 5(M X M) . (10.35) 

Then the corresponding fermionic projector can be introduced perturbatively exactly 
as outlined in §4.41 Before we can make use of the general ansatz (110. 34p and (110. 35p , 
we must specify n(x,y). To this end, we fix the variable C := y + x and consider 
the fermionic projector as a function of the variable ^ only. Then we are again in 
the homogeneous setting, and we can choose the operator n as in (jl0.33p . In order 
to clarify the dependence on the parameter C, we denote this operator by n(A;, C). As 
explained at the beginning of this chapter, we consider y + x as a macroscopic variable, 
whereas k is the momentum of the quantum mechanical oscillations. This motivates us 
to introduce the kernel n(x, y) similar to the transformation from the Wigner function 
to the statistical operator by 

y) = I <k,y + x) e^'^^y--^ . (10.36) 

Although this so-called quasi-homogeneous ansatz is quite simple, it makes it possible 
to satisfy the EL equations to degree four locally around every space-time point, as is 
made precise in the following theorem. 

Theorem 10.5. Suppose that the number of generations g > 2, and that we are given 
an integer L > 2 and a parameter 5 > 0. Then for any given smooth functions h^, /la G 
S{M X M), there is a nonlocal potential n of the form (110. 34p with a kernel n{x,y) G 
S{M X M) of the form (jl0.36p such that for all multi-indices k with 2 < \k.\ < L, 



<5. 



dl[TT{^AP{x,y))-h,{x,y)] + [Tt {-f^^ AP{x,y)) - h^{x,y)] 

x=y x=y 

Here AP denotes the perturbation of the fermionic projector to first order in the non- 
local potential n. 

Proof. For a fixed choice of the parameter and for any given we choose the homo- 
geneous perturbation 5P as in the proof of Proposition 110.41 and rewrite it according 
to (110. 33p as a homogeneous perturbation n{k, () of the Dirac equation. Introducing 
the nonlocal potential by (110. 34|) and (110. 36p . the rapid decay of the functions /iv,^a 
implies that the kernel n(x, y) has rapid decay also in ( (the rapid decay in ^ is ob- 
vious because n(.,C) £ 5(M)). Since the same is true for all derivatives, we conclude 
that n{x,y) G S{M x M). 
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The corresponding perturbation of the fermionic projector AP can be analyzed with 
the methods introduced in [12j. We first pull out the Dirac matrices to obtain 

AP{x, y) = (i^^ + m) (-^^ + [n] (x, y) 7'^ , (10.37) 

where AT„2 is the perturbation of the corresponding solution of the inhomogeneous 
Klein-Gordon equation (see [121 equations (2.4) and (2.5)]) 



AT^2[n](x,y) 



M 



d z\ \ d Z2 



M 



X {Sm^{x,zi)n{zi,Z2)T^2{z2,y) + T^2{x,zi)n{zi,Z2) S^2{z2,y)^ 

(where Ta and Sa are again given by (I4.12|) and (I8.21|) ). We next transform to mo- 
mentum space. Setting 



n{p,q) = / n(p,C)e'2^ d^C 

JM 



with the "macroscopic" momentum vector q, the above formula for AT^2 becomes 
(see [m equations (3.8) and (3.9)]) 



AT^2[n](p+|,p-| 



(P + o ) "(P> (l)Sm?[p 



1 
2 

Now we can perform the light-cone expansion exactly as in |121 Section 3]. This gives 
(cf. [m equation (3.21)]) 

q- 



AT^2[n](p+|,p 2 



n=0 



{-IT 



2\" °o 

^ k 



-n(p, q) 
1 



A: 

^ (2A; + 1)! E 
^ ' M=o 



2k 
I 



2 Sp? 



where the curly brackets are combinatorial factors whose detailed form is not needed 
here (see [12l equation (3.13)]). 

Let us discuss how the Fourier transform of this expansion behaves in the limit 
—00. According to the construction of the functions in the proof of Prop osit ion 1 1 . 4l 
the function n{p,q) is supported in the region p'^ ~ Moreover, the p-derivatives 
of n scale in powers of l/Jlo- This implies that every derivative of the factor T^2 gives 
a scaling factor of J^q^. Since every such derivative comes with factor q"^, we obtain 
a scaling factor {q"^ / QqT'^K Thus in the limit Q,q — )• —00, it suffices to consider the 
lowest summand in n + I, 



Ar„2[n] = -nip,q)Y^ 



k=0 



1 



{2k + l)\ V 2 dp^ 



When transforming to position space, the p-derivatives can be integrated by parts. If 
they act on the function n{p,q), this generates scale factors of the order 0(|g|/|J7o|) 
which again tend to zero as f^o — >• —00. Thus it remains to consider the case when 
these derivatives act on the plane wave e*^^. We thus obtain 
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We now consider a Taylor expansion in ^ around ^ = up to the given order L. This 
amounts to replacing the factor e*^^^^ by its power series and collecting the powers 
of ^. The remaining task is to compare the factors p^, with q^,. This is a subtle 
point, because the fact that p^ ~ Qq does not imply that the inner product is 
large. Indeed, this effect was responsible for the appearance of low frequencies in the 
Fourier integral (|lU.19p . However, in the proof of Proposition I1U.41 we arranged by a 
suitable choice of /i and fg that these low-frequency contributions cancel, so that we 
were working only with the high-frequency terms (|1U.31|) . Restating this fact in the 
present context, we can say that for the leading contribution to AT„2, the factor 
is larger than qS^ by a factor of the order 0(|g|/|r2o|). Thus it suffices the consider the 
summand A; = 0, 

AT„,M(x..)^-/ /|i,„(,,)TS](ri.-¥,o(i-), 

Now we can carry out the q-integration to obtain 

AT^.[n]{x,y) = -16 J ^ n(p, y + x) T« (p) e^P^ + o[^). 

Using this result in ()10.37p . we can carry out the derivatives to recover precisely the 
homogeneous perturbation (jlO.Sp for fixed C- D 

The scaling argument used in the last proof can be understood non-technically as 
follows. It clearly suffices to consider the region where y lies in a small neighborhood 
of X. Thus we may perform the rescaling C — ?• ^^/A with a scale factor A ^> 1, leaving C 
unchanged. This corresponds to changing the momentum scale by ~^ XQq. In the 
limiting case A — )• oo, the fermionic projector depends on on a smaller and smaller 
scale. On this scale, the dependence on the variable x + y drops out, so that the quasi- 
homogeneous ansatz (jl0.36p becomes exact. For this argument to work, one must 
ensure that in the homogeneous setting all frequencies scale like ^Iq, as was arranged 
in (fTa3T]l . 

10.4. Concluding Remarks. We conclude this chapter with a brief discussion of 
our results. We first point out that, since the orders \k\ = and |k| = 1 are excluded 
in Theorem 110. 5( the results of the previous chapters are not affected (see also the 
paragraph after the statement of Proposition 110.41) . In particular, the field equations 
of Theorem 18.21 remain valid. The main conclusion from Theorem 110.51 is that the 
higher orders in ^ of an expansion of the EL equations to degree four can always be 
satisfied by suitable nonlocal potentials. These nonlocal potentials are non-dynamical 
similar as explained in §9.41 for local potentials. Furthermore, in view of the freedom 
in the constructions of Proposition 110.41 these nonlocal potentials are not uniquely 
determined, so that we do not get any physical predictions. Our point of view is 
that the nonlocal potentials give us a lot of freedom to perturb the higher orders in ^, 
explaining why we do not get dynamical field equations for tensor fields of higher rank. 

A surprising conclusion of our analysis is that, although the EL equations can be 
satisfied locally to every order in ^, it is in general impossible to satisfy them globally 
for large ^ (see Lemma 110.31 and the paragraph before this lemma). At first sight, this 
might seem to imply that the EL equations are overdetermined and cannot be solved. 
On the other hand, the general compactness results in [21] indicate that our action 
principle does have non-trivial minimizers, so that the EL equations are expected to 
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admit solutions. Thus there should be a way to compensate the above nonlocal error 
terms. A possible method is to modify the wave functions globally in space-time. 
Whether and how in detail this is supposed to work is a difficult question which we 
cannot answer here. Instead, we explain what this situation means physically: Suppose 
that a physical system is described by a minimizer of our action principle. Then 
the corresponding EL equations to degree four do not only yield the field equations, 
but they give rise to additional conditions which are nonlocal and can therefore not 
be specified by a local observer. In [15l Section 7] such so-called nonlocal quantum 
conditions were proposed to explain phenomena which in ordinary quantum mechanics 
are probabilistic. 

More specifically, the fact that the EL equations cannot in general be satisfied 
globally might explain the tendency for quantum mechanical wave functions to be 
localized, as we now outline. Suppose that the fermionic projector is perturbed by a 
fermionic wave function (jlO.ip . At the origin ^ = 0, this perturbation leads to the field 
equations as worked out in Chapter [8l The higher orders in can be compensated by 
nonlocal potentials. But the contribution for large ^ cannot be compensated, thereby 
increasing our action. Thus seeking for minimizers, our action principle should try 
to arrange that the contribution (jlO.ip vanishes for large ^. This might explain why 
quantum mechanical wave functions are usually not spread out over large distances, 
but are as much as possible localized, even behaving as point particles. This idea is 
elaborated further in 1191 . 



Appendix A. Testing on Null Lines 

In this appendix we justify the EL equations in the continuum limit (j5.29p by speci- 
fying the wave functions ^'i and ^2 used for testing ()5.24p in the setting with a general 
interaction and for systems involving several generations. Our method is to adapt the 
causal perturbation expansion ()4.9p to obtain corresponding expansions for ^'i and ^'2- 
We then consider the scaling of these terms to every order in perturbation theory. We 
rely on results from |14j and [22], also using the same notation. 

We begin with the Dirac equation for the auxiliary fermionic projector of the general 
form (14. Sp . where we assume that !B is a multiplication operator which is smooth and 
decays so fast at infinity that 

d!^'B{x) \ d^x < 00 for all multi-indices / and J with |/| < 2 . (A.l) 

Under this assumption, every Feynman diagram of the causal perturbation expan- 
sion (14. 9p is well-defined and finite (see jl4l Lemma 2.2.2]). As in [141 §2.6], we denote 
the spectral projectors of the operator {i^+'B — mY— /il) by In contrast to (I5.26p . 
where we cut out an tj-strip around the mass shell, it is here more convenient to remove 
neighboring mass shells by setting 

/Am ^ 
P+f, r] dfi , (A.2) 
-Am 

where the tilde again denotes the partial trace over the generations. When taking 
the product P^i, we get cross terms involving different generations. However, as in 
the proof of [14, Theorem 2.6.1] one sees that these cross terms vanish in the infinite 
volume limit. Thus indeed lies in the kernel of the Dirac operator. Moreover, by 
choosing Am sufficiently small, we can make the difference — r] as small as we like. 
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The construction of ^'2 is a bit more involved. In order to get into the framework 
involving several generations, we first extend the wave packet in (j5.27p to an object 
with 4(7 components, 

:= {i^ + mY)e with e = (e~'"^*+^) (/>(t + x - £, z)) • • • ■ (A.3) 

g—l summands 

For ^2 we make an ansatz involving a partial trace over the generation index, 

9 

*2 = J](^/3 + AV^D + AV'|), (A.4) 

/3=1 

where the corrections ijj^ and A'i/'^ should take into account that the auxiliary Dirac 
equation must hold and that the generalized energy must be negative, respectively. In 
order to specify A^^, we first apply the free Dirac operator to ip, 

{i^ - mY) ip = {i^- mY)(i^ + mY) (^e-*^(*+^) ^(t + ^ - y, z) • • • o) 

= -(□ + m^) (^e-'^(*+^) </.(t + X - ^, y, z)) • • • = + _ ^2^ q ^ 

where m = mi is the mass of the first generation. We point out that the obtained 
expression does not involve the "large" parameter |r2|; in this sense it is a small error 
term. In order for ^p + Atp^ to satisfy the auxiliary Dirac equation (|4.5p . the wave 
function A^^''"'^^ must be a solution of the inhomogeneous Dirac equation 

(i^ + S - mY) A^° = - (5^ + 5^ - m^) 6 - 'Bip . (A.5) 

Solutions of this equation could be constructed rigorously with energy estimates (see for 
example [33] )■ But here we are here content with a perturbative treatment. Denoting 
the Green's function of the zero mass free Dirac operator by s, i.e. 

i(^s = 1 , 

we can solve for A^^ in terms of the perturbation series 

00 

AV'^ = -Yl ("^(^ " "'^))'' ^ ii^y + dl-m^)f + ■ (A.6) 

fc=0 

We point out that A-^'^ is not uniquely determined, and this non-uniqueness is reflected 
by the fact that there is the freedom in choosing different Green's functions, like the 
advanced or retarded Green's functions or the Feynman propagator (for details on the 
above operator expansions and the different Green's functions see |14] or [22]). For our 
purpose, it is most preferable to work with the retarded Green's function s^, whose 
kernel s^(x,y) is given explicitly by (see [HI §2.5]) 

s\x,y) = -^ime)Q{-e). (A.7) 

This has the advantage that the support of /^ij)^ lies in the future of ip, and thus it is 
disjoint from the support of ^1 (see Figured] on page I25|) . 

The function ■0 + A^/;^ solves the auxiliary Dirac equation, but it will in general 
have a component of generalized positive energy. This positive-energy contribution 
must be subtracted in order to obtain a vector in the image of P. Formally, this can 
be achieved by setting 

A^/-^ = -(1 -P) (V' + A0°) . (A.8) 
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In order to give this equation a meaning, one must keep in mind that the normahza- 
tion of the fermionic projector involves a ^-distribution in the mass parameters, i.e. 
P+^P+^' = S{fi — iJ,')P-f-fj,. Thus using the formahsm introduced in [22l Section 2], 
we can make sense of (jA.SP as an operator product simply by omitting the resulting 
(5-distributions. With ()A.4p as well as ()A.3p . ()A.6P and (jA.Sp . we have introduced ^'2 
by a well-defined perturbation series. 

We now estimate for large with a similar method as previously used in [131 
Theorem 3.4] for the estimate of the non-causal high energy contribution. 

Lemma A.l. To very order n in perturbation theory and for every v G there is 
a constant C{n,u) such that the wave function Aip^ as defined by (|A.8p satisfies the 
inequality 

C{n,u) 



sup 



< 



(A.9) 



Proof The n*^ order contribution (Aip^)^"^ can be written as a finite number of terms 
of the form 

5:=C„SC„_iS ••• SCo0, (A.IO) 
where every factor Ci is a linear combination of the operators p, k, and s. Here (j) 
stands either for the wave function ip in (lA.SP or for the square bracket in (IA.6p . In 
either case, ip is given explicitly and involves the free parameter Q. It is preferable to 
proceed in momentum space. The regularity and decay assumption (lA.ip implies that 



sup 

fceM 



Setting Fq = (j) and 

m) -- 



< 00 



for all multi-indices I and J with 1/1 < 2 



(A.ll) 



d^q 



2tt) 



- q) Ci{q) Fi^iiq) (where 1 < I < n) , 



(A.12) 
(A.13) 



we can write the Fourier transform of g as 

g{k) = Cn{k) Fnik) . 

It clearly suffices to prove the lemma for v an even number. Let us show inductively 
that the functions Fi satisfy the bounds 

3 



sup (a; 

(aj,fc)eM 



< C{1, u) uniformly in Cl . (A. 14) 



i=0 



In the case / = 0, the claim follows immediately from the explicit form ip. To prove 
the induction step, we use the inequality 



{u - ny < c(i/) ((w - oj'Y + (w' - ny) 



to obtain the estimate 

\{uj -ny Fi{uj,k)\ 
, , dw' 



+ c(i/) 



3 27r 
°° duj^ 
^00 27r 



dk' 
dk' 



{2n 



{oj - Jy -u',k- k') Ci{J, k') Fi^iioj', k') 

3 B(w -J,k- k') Ci{J, k') \{oj' - nyFi_i{J, k') 
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Furthermore, the factors Ci involve at most first derivatives; more precisely, they are 
bounded in terms of Schwartz norms by (see for example [14^ Proof of Lemma 2.2.2]) 

< const ||/||4,i foran/G5. 

Combining these inequalities, we can use the induction hypothesis together with (jA.lip 
to bound the expression — Fi(u},k)\ uniformly in 0, and {uj,k). The expres- 
sion \(u — diFi{u},k)\ can be estimated in exactly the same way if one keeps in 
mind that if we differentiate ()A.12p with respect to k, the derivative acts only on the 
potential S, but not on the factor Fi^i. This proves ()A.14|) . 

We next consider the operators in (jA.lOp in more detail. The factor (1 — P) in (jA.SP 
can be regarded as a projector onto the generalized positive-energy solutions of the 
Dirac equation. The perturbation expansion of the fermionic projector can be ar- 
ranged in such a way that each operator product involves at least one factor p — k 
which projects onto the negative-energy solutions (for details see [22]). Similarly, the 
factor (1 — P) in (|A.8P implies that we can arrange the operator products such that 
every contribution (|A.10p involves at least one factor p + k, being supported on the 
upper mass cone. Thus in the corresponding induction step, we may replace an arbi- 
trary even number of factors (w — O) by factors of jfi). In the following induction steps 
we proceed as in ()A.14p . At the end, we apply ()A.13P to obtain the result. □ 

We can now prove the main result of this appendix. 

Proposition A. 2. Consider a fermion system in Minkowski space with an interac- 
tion S which is a multiplication operator satisfying the regularity and decay assump- 
tions (jA.ip . Assume furthermore that the pole of Q is of order o{\^\'"^) at the origin 
(see Definition \7.2^ . Then the EL equations (I5.20p imply that the operator Q vanishes 
identically in the continuum limit (I5.29p . 

Proof We introduce the wave functions ^'i and ^'2 perturbatively by ()A.2p and ()A.4p 
with ip, Aip^ and Aip^ according to ()A.3p . (|A.6P and (jA.Sp . Evaluating the commu- 
tator [-P, Q] as in ()5.23p gives the condition ()5.24p . Following the arguments in ^5.2\ 
the leading terms give ()5.29p . and thus it remains to consider all correction terms. 
The corrections of ^1 can be made arbitrarily small by choosing the parameter Am 
in ()A.2p sufficiently small. Working in ()A.6P with the retarded Green's function ()A.7p . 
the support of does not intersect the support of rj, so that the corresponding 

contribution to (]5.24p is well-defined in the continuum limit. 

The wave function Atp^ is more problematic, because it will in general not vanish 
on the support of rj. But according to Lemma lA.ll we can make A^^ arbitrarily 
small by choosing \Q\ sufficiently large. This is not quite good enough for two reasons: 
First, the integrand in (j5.24p becomes more and more oscillatory as \ is increased, 
so that the leading contribution to (|5.24p will also become small as \Q\ gets large. And 
secondly, even if Aip'^ is small, it gives rise to a contribution at x = y where Q{x,y) 
is ill-defined. The first problem can be treated by noting that the oscillations in 
the integrand of (I5.24p will give rise to a polynomial decay in 0, (typically a l/fi 
behavior), whereas according to ()A.9p . the wave function A^^ decays in ft even rapidly. 
Thus we can indeed arrange that (]A.9p is arbitrarily small compared to the leading 
contribution in (I5.24p . For the second problem we need to use that the pole of Q is 
of order o(]^]~^) at the origin: Due to this assumption, the integrand in (I5.24p will be 
at most logarithmically divergent at x = y. By modifying ^'2 by a suitable negative- 
energy solution of the Dirac equation (for example a wave packet of negative energy. 
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whose amplitude is fine-tuned) , one can arrange that this logarithmic divergence drops 
out. Then the integrals in (j5.24p become finite, and by choosing sufficiently large, 
we can arrange that the contribution of A-i/;^ to (j5.24p is much smaller than the leading 
contribution which yields (|5.29p . □ 

Before discussing the result of this proposition, we estimate an operator product which 
is similar to (jA.lOp but involves a nonlocal potential as considered in ^10.31 

Lemma A. 3. We consider the expression 

g = CnBn C„_i Bn-i ■ ■ ■ Bi CqcP 5 

where every factor Ci stands for one of the operators p, k, or s. As in the proof 
of Lemma \A.1[ the function cj) is either the wave function ip in ()A.3p or the square 
bracket in (|A.6p . Each factor Bi either stands for the multiplication operator B sat- 
isfying ()A.ip . or else it is a nonlocal operator n of the Schwartz class (|lU.35p . We 
assume that at least one factor Bi is a nonlocal operator. Then for every integer u 
there is a constant C{v) such that 

sup \g{x)\ < ^ . (A.15) 

Proof. As in the proof of Lemma lA.H we proceed inductively in momentum space. 
Suppose that p is the smallest index such that Bp = n. Then for all / < p, only 
the potential 23 is involved, and the functions Fi defined by (|A.12p again satisfy the 
inequalities ()A.14p . In the p^^ induction step, we must replace (IA.12P by 

where h £ S{M x M) denotes the Fourier transform of n{x,y). Using the induction 
hypothesis (IA.14P together with the rapid decay of h{p, q) in the variable q, we obtain 
a factor IJ^I"*^. In the remaining induction steps, we can use the simpler method of [141 
Lemma 2.2.2] to obtain the result. □ 

The setting of Proposition lX^ is too special for our applications, making it necessary 
to consider the following extensions: 

• Taking into account the wave functions of the particles and anti-particles: We 
first note that, being solutions of the Dirac equation, the wave functions of the 
particles and anti-particles in ()4.10p are orthogonal to the wave function ^'i 
(as is obvious from (|A.2p ). Furthermore, by choosing much larger than the 
energies of all particle and anti-particle wave functions, we can arrange that 
these wave functions are also orthogonal to ^'2. Then all the wave functions 
drop out of ()5.24p , so that we are back in the setting of Proposition IA.2I 

• Handling the local axial transformation: Following the constructions in §7.61 
we must apply the local axial transformation (j7.24p to the fermionic projector 
before taking the partial trace. Likewise, we here apply this transformation to 
the wave functions ^'i and ^'2 before taking the partial trace. This additional 
transformation preserves all local estimates, so that all our arguments still go 
through. 

• Arranging the right order of the pole of Q at the origin: As we saw in Chapter[6l 
the operator Q vanishes identically to degree five on the light cone. Thus the 
leading contribution to Q is of degree four on the light cone. Since Q always 
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involves a factor ^ (see (|7.2p ). the pole of Q at the origin is indeed of the 
required order o(|.^|~^). 

One might object that the near the origin x = y, where the continuum limit 
of Q{x, y) is not well-defined, the arguments of Chapter [6] do not apply, and 
thus there might be a non-zero contribution to Q which scales like (j6.8p . thus 
having a pole ~ I'^l"^- However, as explained after (j6.9p . we may assume that 
in the vacuum, the operator Q vanishes identically, even at the origin where 
the formalism of the continuum limit does not apply. Since to degree five, an 
interaction only leads to phase transformations (see (|6.27p ). the operator Q 
will then again vanish identically. As a consequence, the pole of Q will indeed 
scale like even without relying on the formalism of the continuum limit. 

• Handling nonlocal potentials: Proposition IA.2I does not apply to nonlocal po- 
tentials as introduced in §10.31 Another difficulty is that the support argu- 
ment used for Aij:^ will no longer apply. But as shown in Lemma IA.31 any 
contribution to Ai/^^ or A^^ which involves a nonlocal potential satisfies the 
inequality (jA.lSP and can thus be made arbitrarily small by choosing \Q\ suf- 
ficiently large. Following the arguments in the proof of Proposition IA.21 this 
gives us control of all error terms due to the nonlocal potentials, to every order 
in perturbation theory. 

We conclude that with the help of Proposition IA.2( we can justify the EL equations 
of the continuum limit (I5.29P for all fermion systems considered in this paper. 
We finally analyze the scalings in a universe of finite life time. 

Remark A. 4. (A universe of finite life-time) Suppose that instead of Minkowski 
space, we are considering a more realistic universe of finite life time tmax) like a cos- 
mology with a "big bang" and a "big crunch." In this case, the Fourier integral (I5.25P 
still gives a good local description of a Dirac sea (this is made precise in the example 
of a closed FRW geometry in [25^ Theorem 5.1]). However, one can no longer expect 
a continuum of states, and therefore the condition P^'i = can no longer be satis- 
fied by removing an arbitrarily thin strip around the mass shell. More precisely, the 
width Au of the strip in ()5.26p should be at least as large as the "coarseness" of the 
states in momentum space. This gives rise to the scaling (for details in the example of 
the closed FRW geometry see [25l Section 5]) 

Aoo ~ . 

^max 

The corresponding contribution to the Fourier integral (j5.26p scales as follows, 
A'^i := 'I'i(x) - r]{x) ~ sup |??| ^ ~ Aw 5 |??(0)| ~ — ^ l^i(0)l • 

" ^max 

As a consequence, the wave function no longer vanishes on £, but 

T I ^ IT I 

^l|£ ~ sup I Wi| . 

^max 

Furthermore, since A^'i is supported near the lower mass shell in momentum space, 
it decays in position space at infinity like the fundamental solution (pm — km){0,y), 
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smeared out on the scale 5. Combining these statements, we find that the correspond- 
ing contribution to the expectation value (j5.24p scales like 

<A'^i\Q ^'2> ~ sup l^fil sup |^'2| (5" — ^ e^~P . 

where p denotes the order of the pole at the origin being defined as the smallest 
integer p such that 

limsup(|e°| + lll)^"nr/(x,2/)| <oo. 
In comparison, the main contribution on the light cone around the origin scales like 

<^'i|Q ^'2> ~ sup l^'il sup |^'2| ^ ' 



and thus 



— e|ii| — . (A.lbj 



<^i|Q^2> t 

max ^ V ^ / £^max ^ 

This equation involves the fundamental length scale V^^max- The time since the big 
bang is estimated to about 13 billion years, which is the same order of magnitude as 
the size of the visible universe, estimated to 28 billion parsec. Thus it seems reasonable 
to assume that 

^max > 10^° years ~ 10^^ meters . 
Taking for e the Planck length e ~ 10~^^ meters, we obtain 



V £ /max ~ 10 meters . 

It is remarkable that this is about the length scale of macroscopic physics. Thus by 
choosing e\i}\ sufficiently small, we can make the quotient ()A.16P arbitrarily small 
without violating the scalings ()5.28p . We conclude that even if the life time of our 
universe is finite, this would have no effect on the statement of Proposition IA.2I 

Appendix B. Spectral Analysis of the Closed Chain 

In this appendix we analyze how different contributions to the fermionic projector 
effect the EL equations. In particular, we shall give the proofs of Lemmas 17.31 17.41 
17. 5( 17.61 and 17.71 Furthermore, we will analyze a pseudoscalar differential potential 
(see ()B.32p and ()B.33P ) and a scalar /pseudoscalar potential (Lemma IB. ip . 

We first explain the methods and the general procedure. The behavior of the 
fermionic projector near the light cone is described by the light-cone expansion (j4.17p . 
We concentrate on the singular behavior on the light cone as described by the se- 
ries in ()4.17p . disregarding the smooth non-causal contributions P^'^ and P^'^ (for the 
smooth contributions see Appendix |D] and ^8.3| also cf . the end of ^5.ip . The terms of 
this series can be computed as described in |I4t §2.5] (for more details see [I2| and [13]). 

The main task is to calculate the corresponding perturbation of the eigenvalues A_|_ , 
because then the effect on the EL equations is given by Lemma 17.11 In principle, 
the perturbation of the eigenvalues can be determined in a straightforward manner 
by substituting the summands of the light-cone expansion into the closed chain A^y 
(12. Sp . and by performing a standard perturbation calculation for the eigenvalues of the 
(4 X 4)-matrix A^y. However, the combinatorics of the tensor contractions inside the 
closed chain makes this direct approach so complicated that it is preferable to use a 
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more efficient method developed in jj^ Appendix G]. We now outline this method, 
giving at the same time a somewhat difi^erent viewpoint. 

Basically, the method is to work in a special basis of the spinors in which the 
calculations become as simple as possible. More precisely, we choose a double null 
spinor frame {f±^), being an eigenvector basis of the closed chain A^y of the vacuum 
to degree three. Thus following (j6.3p . we introduce the matrix 

2 



10 _ 5 
4 

According to ()6.6p . the corresponding spectral projectors in the formalism of the con 
tinuum limit are given by 

.0 1 



and they satisfy the relations 

FlU = zFl, and fU1=zF^_. (B.2) 

Furthermore, is invariant on the left- and right-handed components, and thus we 

may choose choose joint eigenvectors of the matrices Aq and 7^. This leads us to 

J //? 

introduce the four eigenvectors f_|_ by the relations 



Xci^f fs = f^ (B.3) 



with c € {L, -R} and s G {+,—}, which defines each of these vectors up to a constant. 
For clarity in notation, we write the inner product on Dirac spinors = ^'t^o^ 
as ^^l^*^, and refer to it as the spin scalar product. Then the calculation 

<f+ I f^^ = ^XLf+ I XLf+^ = -<f+ I XRXLf+^ = 

(and similarly for the other eigenvectors) shows that these vectors are indeed all null 
with respect to the spin scalar product. Moreover, taking the adjoint of (jB.ip with 
respect to the spin scalar product, one sees that 

{Fly = F° . (B.4) 

As a consequence, the inner products vanish unless the lower indices are different, for 
example 

I = ^F+f+ I i^+f?^ = -<f+ I F^Flf^y = . 
We conclude that all inner products between the basis vectors vanish except for 
the inner products ^f:^|f^^;^, -<f+|f:^^ as well as their complex conjugates -<f^|f+^ 
and -<f^\f^>-. We assume that all the non-vanishing inner products are equal to one. 



|^f^|f«^| = 1 = l^f^lf^^l 



(B.5) 



In order to specify the phases and relative scalings of the basis vectors, we introduce 
a space-like unit vector u which is orthogonal to ^ and ^. Then the imaginary vector 
V = iu satisfies the relations 

(f,^) = = , {v,v) = 1 and v = —v . (B.6) 

As a consequence, the operator ^ commutes with F^ and F^ , and since it flips parity, 
we may set f:^ = ^f:^. Next, from (jB.ip one derives the identities 

F^_( = iFl and F^_J = ( FI , (B.7) 
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which can be used as fohows. The first of these identities implies that the vectors ^ 
and are linearly dependent. The calculation 



(where we used (IB.Sh and (jB.ip ) shows that the vector ^ is indeed a real multiple 
of f:^^. Hence by normalizing appropriately, we can arrang^ that f;? = ^ f^. Using 
the second identity in ()B.7p , we also find that = ^ f + . Similarly, we may also 
set = ^ f + = ^f+- The resulting relations between our basis vectors are summarized 
in the following commutative diagram: 



- 




J 

- 


T 



(B.8) 



With (IB.3p . (IB.5P and (IB.SP we have introduced the double null spinor frame {f±^)- 
The construction involves the freedom in choosing the operator ^ according to (|B.6P : 
for given ^ the basis vectors are unique up to an irrelevant common phase. 

We next explain how we can represent a given linear operator B on the spinors in 
the double null frame (f^''^)- Following the notation in [141 Appendix G], we denote 
the matrix entry in the column (c, s) and row (c', s') by F^^'{B). These matrix entries 
are obtained by acting with B on the vector f^', and taking the inner product with the 
basis vector which is conjugate to f^, i.e. 



F-,(i?) = ^fl I Sf^,^ 



(B. 



where the conjugation flips the indices according to L R and + <-)•—. Similarly, we 
can also express the projectors XcFs in terms of the basis vectors, for example 



(B.IO) 



For computing (|B.9p . we use the relations in (|B.8j) to express the vector f^, in terms 
of f^, choosing the relations which do not involve factors of ^. Similarly, we express the 
vector in terms of f^^, avoiding factors of ^. Applying (jB.lOp . we can then rewrite the 
inner product as a trace involving the operator F?. More precisely, a straightforward 



Let us explain why we do not consider the opposite sign = To this end, we must show 

that -if+|^f+)- > 0. Since for any given positive or definite spinor ^, the vector X-L-P+C is a multiple 
of f+, it suffices to compute instead the sign of the combination XX-L-F+CI^Xi-P+C^- Applying HB.4|) 
and (|B.7|) . this inner product simplifies to -<C,\xRF-iC,'^- With the help of (|6.20|) and (|6.2|) . we can treat 
the factor ^ as an outer factor. Then our inner product simplifies to the expectation value ^ClXfli^C^- 
This expectation value is positive if we follow the convention introduced before (15. 7[) that > 0. 
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(B.ll) 



calculation yields 

F^^iB) = Tt{FIxlB) , F^^{B) = Tt{fI ^ xl B) 

F^^{B) = Tv{(FlixLB) , F^R{B) = Tv{^ F^ xl B) 

F^^{B) = i Tr(4 ^ ^XL B) , F^^{B) = i Tt{fI i xl B) 
Fht{B) = i Tr(^ Fl ixLB), Fh^{B) = i Tr(^ FIUxl B) 

(see also [141 equation G.19], where these relations are derived with a different method). 
Indeed, it suffices to compute the given eight matrix elements, because the other eight 
matrix elements are obtained by the replacements L -H- R. Moreover, the matrix 
elements of the adjoint (with respect to the spin scalar product) are obtained by 



(B*) = I B* = -<ff, I i? fi^ = F|£(i?) . 

Now we can describe our general procedure for the spectral analysis of the closed 
chain. We first perform the light-cone expansion of the fermionic projector. Then we 
calculate the matrix elements of the fermionic projector according to (jB.lip . Trans- 
forming to the double null spinor frame at such an early stage has the advantage that 
the contractions of the tensor indices (which come up by taking traces of products of 
Dirac matrices) are relatively easy to compute. After forming the closed chain Axy in 
our double-null spinor basis, we can compute the eigenvalues of A^y with a standard 
perturbation calculation (see \14:\ Appendix G.l] for a formulation with contour inte- 
grals). As the unperturbed operator we choose the closed chain (j6.25p which involves 
the axial phases. This is particularly convenient because the unperturbed operator is 
diagonal in our double null spinor basis, and moreover the unperturbed eigenvalues are 
non-degenerate according to (j6.27p . Thus it suffices to use simple perturbation theory 
without degeneracies. Next, it is useful that the unperturbed eigenvalues (16.270 form 
two complex conjugate pairs. This will remain true if perturbations of lower degree 
are taken into account, so that = A§. Therefore, it suffices to consider the eigen- 
value A^. The eigenvalue A;^ is then obtained by the replacement L ^ R, whereas the 
J //? 

eigenvalues A_ are obtained by complex conjugation. Expressing the perturbation 
calculation for A;^ in terms of the traces (jB.lip , one finds that to the considered degree 
on the light cone, the vector v drops out. 

To avoid calculation errors, the light-cone expansion was carried out with the help of 
the C-|— |- program class_commute, which was originally developed for the calculations 
in [12] and The traces in (jB.lip . which involve the contractions of tensor indices, 
are also computed with the help of class_commute. The resulting matrix elements of 
the fermionic projector are exported to the computer algebra program Mathematica 
(from this moment on, the tensor indices are simply treated as fixed text strings). The 
perturbation calculation as well as the expansions around the origin are then carried 
out by an algorithm implemented in Mathematica. This also has the advantage that 
the standard simplification algorithms of Mathematica and the comfortable front end 
are available. The C-|— |- program class_commute and its computational output as well 
as the Mathematica worksheets are available from the author on request. 

We now list the results of these calculations, also giving some intermediate steps. 
We note that some of these results are already obtained in [141 Appendix G.3], however 
without using the algorithm implemented in Mathematica. 
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Proof of Lemma \7.3[ Using for line integrals as in [13] the short notation (j6.24p and 

r[p,q\r]f:= [ {1 - {a - a'^Y f{ay + {1 - a)x) da , (B.12) 

Jx Jo 

the relevant contributions to the light-cone expansion can be written as (cf. |14t Ap- 
pendix G.3] and [HI Appendix A]) 

XLP(x,y) = ^XLe-^^''^T(-i) (B.13) 

-^XL^e.£[0,0|l]jlTW (B.14) 

+ lxL^ j'ril^l.T^'^ (B.15) 

r[0,l\0]F}^^jT^^^ (B.16) 

Jx 

r[0,l|l]^jir« (B.17) 

J X 



XL 
XL 



XL 



[0,2|0]jl7ir(i) (B.18) 



iXL 



y 



ii / y^^r(°) (B.19) 



X 

y 



+ V V. ^^^^ " ^^^^ ^^'^"^ 

+ irnxL^i r[0,0|l]yj)^r« (B.21) 

J X 

r[l,0|0]yF]j^7i7iTW (B.22) 

J X 



im 

-^XL 



2 XL I [G,mF]^^ajYT^^^ (B.23) 



+ imxL j [0, 1 1 0] \Y{djA^^) - {djA{) yj T^^) (B.24) 

+ '^XLHiljho\ 0] yy^5, r(o) (B.25) 

+ ^ XL ^ £[0, 1 1 0] A'lYY r(0) (B.26) 

- XL r% I 1] yyF^^j T(i) (B.27) 

- XL Ho, 2 I 0] Fl^^j YY T^^) (B.28) 

J X 

+ m^XL r[i,o|o]yy4^rW (B.29) 

J X 

^[o,o|o]y4^yr(i) (B.30) 
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+ XL 



/ [0, 1 1 0] 4^ yy r(^) (B.31) 

J X 

i (deg < 1) + (deg < 0) + 0(^i/K) , 



where fI^ = d^A^ — d^Ai is the chiral field tensor and = d'jAi — OA^ is the 
corresponding chiral current. The term ()B.13P is our unperturbed fermionic projec- 
tor ()6.23p : all other summands are our perturbation. We take the partial trace, expand 
in powers of ^ and compute the matrix elements (jB.lip . In these computations, the 
error terms of the form ^ (deg < 1) are contracted with another factor ^ or ^, giving 
rise to a term of lower degree. Likewise, the higher orders in A either give rise to terms 
either of lower degree on the light cone or of higher order in the expansion around the 
origin. Computing AX^ by a first order perturbation calculation gives 



.2 , , — r 0^2 



AxL .i (1) (-1) ig i (0) io) 

^■^+ 3l [0] [0] ^» [0] [0] 



t^'M';'Ht^!>M^ + C.C. ]-c.c. 



- 2im At. fi Y =L= h (deg < 2) + o( £ ) . 

T^(o)^(-i) _ ^(-1)^(0) ^ ^ ' ' 

The other eigenvalues are obtained by the replacement L ^ R and by complex conju- 
gation. Substituting the resulting formulas into (jT.ip gives the result for the current 
terms. 

For the computation of the mass term, we must consider another contribution where 
the mass expansion of the fermionic projector of the vacuum 

P(x, y) X mY t(°) + ^ rr?Y^ i r(°) 
is taken into account in a perturbation calculation to second order. 



AA^ X y2 r^'J^r^'J) + gm^ YY T^JJ^T^';) 

\L-\R V [1] M M [1] ; V [0] [1] [1] [0] 

The analogous formulas for the other eigenvalues are obtained by the replacements L o 
R and -H- 17 as well as by complex conjugation. In order to compute |AA:^|, we use 
the identity 

a|a^:|x^(a^aa^: + a^:aa^) 

and express the unperturbed eigenvalues in terms of (j6.27p . Expanding the phase 
factors V and u in powers of the axial potential (see (j6.26p ) and keeping the linear 
term in j4a, we obtain 

lA+l [0] ^[0] Vn ^[2] ^[2] ^[0] 
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and similarly for the other eigenvalues. Substituting these formulas into (j7.ip com- 
pletes the proof. □ 



Proof of LemmaVrA According to (|7.8|) and (|4.1U|) . the perturbation AP(a;, y) by the 



Dirac current satisfies the relations 



and thus 



Tr(7^ AP(x,y)) =-^4, Tr (7V A^^, y)) = ^' 



AP(x,x) = -^7,4 + ^7S,J^ 



The corresponding perturbation of the eigenvalue AA^l is computed to be 



The formula for TZ follows by direct calculation. □ 

We now come to the pseudoscalar differential potential (j7.14p . The contribution to 
the fermionic projector linear in v has the form 

P{x, y) X i 75^ J' ^v' r(-i) + 7' C^ v\x) r(-^) + (deg < 2) . (B.32) 

After applying the relation 2^^v^ = 2S^^djV^ + [^,(^v^], we can in the first term inte- 
grate by parts to obtain ()7.15p . The light-cone expansion to lower degree involves 
many terms, which we do not want to give here. To higher order in the mass, the 
contributions become less singular on the light cone. In particular, the leading term 
cubic in the mass takes the form 

3 ry 

P(x, 2/) X - ^ 7^^ / [0, 1 1 0] (yv'YY - v'YYY^ t(°) 
l^i y"^[l, I 0] (yYYv' - YYv'Y^ r(°) 



2 



+ im^-i^ /^[0, 1 1 0] (Yi^YY - jYYY\ T^^) 

+ im^-f^ y^[i, 1 0] (yyy^ - yy^y) r(^) 



(B.33) 



+ i (deg < 1) + (deg < 0) . 

Expanding in powers of ^, we obtain ()7.16p and (I7.17p . 

We point out that for the scalar differential potential, the higher orders in pertur- 
bation theory are difficult to handle because they are not of lower degree on the light 
cone. Moreover, a resummation procedure similar to that for chiral potential does not 
seem to work. For a constant potential, this problem corresponds to the effect of the 
"deformation of the light cone" as discussed after (j7.20p . In the more general setting 
here, this problem means that the scalar differential potential cannot be treated per- 
turbatively in a convincing way. This serious difficulty was our original motivation for 
introducing the vector differential potential (I7.22[) . and to rewrite the combination of 
these potentials by the local axial transformation (j7.24p . 

For the local axial transformation ()7.24p . the fermionic projector can easily be com- 
puted non-perturbatively. To first order in v, we obtain the contribution (j7.30p (for a 
non-perturbative treatment see Appendix [C]) . 
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Proof of Lemma 7.5, Since the commutator in (j7.30p vanishes if Pj^ is a multiple of 
the identity matrix, the contribution by the parameters in (j7.30p is the term 

g 

P{x, y) X 7^ v,e E (^[0] + "^1 ^[2?) + (deg < 1) • (B.34) 

/3=1 

Denoting this term by AP(x,y), we must compute the corresponding first order con- 
tribution to the eigenvalues A^. For the unperturbed fermionic projector PQ[x,y) we 
now choose the fermionic projector of the vacuum (|6.ip . The first order perturbation 
calculation involves traces of products of the matrices PQ{x.,y)^ Po{x,y)* = Po{y,x), 
AP{x, y) and (AP(x, y))* ■ More precisely, omitting the arguments (x, y), to first order 
in AP we obtain terms of the form 

Tr((AP)Po*---i^oi^o) and Tr((AP)*Po • • • ^'o ^'o) , (B.35) 

which involve one factor AP or (AP)*, and where the matrices and their adjoints 
alternate. According to ()B.34p . the matrix AP (and similarly (AP)*) is a multiple 
of 7^. But Pq and Pq only involve the Dirac matrices <f and ^. As a consequence, the 
traces in (|B.35j) vanish. □ 

Proof of Lemma 7.6, In view of ()7.33p . the contribution of the parameters djs in (I7.30p 
to the fermionic projector is given by 

a 

P(x, y) X 2 E d/s 7^ ^^"^ + ^ (deg < 2) + (deg < 1) . 

13=1 

The corresponding perturbation of the eigenvalue AX^ is computed by 



AA^ X 2ig v,i^ Yl "^/^ dp T^^^Tjj^j + (deg < 3) . 

/3=1 

The other eigenvalues are again obtained by the replacement L o P and by complex 
conjugation. Substituting the resulting formulas into (|7.ip gives the result. □ 

Proof of Lemma \7. 7\ In view of (|7.33p and (|7.34p , the contribution of the parame- 
ters dp in (|7.30p to the fermionic projector is 

a a 
P{x, y)^'-Y,mldp 7'[^, i] + 2 ^ 7V T^^^ + (deg < 0) . 

/3=1 /3=1 

The corresponding perturbation of the eigenvalue AA^ is given by 

a 



/9=1 



+ 2i Vj^^ mY rrip d^ 



(-1) (0) (0)^(0) _ 
^™ ^[0] -'[I] ^[2] 



— 7.{0)^{-l) _ ^(-1)7.(0) 



(0) (-1) /^(O)^(O) _ ^(0)^(0) 
+ 2i VjH mY > rrif, da =^ + (deg < 2) . 
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The other eigenvalues are again obtained by the replacement L ■H' R and by complex 
conjugation. Substituting the resulting formulas into (jT.ip gives the result. □ 

We finally consider the perturbation of the fermionic projector by the scalar and 
pseudoscalar Dirac current (j9.ip . According to (j4.10p , the corresponding perturbation 
of the fermionic projector is given by 

AP(x, y) = -i- ( J, + ij^ J,) + o(|ef ) . (B.36) 

Lemma B.l. The first order contribution of the perturbation ()B.36p to the opera- 
tor Q{x, y) is of degree two on the light cone. 

Proof A first order perturbation calculation yields 



+ 47r T.{0)^{-1) _ ^.(-1)^(0) 

Note that the pseudoscalar current dropped out. This cancellation can be understood 
exactly as in the proof of Lemma 17.51 by considering the traces (|B.35p and noting that 
the corresponding AP is again a multiple of the matrix 7^. 

The first order contribution to the EL equations is obtained by considering first 
variations of the Lagrangian in the from (I6.29P and substituting the formulas for AA^. 
Counting degrees, this contribution is expected to be of degree three on the light cone. 
However, since ()B.37p has even parity, whereas the first variation of the Lagrangian 
involves factors (Af — |Af ) with odd parity, this expected contribution of degree three 
vanishes. □ 



Appendix C. The Local Axial Transformation to Higher Order 

We shall now analyze the local axial transformation of the fermionic projector of the 
vacuum, taking into account the higher orders in v which were disregarded after (I7.27p . 
In order to understand why these higher orders in v are a potential problem, one should 
recall that in §7.61 we satisfied the EL equations to degree five by a special choice 
of the local axial transformation (see (|7.34p ). Then the local axial transformation 
still affected the EL equations to degree four, making it possible to compensate the 
logarithmic poles on the light cone (see Lemma l7.7p . If the higher orders in v gave 
contributions to the EL equations to degree five, these contributions would be by an 
order larger than the contribution computed in Lemma |7.7| and would thus not 
be negligible (even if the expansion parameter qv were very small). 

In order to rule out this potential problem, we must show that the EL equations to 
degree five can be respected by the local axial transformation even non-perturbatively. 
This is done in the following theorem. 

Theorem C.l. Suppose that the number of generations g > 3. Then for any vector 
field u{x) satisfying the condition 

{u(x),u(x)) > -—— max (ruq - rna)'^(rna - rniS^irni + rn^ + rn„)'^ (C.l) 

25d ae{2,...,g-l} 

(where the masses are in increasing order (|3.2p ) there is a local axial transformation of 
the form ()7.24p such that at every space-time point x and for all Minkowski vectors ^ 
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the following relations hold: 




constants c\ and C3 may depend on x, and cq may depend on x and ^. 

The relations ()C.2p and (|C.3P imply that the local axial transformation does not affect 
the EL equations to degree five. The contribution to the fermionic projector corre- 
sponding to (|C.4p makes it possible to compensate logarithmic poles on the light cone, 
(exactly as explained in ^8.ip . The inequality (jC.ip implies that the vector field ap- 
pearing in the logarithmic pole must not be spacelike and too large. However, the 
constraint (jC.ip does not seem to be physically relevant because under realistic condi- 
tions, the vector fields rn^Aa and ja are very small on the scale ~ . For notational 
convenience, the matrix m'^U^Y'^U~^ , which is disregarded in the above theorem al- 
though it enters the EL equations to degree four (cf. the contribution (I7.35P ). will be 
analyzed separately in Proposition IC.6I below. 

The remainder of this appendix is devoted to the proof of this theorem and to the 
statement and proof of Proposition [C^ First, since all constructions will be performed 
locally at a given space-time point x, we may fix x throughout and omit the space- 
time dependence. Moreover, we always choose the vector v in (I7.24p as a multiple 
of the vector u which appears in the statement of the theorem. Since in (I7.24p we 
may multiply g by a positive number and divide v by the same number, it suffices 
to consider the three cases = 0, = 1 and v'^ = —1, where the vector v is null, 
timelike or spacelike, respectively. In the case v'^ = 0, the higher orders in v vanish, so 
that (j7.27p and the constructions thereafter become exact. Thus it remains to consider 
the cases = 1 and = —1. According to ()7.25p . the local axial transformation in 
these two cases involves hyperbolic and trigonometric functions, respectively, 

[/ = coshg — ^7^^ sinhg , C/^^ = coshg + 17^^ sinhg if = 1 (C5) 

[/ = cosg — ^7^^ sing , [/^^ = cos g + ^7^^ sing if = — 1 . (C.6) 

In view of this difference, it is preferable to treat the cases v'^ = 1 and = —1 sepa- 
rately. 

First case: = 1. The first step is to rewrite the conditions (IC.2P and (IC.3P in 
terms oi gx ^-matrices acting on the generations. The condition (1C.2P should hold for 
any and by linearity it suffices to consider the two cases v) = and ^ = v. In 
the first case, the factor ^ commutes with U . Thus U^U^^ = ^, and ()C.2p is trivially 
satisfied. In the second case, the factors ^ and 7^^ anti-commute, and using ()C.5P we 
find that U^U~^ = U"^^- Again using (fasjl . the condition ([02]) can be written as 



(C.7) 



Similarly, the relations (jC.SP and ()C.4p can be restated as 



UYU~^ - U'^YU = 



(C.8) 
(C.9) 
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Our procedure is to first satisfy ()C.7p and (jC.SP : then we will analyze for which vec- 
tors u we can arrange ()C.9P . The main advantage of ()C.7p and (jC.SP over the equivalent 
conditions ([02]) and ([03]) is that in ([07]) and ([OSP only the factor ^^^^ in [7 and C/"^ 
involves Dirac spinors. Thus diagonalizing this factor and restricting attention to the 
respective eigenspaces, we obtain matrix conditions on C^. More precisely, using the 
relation (17^^)^ = 1, one readily sees that the matrix 17^^ has eigenvalues ±1. We 
can thus rewrite ([07p and ([OSP in the equivalent form 

{l\V^l) = {l\V-^l) , H\VYV-'^l) = {l\V-^YVl) , (C.IO) 

where we introduced the matrix V and the vector [ by 

y = coshs + sinhg = e0 G Mat(C^), [ = (1, . . . , 1) G , (Oil) 

and (.|.) denotes the standard scalar product on C^. Note that the matrix V is positive 
definite. In what follows, we can work only with this property of V, because the 
representation ([Clip can then be obtained by setting g = logV. Next, we want to 
exploit the fact that only certain combinations of the matrix elements of V and 
enter ([C.lOp . This can be achieved by setting 

m = V-^l and n = Vl, (C.12) 

giving the equivalent conditions 



(n|n) = (m|m), {n\Ym) = {m\Yn) , (013) 



which instead of matrices merely involve the two unknown vectors m, n G C^. 

Before we can work exclusively with the vectors I, m and n, we need to work out geo- 
metric conditions which allow the representation ([012p . In order to avoid trivialities, 
we assume that the vectors [, m and n are linearly independent. 

Lemma C.2. For given linearly independent vectors I, m, n G there exists a positive 
definite matrix V satisfying (I012P if and only if the following relations hold: 

(m|n) = ||[f (014) 

(m|[)>0, ([|n)>0 (015) 

(m|[)([|n)>||[r. (016) 

Proof We first verify that the relations ([OT4P - ([OT6]) are necessary. Thus assume 
that (|C.12p holds for a positive definite matrix V. Then (|C.14p follows immediately 
from the calculation 

||[||2 = {[\VV-H) = {Vl\V-^l) = (n|m) . 
The inequalities ([C.ISP are readily obtained from the estimates 

< = ([|n) and < ([ | F^^l) = | m) . 

We now form the matrix of expectation values of V with respect to the vectors [ and m, 

{l\Vl) {l\Vm)\ _ f{l\n) 



{m\Vm)) (m|[)^ 

Since V is positive definite, the determinant of the last matrix must be positive, giv- 
ing dOTH]) . 

It remains to show that the relations ([OT4]) - ([OT6]) are sufficient. Thus assuming 
that these relations hold, we must construct a corresponding positive definite matrix V 
satisfying (|C.12p . It suffices to construct a positive definite matrix V on the subspace 
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spanned by [, m and n, because extending V on the orthogonal complement of this 
subspace by the identity gives the desired matrix V. Since we want to satisfy ()C.12p . 
the corresponding matrix of expectation values must be of the form 

{l\Vl) {l\Vm) {l\Vn)\ /([|n) (n|n)\ 
{m\Vi) {m\Vm) (m|yn) = (m|n) (m|[) (l|n) , (C.17) 
(n I VI) (n I Vm) (n \Vn) J \{n\n) (n 1 1) C J 

involving only one unknown parameter C := (n | Fn) > 0. According to (IC.lSh 
and (fcn6]l . the upper left 2x2 submatrix is positive definite. Thus by choosing C suf- 
ficiently large, we can arrange that the whole matrix in ()C.17|) is positive. Then (1C.17P 
defines a positive definite matrix V with the required properties. □ 

It remains to construct vectors m and n satisfying the relations ()C.13p - (|C.16p . We 
first note that the first relation in ()C.13P can be arranged by the scale transformation 

tn— ;>Am, n ^ with A > (C.18) 
A 

without affecting all the other relations. Thus the first relation in (jC.lSP can be 
disregarded. Moreover, the first inequality in (jC.lSP can be arranged by the phase 
transformation 

m^e^'^m, n^e^'^n with G [0, 27r) , (C.19) 

leaving all the relations except for ()C.15P unchanged. If the first relation in (jC.lSP 
holds, the second follows immediately from (]C.16p . Thus we may also disregard (jC.lSp . 
leaving us with the set of conditions 

(n|ym) = (m|yn) , {m\n) = g (C.20) 
(m|[)([|n) >ff2. (C.21) 

Introducing the real parameters and by 

n^mp = + 2id/3 , (C.22) 

we can write the relations (|C.20p as 

a 9 9 

13=1 13=1 13=1 

This notation gives a close connection to the analysis in ^7.61 because in a first order 
expansion in v one easily verifies that = 1, whereas dp coincides with the corre- 
sponding parameters in (j7.29p . Hence the last two equations in ()C.23p correspond 
to (|7.33p and (j7.34p . whereas the first equation in ()C.23p gives an additional condition 
to be fulfilled to higher order in v. 

We now give a procedure to satisfy the condition ()C.2ip . 

Lemma C.3. For a given solution of (|C.23p . there are vectors m and n which satisfy 
both (IU:20D and (fCTl . 

Proof. For given parameters Tp and dp satisfying ()C.23p . we can fulfill the condi- 
tions ()C.20p by choosing any complex solutions m.p and np of ()C.22p . This choice 
involves the freedom of rescaling the parameters and np by 

mp ^ zpmp, np^ (z^)"^ np , with zpeC\{{)] . (C.24) 

We want to use this freedom to arrange ()C.2ip . 
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We begin with the case that one of the components of m or n vanishes; by symmetry 
we can assume that mi = 0. Furthermore, we may assume that ni 7^ 0, because 
otherwise we can remove the first component and repeat the proof for g replaced 
by 5 — 1. Then by a suitable choice of the parameter zi in (10241) close to zero, we 
can make the component txi positive and arbitrarily large. Thus ([|n) can be made 
arbitrarily large without affecting (f|m). Next, by a suitable choice of the phases of 
the remaining parameters Z2, ■ ■ ■ ,Zg, we can arrange that the components of m are all 
positive, so that (I | m) > 0. In this way, we make the product (m | [) ([ | n) as large as 
we like, thus satisfying ()C.2ip . 

It remains to consider the case that the components of m and n are all non-zero. 
In this case, we can use the transformation ()C.24p to arrange that m/3 = 1 for all /? = 
1, . . . ,g. Then m = I and thus 

{m\l) H\n) = g {m\n) = , 

where in the last step we used the relation on the right of (lC.20p . This is not quite 
good enough, because in (IC.2ip the strict inequality appears. Therefore, we now use 
the transformation (lC.24j) to vary m. Considering first order variations, we have 

m/3 = 1 + (5m/3 , 6ni3 = -5m/3 

and thus 



S{{xn\l){l\n)]=Y,Sra^ 

13=1 



a 

^na-gn,3 

0=1 



(C.25) 



The square bracket vanishes identically only if the components are all equal. In 
view of ()C.22p and ()C.23p . this corresponds to the trivial case 

Tg = 1 and = for all /3 = 1 , . . . , 5 , 

which can clearly be arranged by choosing V = 1. In all other cases, we can choose Jm^ 
such that (1025]) becomes positive. Having shown that ()C.2ip is satisfied by a suitable 
linear perturbation, this condition clearly also holds for a corresponding nonlinear 
perturbation, provided that this perturbation is sufficiently small. □ 

We have thus reduced the problem to the analysis of the relations ()C.23p . Clearly, 
these equations have non-trivial solutions, and by a suitable choice of the parame- 
ters di3, we can give the combination 

g 

((n I Y^m) - (m | Y^n)) = 4i ^ m| df} (C.26) 

/3=1 

an arbitrary imaginary value. Using this fact in ()C.9p . we see that the vector u can be 
an arbitrary multiple of v. This concludes the proof of Theorem lC.il in the case v"^ = 1. 

Second case: v'^ = —1. Proceeding as in the first case, we can rewrite the condi- 
tions ()C.2p and ()C.3P again in the form (jC.lOp . where V now is the unitary matrix 

y = cosB + ising = . (C.27) 

Again introducing the vectors [ = (1, . . . , 1) and m, n by ()C.12p . we can write the 
conditions (jC.lOp in analogy to (IC.ISP as 



(m|n) = (n|m), {m\Ym) = {n\Yn) . (C.28) 



The next lemma is the analog of Lemma IC. 21 in the case that V is unitary. 
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Lemma C.4. For given linearly independent vectors [, m, n G there exists a unitary 
matrix V satisfying ()C.12p if and only if the following relations hold: 

||[|| = ||m|| = ||n|| and (m|[) = ([|n). (C.29) 

Proof. The necessity of the relations (|C.29p is obvious from the calculations 

(m|m) = {V-H\V-H) = = {Vl\Vl) = (n|n) 

(m| [) = {U-H\l) = H\Ul) = ([|n) . 

In order to show that these conditions are also sufficient, we make the ansatz 

V = V2V1 

with two unitary transformations Vi and V2. We choose Vi such that it maps m to I 
(this is possible, because according to (|C.29P both vectors have the same norm). Then 
the unitary transformation V2 should map the vector [ := to n, while leaving the 
vector [ invariant. Using the right equation in ()C.29p . we find 

([|[) = ([|Ffi[) = (I|m) = (n|[), 

showing that the the vectors I and n have the same projection on [. According to the 
left equation in ()C.29p . they also have the same norm. Thus we can choose V2 as a 
unitary transformation of the orthogonal complement of [ which maps [ to n. □ 

The remaining task is to construct vectors m and n satisfying the conditions ()C.28P 
and ()C.29p . The relation on the left of ()C.28P can be arranged by the phase transfor- 
mation 

m^e^'^m, n-fe'^'^n with G [0, 27r) , (C.30) 
leaving all other conditions unchanged. Introducing the parameters and dp by 

rp = \{\M' + \mp\') and dp = \{\np\'' - \rap\') , (C.31) 

the conditions on the right of ()C.28P and the left of ()C.29P again take the form ()C.23p . 
but now with the obvious additional constraint 

\dp\ < ^ (C.32) 

(again, our definition of dp agrees with (j7.29p to first order in v). The condition on 
the right of ()C.29P can be arranged by the following construction. 

Lemma C.5. Suppose that the parameters Tp and dp satisfy the relations (IC.23P 
and (|C.32p . Then there are vectors m and n which fulfill all conditions in (IC.28P 
and (ICIMI) . 

Proof. For given parameters Tp and dp satisfying ()C.23p and (lC.32p . we choose solu- 
tions mp and np of (IC.3ip . This choice involves the freedom to change the phases 
according to 

m;3 e'^ff mp, np ^ e'^P np with ipp,^p G [0, 27r) . (C.33) 

Our goal is to show that by choosing these phases appropriately, we can satisfy the 
condition on the right of ()C.29p . In order to get a more convenient notation, we take 
the absolute values of the components |m^| and |n^| and bring them in increasing 
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order to obtain the parameters pi, . . . , p2f- Then our task is to show that there are 
angles (pi, . . . ,(j)2f such that 

2/ 

^e''^^Pf, = 0. (C.34) 

13=1 

The relations on the left of ()C.29p imply that X]^=i p"^ ~ ^^'^ moreover we know 
that there is a partition of the pfs into two subsets such that the sum of each subset 
equals g. This implies that 

P2/<Pi + ---P2/-i- 

Taking the square root and inductively applying the inequality Vab < i/a + Vb (with 
a,b > 0), we find that 

P2f - P2f-i < Pi-\ h P2/-2 • (C.35) 

We claim that this inequality allows us to choose sign factors £ /3 = 

1, ... ,2/ - 2, such that 

Z := Sipi H h S2/-2 P2/-2 G [P2f - P2f-1, P2f + P2f-l] ■ (C.36) 

This claim can be verified inductively as follows. We choose si = 1, S2 = 1, etc. 
until the partial sum S{k) := sipi + • • • + s^Pk exceeds P2f — P2f-i (which necessarily 
happens according to ()C.35P ). Then due to the increasing ordering of the p^, the 
partial sum lies in the required interval in ()C.36p . Moreover, we can choose the 
following sign factors s^+i, ■ ■ ■ , S2/-2 inductively such that the partial sums stay inside 
this interval. This gives the claim. 

Using ()C.36p . it remains to show that there are angles ip2f and (p2f-i such that 

e''P^fp2f + e''P^f-'p2f-i = -z. 
This verified by an elementary consideration, proving (IC.34p . □ 

It remains to analyze the conditions (1C.23P under the constraint (IC.32p . Obviously, 
there are non-trivial solutions. We must show that in (IC.4P we can realize every 
space- like vector u which satisfies (IC.ip . Evaluating (IC.9P on the eigenspaces of the 
operator i^^-^, we find that 

u = -m^ (([ I yy V"H) - ([ I V-^Y^Vl)) V . 

Using ()C.12p together with the fact that V is unitary, we conclude that 

u = ^6v, (C.37) 

where 

9 

6 := - W ((m I y^m) - (n | Y^n)) = ^ m| , 

/9=1 

and in the last step we applied ()C.3ip . We need to determine the possible values 
of the functional 6. Since multiplying the parameters by a number of modulus 
smaller than one preserves the conditions ()C.23P and ()C.32p . it is obvious that the 
possible values of & form an interval [—©max, ©max], and thus it suffices to compute 



AN ACTION PRINCIPLE FOR A FERMION SYSTEM AND ITS CONTINUUM LIMIT 105 



the maximal value ©max- We first note that we can eliminate the parameters by 
rewriting the constraint ()C.32p as 

E|d/3|<f. (C.38) 

13=1 

Moreover, since scaling the parameters djs by 



P ' 9 



— 7- — with p = - y \d(^\ < 1 



a=l 



changes the inequality in ()C.38P to an equality and at most increases (3, we see that 
for finding the maximum of (3 we can assume that equality holds in (jC.SSp . Thus our 
task is to maximize (3 under the constraints 

9 9 9 

^^^ = 0, ^m/3(i^ = and ^|(i^| = -. (C.39) 

(3=1 13=1 13=1 

Treating the constraints with Lagrange multipliers A, fj, and k, we are led to searching 
for critical points of the functional 

9 r 9 

E (m| - - /x) d/3 - K ^ {df^l - - 

13=1 ^13=1 

Taking the distributional derivative with respect to da, we obtain the conditions 

— A iria — fi = K if da > 

w?a - Xnia - p- = - K if da < (C.40) 

\'W?a — Xrua — if do = . 

In order to satisfy the first two conditions in ()C.39p . there must be indices ai < 02 < 0:3 
such that da^ and have the opposite sign as da2- This is compatible with ()C.40p 
only if A > 0, so that m^^ is on the decreasing and niao, on the increasing branch of 
the function nv' — Am. The last condition in ()C.40p implies that 

di,dg > and d2, ■ ■ ■ , dg-i < . 

Comparing the first equation in ()C.40p for mi and nig allows us to determine A, 

A= ^ V (C.41) 

mg — mi 

Next, we multiply the equations in ()C.40p by da and sum over a. Using the con- 
straints ([a39]l . we obtain 

6 = Y • (C.42) 

Thus to maximize S, we must make k as large as possible. Since at least one of the 
parameters d2, ■ ■ ■ , c?/-i must be negative, we find from (IC.40p that 



2k < max 

ae{2 



,x ({ml - Ami) - {iria - Am^) ) . (C.43) 
,9-1} ^ ^ 

This upper bound of k can indeed be realized by solving ()C.39P for a configuration 
where only three parameters djs are non-zero: di,dg > and da < for an a where 
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the maximum in ()C.43P is attained. A short calculation using ()C.4ip and ()C.42p shows 
that 

©max = 7 max {mg-ma){ma-mi){mi+ma + mg) . 

4 oe{2,...,g-l} 

Substituting this result into (|C.37p gives (jC.ip . This concludes the proof of Theo- 
rem [C]TJ 

We note that our method of proof is constructive in the sense that for a given 
vector u in ()C.4p we could compute the corresponding transformation U explicitly. 
More precisely, if u is timelike, we set v = u/ {u, u) to be in the first case v'^ = 1. 
We choose parameters and which satisfy ()C.23p . Then the construction of 
Lemma IC.31 gives us parameters rri/j and which fulfill (|C.20p and (|C.2ip . Next, the 
transformations (jC.lSp and (|C.19p allow us to satisfy the conditions (|C.13p - ()C.16p . 
The corresponding positive definite operator V can then be chosen according to (1C.17P 
for a sufficiently large parameter C. According to the definition of V, (jC.lip . we know 
that coshg = {V + V~^)/2 and sinhg = {V — V~^)/2, and substituting these formulas 
into (|C.5p gives U. Likewise, if u is spacelike, we set v = u/y/\{u \ u)\ to be in the 
second case = —1. We choose parameters and dfs which satisfy (|C.23p as well 
as the constraint (|C.32p . The construction of Lemma IC.5I gives us solutions of (jC.28p 
and ()C.29p . Then the unitary transformation V can be constructed as in Lemma IC.4i 
According to ([QT]) . we know that cosg = {V + V'^)!'! and sing = - F"^)/(2i), 
and using these relations in ()C.6P gives \J . 

We finally analyze the matrix C/^y^f7~^. 

Proposition C.6. Using the above notation, 

9 

m^U^Y^U-^ = C2i mHy + C2 ij^ (C, - 7' ^] ^ ^ , (C.44) 

/3=1 

where the constants C2 and C2 depend on x and ^. 

Proof. Since the term involving the inner product v) in (IC.44P involves a free con- 
stant, we may assume that v is orthogonal to Then we can rewrite the claim ()C.44p 
similar to (jClSl) and (IClOl) as 

g 
/3=1 

Both sides of this equation commute with the operator ^7^^, and thus we may again 
restrict attention to the respective eigenspaces. We thus obtain the equivalent relation 

g 

13=1 

where we again compared ()C.5p and (|C.6p with (IC.llh and (|C.27j) . This relation 
is verified in the case t;^ = 1 by using that V is Hermitian and applying ()C.12p 
and (IC.22P , whereas in the case = —1 we use that V is unitary and apply ()C.12p as 
well as (IClMll . □ 

The point of interest is that the contribution which enters the field equation can 
again be expressed in terms of the parameters dp. We conclude that the analysis 
in ^8.11 remains valid without changes even if the local axial transformation is treated 



J Ai if = 1 
1 -4 liv"^ = -1 , 



AN ACTION PRINCIPLE FOR A FERMION SYSTEM AND ITS CONTINUUM LIMIT 107 



non-perturbatively. In particular, we learn that, in contrast to what one might have 
expected naively, the higher orders in qv do not give rise to higher order corrections 
to the field equations. 

Appendix D. Resummation of the Current and Mass Terms at the Origin 

As pointed out in ^4.41 the distribution Ta is not a power series in a, and thus it 
cannot be expanded in a Taylor series around a = (see (j4.13p and the explanation 
thereafter). The method of subtracting suitable counter terms (|4.14p has the short- 
coming that the subsequent calculations are valid only modulo smooth contributions on 
the light cone. This method is suitable for analyzing the singularities on the light cone, 
but it is not sufficient when smooth contributions to the fermionic projector become 
important (cf. the discussion after (jS.lip and the beginning of ^8.ip . We now present a 
convenient method for computing the smooth contributions to the fermionic projector. 
Our method is based on the resummation technique developed in [12^ Section 4] and 
is outlined as follows. We first perform the mass expansion not around zero mass, but 
around a given mass parameter a > 0. Then according to ()4.13p . the distribution Ta 
is smooth in a, and we may set 

T<") = (A)"r„. (D.i) 

Adapting the method of the light-cone expansion, we can express any Feynman tree 
diagram as a sum of terms of the form 

oo 

P^'^^{x,y) = ^ mP'' (phase-inserted nested line integrals) x T^"'\x,y) , (D.2) 

n=— 1 k 

where for each n, the fc-sum is finite, whereas the n-sum is to be understood as a formal 
power series. Note that, in contrast to the series in (|4.17p . the infinite sum in ()D.2p is 

not a light-cone expansion in the sense of Definition 14. H because the distributions ri"^ 
all involve smooth contributions and are thus only of the order 0{{y — x)^). We proceed 
by partially carrying out the series in ()D.2p , giving rise to explicit smooth contributions 
on the light cone. After this resummation has been performed, we recover (j4.17p . but 
now with an explicit formula for P^'^{x,y). 

For simplicity, we develop the method only for the contribution to the fermionic 
projector needed here: the vector and axial components of the fermionic projector 
perturbed by chiral potentials to first order. But the method generalizes in a straight- 
forward way to arbitrary Feynman tree diagrams. Furthermore, we begin by con- 
sidering a single Dirac sea (the generalization to several generalizations will then be 
straightforward; see the proof of Lemma lS.ll below). We thus consider the contribution 
to the fermionic projector 

AP = -SmiXL^R + XR4-L)'^m - t^iiXL^R + XiJ^L)*^ (D-3) 

with the spectral projector tm and the Green's function Sm as in ()8.26|) . In order 
to concentrate on the vector and axial components, we want to consider the expres- 
sion Tt{^xl AP[x,y)), being a well-defined distribution. The singular part of this 
distribution on the light cone can be computed by inserting the formulas of the light- 
cone expansion (]B.13j) - (lB.3ip and using the contraction rule 

i^T^'^\x,y) = -4nr("+^)(x,y) + (smooth contribution) , n G {-1,0} . (D.4) 
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(which is immediately verified from the exphcit formulas (|4.13p - ()4.16p l. We thus 
obtain 

\Tv{ixL AP{x, y)) = 2j\k At TjiJ) {x, y) (D.5) 

- 2 j\a -a^)ik A (x, y) + j' [a^ - r« (x, y) (D.6) 

+ e'/fc(x,y) + (deg<0), (D.7) 

where we added subscripts [.] in order to indicate how these factors are to be reg- 
ularized (although we do not need a regularization at this point), and f{x,y) is a 
yet undetermined smooth function (clearly, the summand in ()D.5p is the gauge term 
as discussed in §6.21 whereas the summands in (ID.6P correspond to the current and 
mass terms considered in ^7.ip . Our goal is to compute the functions fk{x,y) at the 
origin x = y. 

Our first step is to perform a mass expansion of the Feynman diagram (ID.3P around 
a given a ^ 0. To this end, we need suitable calculation rules which are derived in the 
next lemma. 

Lemma D.l. The distributions (ID.ip satisfy for all n G No the calculation rules 
(-a - a)rW(x,y) = nTi"-i)(x,y) (D.8) 
ATi"+i)(x,y) = iefcrW(x,y) (D.9) 
e2ri")(x,y) = -4nri-+i)(x,y) -4ari"+2)(^^y) . (d.iq) 

In the case n = -1, the rule i^PM can be used to define the distribution ^^ri ^\ Using 
this definition, the rule (jP.lOp also holds in the case n = — l. 

Proof. The relations (jP.Sp and (|D.9p were already derived in [13] (see [131 equa- 
tions (3.5) and (3.6)]). For self-consistency we here repeat the proof. Clearly, Ta 
is a distributional solution of the Klein-Gordon equation, 

(-□^ - a)Ta{x,y) = . 

Differentiating n times with respect to a gives (|D.8p . Next, we differentiate the identity 
in momentum space 

Taip) = Kp" - a) 0(V) 

with respect to pk to obtain 

A. 

Using that differentiation in momentum space corresponds to multiplication in position 
space and vice versa, we find 

Cfcr,(x,y) = 2AT«(x,y). 

Differentiating n times with respect to a gives ()D.9p . 

To derive (jD.lOp . we first combine (ID.OP with the product rule to obtain 



-^Ta{p) = 2pkTi^^\p). 
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On the other hand, we know from (jP.SP that 

a^T^ = -rw - ar(i) . 

Solving for Ta ^\ we obtain 

We finally differentiate this relation n + 1 times with respect to a, giving (|D.10p . □ 

Alternatively, this lemma could be proved by manipulating the series representa- 
tion ()4.13p . We also remark that in the case n = —1, the rule ()D.9p is consistent 
with our earlier definition (j4.16p . 

Using the relations (jP.Sp and ()D.9p , the mass expansion of the first order Feynman 
diagram ()D.3p was first performed in [12] (see \12\ Theorem 3.3], where the mass 
expansion is referred to as the "formal light-cone expansion" ) . More generally, for the 
advanced and retarded Green's function, we have the expansion 

(5«y5W)(x,y) 

which is proved exactly as in the case a = (see [I3l Lemma 2.1] or [Ml Lemma 2.5.2]). 
The residual argument (cf. [131 Section 3.1]) also generalizes immediately to the case a > 
0, making it possible to extend ()D.lip to the so-called residual fermionic projector (the 
non-residual part of the fermionic projector is precisely the non-causal high energy con- 
tribution, which can be analyzed with different methods as indicated in §8.3p . Applied 
to our problem, we obtain the expansion 

{Sa VTa+TaV 5J(x, y) = Y.-, /V " «')" (D"^) da Ti"+i)(x, y) , (D.12) 

n=0 

where for the line integrals we again used the short notation ()B.12p , and Sa is the sym- 
metric Green's function (|8.2ip . The mass expansion of ()D.3P is now readily obtained 
by applying the differential operators (i^+m) and simplifying the Dirac matrices using 
the rules (jP.OP and (jP.lOp . Multiplying by ^xl and taking the trace, a straightforward 
calculation using again (jP.lOp yields 

\Tv(ixLP{x,y)) X 2f;l l\a-a^Ti,{u-A\)T^^Xx.y) (D.13) 

n=0 ■ •'^ 

-'^T.^l r^^a - 1) (a - a^T r^r'^(x,y) (D.14) 

n=0 ■ •'^ 

-"^'E^ r(a-a2)"a(n"Ai + nM|)r^"+^)(x,y) (D.15) 

n=0 ^' 

+ 2m2f;l f\2a-l){a-a^T(u^d,Al)T^^^^\x,y) (D.16) 
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(this result was again obtained with the help of class_commute; see page [93]) . Integrat- 
ing the line integrals by parts, 

l\2a - 1) (a - a'r (□"5. A^,) = -1^ jja - a^r+' ^ (□"S.^^l) , 

the divergence terms can be rewritten to recover the chiral currents. In particular, in 
the case Al = Ar of a vector potential, one immediately verifies that (ID.13|) -( lD.16j) 
has the correct behavior under gauge transformations. Furthermore, one readily sees 
that the expansion (|Dl3]) - ([m6l) is compatible with (iD^jl - ljPTl) in the sense that the 
singularities on the light cone coincide. We now subtract ()D.13p -( )D.16P from (jP.Sp - 
()D.7p and solve for ^^fk{x,y). In order to compute fk{x,x), it suffices to take into 
account the constant counter term in (|4.14p . as can be done by the replacement 

r,(x,y) ^T-s + Ar(a) with N{a) := ^ alog\a\ , (D.17) 

and similarly for the a-derivatives. Moreover, in the line integrals we may set x = y. 
In order to keep the formulas simple, we also specialize to the situation where only the 
axial potential in (|6.22|) is present. We thus obtain 

f\x,x) = 2Y,- {a-a'T{u^Al{x))N^^\m')da 

n=0 •'^ 
°° 1 /<1 

_2j]^^-^y (a-Q2r+i(m^,4(x)) (iV("+i)(m2)-m2iV("+2)(m2))da. 

By linearity, it suffices to consider the case that is a plane wave of momentum q, 

Then the above sums are recognized as Taylor series, 

y^ATW (^2)^^(0) (^2^;^) 

^ n\ 

"■—^ fD 18) 

n=0 

where we introduced the abbreviation A = —(a — a^) (the second equation in (jP.lSP 
can be derived from the first by integration over A). We thus obtain 



f{x,x) 



da 

u=0 



2 Al{x)N[m^ + \^da 

-jj\a- a^) d''iAi{x) (atW _ ^^2^(1)^ (rn^ + ^) 
1 

= / ^a(^) ("i^ + A) log + A| da 

^ {a- Q^) d'^iAUx) log Im^ + X\ da , 

16tt'^ Jo 

where in the last step we substituted the explicit formula for N(a) in ()D.17p . In order 
to rewrite the last result in terms of the axial current, we use the identity Avl^ = 
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(a — a^)n^^ to conclude 



m 

167r3 



log 



m 



[a 



167r3 



Ja(^) / (a - a^) log 



a )q 



m 



da 



[a 



a )q 



da . 



Substituting this result into the light-cone expansion (|D.5p ~ ()D.7p evaluated at the 
origin and using that f^ia — a^) = 1/6, one sees that the term ^kf^{x,x) can be 
incorporated into the formulas of the light-cone expansion by the replacements 



(1) ,^(1) logM 
log(m^) 



n(l) 
[2] 



^(1) , 
^[2] + 



327r'^ 



+ 



+ 



6 



327r3 
1 

327r3 



(a — Q^) log 1 — (a — a^) — 2 



da 



log 



1 — (a — a^) 



(D.19) 



da 



To clarify the above construction, we point out that the radius of convergence of 



the Taylor series in (jP.lSP is |A| 



m 



Thus in the case |A| > m , these series do not 



converge absolutely, so that (|D.18P can be understood only on the level of formal Taylor 
series. For the reader who feels uncomfortable with formal power series, we remark that 
all formal expansions could be avoided by regularizing the distribution according 
to (j4.14p before performing the light-cone expansion, making a later resummation 
unnecessary. However, this method seems technically complicated and has not yet 
been carried out (see also the discussions in [131 Section 3.3] and after (18.250 ). In this 
paper, we will be content with the formal character of (jP.lSp . 
We are now ready to prove the main result of this appendix. 



Proof of Lemma \8.1\ We return to the situation with three generations and a general 
axial potential A^{z). As the axial potential is diagonal on the generation index, the 
auxiliary fermionic projector splits into the direct sum of three fermionic projectors, 
corresponding to the Dirac seas of masses mi, m2, and 777,3. Thus the partial trace (14.40 
reduces to a sum over the generation index. Decomposing into Fourier modes, 

for every A^{q) and for every generation we may apply the replacement rules (|D.19p . 
Rewriting the multiplication in momentum space by a convolution in position space 
gives the formulas (;5T^ - (f5T^ . 

In order to check the prefactors, it is convenient to verify whether the arguments of 
the logarithms can be combined to give dimensionless quantities. This is indeed the 
case with the expressions 

1 ^ 



logir 



1 



2 t2| 



3 ^ 

" =1 



3 ^ 

13=1 



logical + 



ra 



log (777^) 



m 



2YY 



log Irri^^^ 

/3=1 



explaining the prefactors in (I8.15P and (|8.16p relative to those in (18. 4p . 

We finally need to verify that the smooth contributions which were disregarded in 
the formalism of ^S.ll reallv enter the EL equations according to the simple replacement 
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rules (|D.19p . The subtle point is that the contraction rule in the continuum limit (j5.5p 
is not the same as the corresponding distributional identity ()D.4p . and this might give 
rise to additional terms which are not captured by ()D.19P . Fortunately, such additional 
terms do not appear, as the following consideration shows: To degree four on the light 
cone, the smooth contributions to P{x, y) enter the EL equations only if the smooth 
term is contracted with a factor ^ without generating a factor S,"^ (the contributions 
involving are of degree three on the light cone). Thus for the smooth contributions, 
the contraction rule (|5.5p is not applied, and therefore it could here be replaced by the 
simpler distributional identities ()D.4p and (jD.lOp . □ 

We finally carry out the a- integrals in ()D.19P in closed form and discuss the result. 
This result will not be used in this paper. But it is nevertheless worth stating, because 
it gives more explicit information on the structure of the non-causal correction terms. 

Lemma D.2. The functions /j^j defined by (|8.17p and (jS.lSp can be written as 



lim,[,] 



4m^ 



(D.20) 



where the functions g^p-^ (z) are defined in the upper half plane by 



5[0]( 



9[2]{z) 



3 + 5z ^ l + z 



2z' 



3z ' 2z ^z{z 



1 L 



log 1 - 2z + - ivr G(z - 1) 



11 



log 1 - 2Z + yjz{z - 1)) - iTT Q{Z - 1 



where the logarithm in the complex plane is as usual cut along the ray — 
is the Heaviside function, extended continuously to the upper half plane). 

Proof Writing the logarithm of the absolute value for any x E M as 



(and Q 



log |1 — x| = lim ( log (l — (x + i6)) 
we obtain the representation (ID.20p with 



ITT 



e{{x + i6) - 1 



m (z) = 6 



{a-a^)(^ log (l - 4(a - a^)z) + ivr G (4(a 



a'^)z 



1) da 



9l2]{z) = (\og{l - 4:{a- a^)z) +i-KQ{A{a 



a^)z 



I]] da. 



It remains to calculate these integrals for z in the upper half plane, thus avoiding the 
singularities on the real line. The term involving the Heaviside function is readily 
computed in closed form. Thus it remains to consider for ^ = 0, 1 the integrals 

^ {a-a^f\og(l-A{a-a^)z^da = ]^j^ log (l - xz) | (^^^ -^=L=| dx , 

where in the last step we transformed to the integration variable x := 4(a — a^). After 
computing the indefinite integral of the expression inside the curly brackets, we can 
integrate by parts. Then the logarithm in the integrand disappears, and the calculation 
of the integral becomes elementary. □ 
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4 m} 



Figure 6. The functions /^j (red, solid) and f^2] (blue, dashed). 

In Figure [6] the functions /jqj and f^^^] plotted. One sees that these functions 
attain their minimum if q"^ = 4m^, and for this value of q'^ the function has a cusp. The 
asymptotics for large is obtained by dropping the summand one in the argument 
of the logarithm in (IHTTl) and (IHIHI) . 



/fni(^) ~ 6 1^ (a-a^) log 



[a — a ) — 2 



m 



da 



3 Vm| 



1 

~ / log 





{a — a^) 



m 



da 



-2 + log 



/3 

revealing a logarithmic divergence as q^ 
have the asymptotics 

^2 

,4\ 



m 



iboo. For small momenta, the functions 



5m^ 



+ 0{q' 



6m^ 



+ 0{q'), 



describing a non-trivial low energy effect. 



Appendix E. The Weight Factors pjs 

In [18] the ansatz for the vacuum (j3.ip was generalized by introducing so-called 
weight factors pjs for the Dirac seas, 

P(x, y) = j2ppj + mp) 6{e - m}) G(-fe°) e'^'^^-y^ . 

This generalization turns out to be useful when considering an action principle for the 
masses of Dirac particles [23]; for a physical discussion see [181 Appendix A]. All the 
constructions of the present paper could immediately be extended to the setting with 
weight factors, as we now explain. 

The weight factors are introduced into the auxiliary fermionic projector of the vac- 
uum (14. ip by the replacement 

0=1 8=1 
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For our systems, the causality compatibility condition (see [TSl equation (A.l)]) does 
not cause problems, because all our potentials are either diagonal in the generation 
index, or else they can be described by a local axial transformation (see ^7.6p . in which 
case the causality compatibility condition is irrelevant. We conclude that the causal 
perturbation series as well as the light-cone expansion remain well-defined. The weight 
factors are taken into account simply by inserting them into the resulting formulas. 
More precisely, the number of generations is to be replaced by the sum of the weights, 

9 
13=1 

Moreover, the weights must be introduced into the partial traces by the replacements 

9 9 




ot,f3=l p factors y a,/3,7iv>7p-l=l 

or more generally using the rule 

9 

p=i 

(to avoid confusion, we note that due to the causality compatibility condition, the 
weight factors could just as well be inserted at the last instead of the first summation 
index). When considering the axial transformation, one must be careful to first multi- 
ply by the weight factors, then one performs the local axial transformation, and finally 
one takes the partial trace. 

The only place where the modifications caused by the introduction of the weight 
factors are not quite obvious is the construction of the local axial transformation 
in ^7.61 and Appendix [O More precisely, in ()7.3ip the factors ^ must be replaced by 
the matrix X^, where X = diag(yOi, . . . , pg) is the asymmetry matrix of the weight 
factors. As a consequence, the argument leading to (|7.33p no longer applies, making 
it necessary to arrange by the the additional condition 

9 

^=1 

that (|7.3ip vanishes. Then the subsequent analysis goes through if in (j7.35p and (|7.36p 
we apply (jE.ip and insert factors of into the sum. Likewise, in the non-perturbative 
treatment of Appendix [Cl we must insert the asymmetry matrix X into ()C.2p 

U^XU-^ = co^, 
and we must change the definition of powers of the mass matrix to 

YP = ^p^my (E.2) 

/3=1 

The subsequent analysis goes through with the following minor modifications. The 
first equatio n in ^ JO^ chang es to {l\VXVl) = {l\V-^XV-'^l) . Likewise, the first 
equations in (IC.13P and (|C.28P must be replaced by the conditions (n | Xn) = (m | Xn) 
and (m I Xn) = (n|Xm), respectively. These two conditions can be satisfied by the 
transformations (jC.lSP resp. ()C.30p . Then all results up to ()C.37p remain true. In 
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particular, Theorem IC.ll still holds for any vector field u which satisfies the condi- 
tion {u{x),u{x)) > —£ for some e > 0. In order to determine the maximal value of e, 
one must again find the maximum of the functional S. Here the weight factors make 
the analysis of the system ()C.40p a bit more complicated, and we leave the details to 
the reader. Proposition IC. 61 also remains true if we use ()E.2p and insert factors pp into 
the sum. 

After these modifications, all our formulas and results remain valid. It seems a 
promising strategy for the construction of realistic physical models to choose the 
fermion masses and the weight factors according to state stable vacuum configura- 
tions as exemplified in [23j . 

Acknowledgments: I would like to thank Andreas Grotz, Christian Hainzl, Joel Smoller 
and Alexander Strohmaier for helpful discussions and valuable comments on the man- 
uscript. 
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Notation Index 



(M, (.,.)) - Minkowski space, [5] 
P - fermionic projector, [S] 
P{x, y) ~ its integral kernel, [5l [8] 
1^,1^, .. . - slash, [8] 
^ - adjoint spinor, [5] 
< . I . > - inner product on wave func- 
tions, [5] 
Axy - closed chain, [H [27] 
I . I - spectral weight, [6] 
C - Lagrangian, [6] 
T - constraint, E] 
5 - action, [6] 

5 - number of generations, [8l [19] 
G - Heaviside function, [8] 

- Planck length, [8] 
E'p - Planck energy, [8] 
■^macro " macroscopic length scale, [8] 
£ - regularization length, [8] 
pe _ regularized fermionic projector, [8] 
^ - vector y — X, [3 [26l 
P^"^ - auxiliary fermionic projector. Illl 
[26] 

y - mass matrix, [11] 
m - mass parameter, [TT] 
S - perturbation operator, [TT] 
psea _ (describes filled Dirac seas, [121 [TJ] 
<.|.> - inner product, [12] 
^1, . . . , '^nf - particle states, [13] 
$1, . . . , - anti-particle states, [l3] 
rif - number of particles, [13] 
Ua - number of anti-particles, [13] 
0((y — - order on light cone, [l3] 
Ta ~ lower mass shell, [14] 
e - sign function, [TT] 
Ta^^ - infrared-regularized lower mass 
shell, M 

r^") - term of mass expansion, [TT] 

P^^ - non-causal low energy contribu- 
tion, [Ml [Ml [56] 

P^*^ - non-causal high energy contribu- 
tion, [12 [Ml [56] 

e - coupling constant, [El [65] 

■^L/R ~ chiral potentials, [151 [31] 

Xl/r ~ chiral projectors, [15] 

7^ - pseudoscalar matrix, [15] 



$ - scalar potential, [16] 
H - pseudoscalar potential, [T6] 
Tj^^ - ultraviolet regularized r^"^ [T9] 
\ " - notation for partial trace, [T9] 
^[p^^ - regularized ^,[20] 

Zj^^ ~ abbreviation for (^[jj^)^, [20] 

r|pj^ - ultraviolet regularized r("),[20] 

ri"^ - stands for T^^ or T^Jj^]!! 
deg - degree on light cone, [2Q] 
L - degree of simple fraction, [21] 
Creg - regularization parameter, [21] 
V - derivation on simple fractions, [2T] 
Sfj, - action involving Lagrange multi- 

pher, [23] 
( - adjoint of ^.[271 

A± - eigenvalues of closed chain in vac- 
uum, [27] 

F± - spectral projectors of closed chain 
in vacuum, [27] 

Q(x, y) - composite operator in EL equa- 
tions, [231 [281 [30] 

D - partial derivative of , [30] 

- vector potential, [3T] 

- axial potential, [3l] 

A^J^ ~ integrated chiral potentials, [32] 
^L/R ~ chiral phases, [32] 
A^^^ - eigenvalues of closed chain, [32] 
^L/i? _ gpg(.^j.g^j projectors of closed chain, 

m 

TZ - appears in EL equations to degree 

four, [33] 
o(l^l^) - order at the origin. [351 
ja - axial current, [35] 
c.c. - complex conjugate fraction, [55] 
•^L/iJ " chiral Dirac current, [36] 
Ja - axial Dirac current, [3U] 
X - denotes a contribution, [36] 
- generation mixing matrix, [TU] 
fj^j - Minkowski metric, [TO] 
U{x) - local axial transformation, [TTl 

[121[M] 

Ca,da - partial trace of g, [Ml 
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S[p] - smooth contribution to T^^^ , HU] 

/j^j - contribution to S[p]. H8H50] 
* - convolution, US] 

/j^j - Fourier transform of /j^j , HHl 11121 

Co 5^72 ~ regularization parameters, 09] 

Sa ~ Klein-Gordon Green's function, [50] 

Sm - Dirac Green's function, [51] 

tm ~ vacuum Dirac sea, [51] 

A - auxiliary potential, [53] 

r - radius |.^|, [M] 

tm - upper Dirac mass shell, [57] 

H - full Hamiltonian, [58] 

Hq - free Hamiltonian, [58] 

i? - perturbation of Hamiltonian, [58l 

-Bint - perturbation in interaction pic- 
ture, [5U] 

Fjj - field tensor, [US] 

M - bosonic mass, [63l [65] 

Js - scalar Dirac current, [66] 

Jp - pseudoscalar Dirac current, [66] 

Jb - bilinear Dirac current, [66] 

P± - half filled Dirac sea, [72] 

5(M) - Schwartz functions in momen- 
tum space, [72] 

n(a;, y) - nonlocal kernel, [ST] 

5(MxM) - Schwartz kernel in Minkowski 
space, [ST] 

f_[. - double null spinor frame, [91] 
- unperturbed closed chain, [21] 
- unperturbed spectral projector. [9T] 

-<.|.> — spin scalar product, [91] 

pec _ matrix elements in double null 
spinor frame, [U5] 

f§ - conjugate spinor frame, [92] 

JJ[p, (7 1 r] - short notation for line inte- 
grals, [Ml 

^I^IR ~ ^^"^^ fi^ld tensor, [95] 

jL/i? ~ chiral current, [95] 

y - axial transformation on C^. llOOl 

[, m, n - describe the axial transforma- 
tion, [Toni 

(.|.) - scalar product on C^. llOOl 
6 - functional describing axial trans- 
formation, 11041 
ri""^ - term of mass expansion for a > 0, 



f{x, y) - smooth contribution to xlP{x-, y), 
[TOSlfTTO] 

N{a) - smooth contribution to IllOl 
PfS - weight factors, 11131 
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Subject Index 



action, [U] 
adjoint spinor,[5] 
anti-particles, [T2] 

Cauchy problem, \UE\ 

causal perturbation expansion, [12] 

causality violation, [501 ED 

for spacelike distances, [57] 

in experiments, [52] 
closed chain, [5] 
continuum limit 

analysis in the, [21] 
contraction rule, [20l 
current 

axial, [35] 

axial Dirac, [36] 

bilinear Dirac, [66] 

chiral Dirac, [Ml [96] 

pseudoscalar Dirac, [66l 1981 

scalar Dirac, [661 EH] 
current term, [35] 

degree, [20] 

Dirac-Maxwell equations, [58] 

modified by convolution terms, [52] 
Dirac- Yang/ Mills equations, [57] 
Dirac- Yang/Mills-Higgs action, [63] 
Dyson series, [60] 

Euler-Lagrange equations, [23] 

in the continuum limit, [26] 
external field problem, [T2l [T6] 

fermionic projector, [5] 

auxiliary, [TT] 

homogeneous, [8] 

smooth contributions to, [35] 
Feynman diagram, [57] 

bosonic loop, [60] 

fermionic loop, [61] 

tree diagram, [T2] 

field 

axial conformal, [68] 
axial gravitational, [68] 
electromagnetic, [67] 
external, [121 [E] 



gravitational, [15] [67] 

Higgs, [63] 
Fock space, [El [T6] 
frequency 

constraint of negative, [TS] [77] 

gauge symmetry, [32] 

gauge transformation, [32l [35] 

axial, [63] 
generation mixing matrix, [30] 
generations 

at least three, [33 

several, [39] 

three, [33 
Green's function, [12] 

of the Klein-Gordon equation, [50] 

retarded, [59] 

Hadamard condition, [T7] 
Higgs mechanism, [35l [62] 
homogeneous perturbations, [75] 

inner factor ^, [20] 

inner product 

Minkowski_(., .),[5l[3Q] 
on spinors ^<I> or -<.|.;^. [5] [9T] 
on wave functions <.|.> , [5] [T2] 
scalar product (.|.) on O. llOOl 

integration-by-parts rule, [21] 

interaction picture, [59] 

Lagrangian, [6] 

light-cone expansion, [13l [90] [93] 
local axial transformation, 1^2] [88l [96l 
[98] 

logarithmic pole on the light cone, [37] 

mass cone, [9] 
mass matrix, [11] 
mass shell, [9] 
mass term, [35] 
minimizer, [7] 

modified Dirac-Maxwell equations 
in variational form, [23 

non-causal correction 

by convolution terms, [50] 
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of higher order, [56] 
non-causal low and high energy contri- 
butions, O SSI [56] 
noncommutative geometry, [TS] 

order 

at the origin, [35] 

on the light cone, [l3] 
outer factor ^, [20] 

partial trace, [Til [191 
particles, [l2] 
Planck energy, [8] 
Planck length, [8] 
potential 

auxiliary, [53] 

axial, [31] 

chiral, US [SH 

non-dynamical, [701 

nonlocal, [TQl [Ml [89] 

pseudoscalar, [HI [68] 

pseudoscalar differential, [TH [38l [96] 

scalar, [E] [Ml [68] 

vector, [3T] 

vector differential, [JT] 

quantum corrections, [571 
quasi-homogeneous ansatz, [81] 

regularization, [S] HT] [19] 

by cutoff, [65] 

by exponential factor, [M] 

special, [ST] [78] 
regularization parameter, [21] 

basic, [2T] 
regularized fermionic projector, [8] 
renormalization, [61] 

by counter terms, [T7] 

point splitting method, [T71 
resummation 

of current and mass terms, 11071 

of light-cone expansion, [TH 

simple fraction, [20] 
spectral weight, [6] 
spin scalar product. [911 
spinor frame 

double null, [HI 
spontaneous symmetry breaking, [35l [63l 



testing on null lines, [24l 

unitary in a compact region, [6] 
units, [51 [13 

weak evaluation on the light cone, [20] 
weight factor, 11131 
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